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ABSTRACT 



This thesis investigates the effects of atmospheric 
refraction on military ground-force tactical systems which 
utilize EH energy propagating through the troposphere. 
Anomalous propagation, mirages, and atmospheric turbulence 
are the three major atmospheric refractive conditions 
discussed. EM energy from visible light to VHF radio is 
considered . 

The refractive characteristics of the atmosphere and the 
conditions that lead to significant refractive effects are 
reviewed. The means of predicting atmospheric refraction 
and the usefulness of such predictions, along with the 
probability of occurrence of atmospheric refraction is 
also discussed. 

Specific examples of atmospheric refractive effects are 
provided in two forms : computer simulations and qualitative 

studies. A scenario is simulated showing the wide range 
of atmospheric refractive effects a typical ground force 
might anticipate in a desert environment. Finally, the 
thesis provides conclusions and recommendations for further 
studies . 
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X. INTRODUCTION 



In order to enhance and extend surveillance and 
communications capabilities, modern ground warfare has 
become increasingly dependent on the entire electromagnetic 
(EM) spectrum. Uses of electromagnetic energy include 
airborne and ground target detection, weapons system 
tracking and guidance, command and control, and electronic 
warfare. Utilization of this energy requires an under- 
standing of the phenomena that affect the propagation of 
EM waves as they travel through the atmosphere. After 
nearly 50 years of study, it is now widely accepted that 
dispersion, extinction, refraction, reflection, and 
diffraction all affect EM propagation significantly. 
Dispersion, extinction, and refraction are properties of 
the atmosphere, while reflection and diffraction are 
properties both of the atmosphere and of the terrain of 
the earth’s surface. For radio/radar frequency (RF) and 
electro-optical (EO) systems, the physical factors listed 
impact on system performance above and beyond the more 
obvious effects of precipitation, wind, and temperature 
extremes . 
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The primary region of Army interest is the troposphere- 
the lower 10 to 20 kilometers of the atmosphere.^ Within 
the troposphere, the main factor affecting RF propagation is 
refraction. EO propagation, while greatly affected by 
aerosols and absorption in the troposphere, is also affected 
by atmospheric refraction. Thus, in order to best employ 
his assets, it is important that the ground forces commander 
be able to forecast and predict refractive conditions that 
degrade or enhance RF and EO system performance. 

It is the purpose of this thesis to examine the effects 
of atmospheric refraction on ground forces' tactical systems. 
Specifically, it will address the environmental conditions 
leading to the three major forms of atmospheric refraction 
(anomalous propagation, mirage formation, and atmospheric 
turbulence), and discuss their effects on EM waves 
propagating through the atmosphere. It will also address 
prediction and probability of occurrence of these three 
forms of refraction and it will end with a discussion of 
their effects on existing ground systems, (either through 
simulation or by presenting empirical test results). 
Recommendations and direction for future research directions 
will be found in the conclusion of the thesis. 



L The specific height of the troposphere is a function of 
location and atmospheric conditions , on the average it is 
the thickest at the equator and narrowest at the poles. 
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II. ATMOSPHERIC REFRACTIVE EFFECTS 



In "this section "the basic theory of EM propagation will 
be introduced. Since the purpose of this thesis is to 
describe the effects of atmospheric refraction on ground 
forces, this discussion will be limited to the troposphere. 

A. EM WAVE PROPAGATION IN THE ATMOSPHERE 

1 . Wavelength Dependence of EM Propagation Effects 

The electromagnetic wavelengths of interest to U. S. 
ground forces extend from those only micrometers (10” meters) 
in length to those several tens of meters long. This 
encompasses the full spectrum, from HF to visible - an 
EM region which is described in Figure 1. 

Electromagnetic waves travel in a straight line in 
free space (vacuum) , but in the atmosphere they are bent by 
several forces. The specific atmospheric effect and layer 
of the atmosphere that dominate in causing this bending 
depend on the propagating frequency. At wavelengths longer 
than very high frequency (VHF), approximately 30 MHz, the 
ionosphere is the primary influence. Since this thesis is 
concerned with the lower atmosphere, HF frequencies will, 
therefore, not be specifically addressed. Above 30 MHz, 
refraction in the troposphere is the controlling influence, 
and it will be these frequencies, from VHF to visible, that 
this thesis will address. 
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Figure 1. Electromagnetic Spectrum 
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There is also a relationship between regions of the 
EH spectrum and the forms of refraction dominant in each 
region. The dominant effects at microwave frequencies are 
long range changes in refraction. At optical frequencies, 
absorption and small-scale refractive changes, in the form 
of mirage formations and atmospheric turbulence, are the 
dominant effects. In the millimeter wave (MMW) region the 
atmosphere is largely opaque due to molecular absorption by 
carbon dioxide and water molecules, except for a few regions 
(or "windows") of higher transmittance. Although little 
used until recently, much current research is focused on 
expanding the use of this region and in understanding the 
effects of atmospheric refraction on MMW propagation. 

2 . Index of Refraction (n) 

Electromagnetic waves travel through matter via 
absorption and re-emission of electromagnetic energy by the 
atomic and molecular constituents of the medium. Thus, the 
speed at which these waves travel is affected by the density 
of the medium - the denser the medium, the slower the speed 
of EM propagation through it. 

The "index of refraction (n)" is the ratio of the 
speed of electromagnetic waves in free space to their speed 
in the medium. Thus, 
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where c is the speed in vacuum (speed of light) and v is the 
speed in the medium. We note that v is dependent on the 
constituents of the medium - in this case, the atmosphere. 
Note also that n = 1 in vacuum, and, therefore, n is greater 
than or equal to 1 for all media. For air, 



n AIR = 1 * 0003 



In the troposphere, n^^ is a function of temperature (T), 
water vapor content (e), and pressure (P), all of which are 
related to the density of the atmosphere and all of which 
vary systematically with climate and elevation. In general, 
the index increases with increasing density. 

3 . Snell’s Law 

The deflection of an electromagnetic wave across a 
discontinuity or boundary separating regions of differing 
values of refractive index is shown in Figure 2, and is 
described by Snell’s law: 



Sin 9^/Sin 0 2 = n 2 /n^ 

Note in Figure 2 that the ray is bending (retarded) towards 
the higher n value, so that refractive bending is toward 
media of higher density. 

Since the density of air generally decreases with 
altitude, the value of n will similarly decrease. Thus, in 
a "standard atmosphere" (a term used to describe average 
meteorological conditions), EM waves do not travel in a 
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straight line but bend slightly downward towards the earth's 
surface, the region of higher n. A list of the factors 
which, according to the National Oceanic and Atmospheric 
Administration (NOAA) [Ref. 1: p. 20], constitute "standard" 
conditions is provided in Table 1, while Figure 3 is a 
graphic illustration of the effect of index of refraction 
variation, with altitude, on EM propagation. When 
discussing changes of the index of refraction in the 
atmosphere, it is important to realize that there are no 
abrupt boundaries, as is implied in Figure 3. Changes in n 
generally occur over relatively large distances, when 
compared to the wavelengths affected. It is these gradual 
changes that create the smooth bending effect shown in 
Figure 4, as opposed to the sharp angles shown in Figures 2 
and 3. Figure 4 also shows another of the effects of 
refraction in the troposphere, i.e., extension of the 
horizon. In essence, an EM ray directed above the geometric 
horizon is refracted and bent slightly around the surface of 
the earth, striking the surface at a distance greater than 
the geometric horizon. For example, J. Lake [Ref. 2: p. 22] 
states that, for microwave rays propagating in a standard 
atmosphere, the horizon is extended beyond the geometric 
horizon by approximately 30%. 

Should the decrease in n, with altitude, be 
sufficiently great, a non-standard atmosphere will exist, 
where EM rays will bend downward more strongly than in a 
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Figure 2. Snell’s Law. 



TABLE 1 

METEOROLOGICAL CONDITIONS FOR THE STANDARD ATMOSPHERE 

At altitude z = 0 (sea level) 

Pressure(P) = 1013.25 mb 

Temperature (T) = 15°C (288. 15K) 

2 

Acceleration (g) = 9.8 m/s 
Density (p) = 1.225 Kg/m 3 
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Refracted 

Path 



Earth 



Ray 



Figure 3. Path of microwave through the troposphere. 

The radar transmitter is located on the 
earth's surface and transmitting along path 
A. The atmosphere is assumed to consist of 
concentric shells each with a constant 
refractive index. 




Figure 4. Extension of the radar horizon in the 

presence of refraction. A is the radar hori- 
zon in the absence of refraction (the 
geometric horizon) and B is the radar 
horizon in the presence of refraction. 
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standard atmosphere. This can lead to the most dramatic 
form of atmospheric refraction, called "trapping”, where a 
horizontally launched ray is bent sufficiently so that it 
remains parallel to the earth's surface, extending the 
horizon to great distances. If we consider the atmosphere 
to consist of layers of varying refractive index, n, there 
is a critical angle of incidence needed to refract an EH 
wave completely back from a layer of lower n value to a 
layer of higher n value. At all angles of incidence smaller 
than the critical angle, EM rays are "trapped" by the layer 
of higher n. For all angles larger than the critical angle, 
the rays will be bent, but not so sufficiently as to be 
trapped. Thus, a region just above the trapped rays is 
formed, where fewer EM rays exist than if standard atmos- 
pheric conditions prevailed. This region is referred to 
as the radio or radar hole, the basic concept of which is 
shown in Figure 5. According to John Beach [Ref. 3: p. 82], 
the use of the word "hole" is not entirely accurate, since, 
under normal circumstances, there is no distinct critical 
angle, and leakage of EM energy into this region can occur 
by means other than direct propagation of the signal. There 
is, however, significantly less energy present in these 
radar/radio holes than would be present if no variation in n 
existed. 

The critical angle for trapping, 0 , as derived from 

Snell's law, is that angle, 9, which causes 02=90. Thus, 
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Figure 5. Radar/Radio Hole. 
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(2.3) 



Sin 0 = n„/n. 

c 2 1 

Since n only varies from about 1.00025 to 1.004 in the 
troposphere , n 2 /n^ = . 9998 and 9^ ~ 89°. Hence, as shown in 
Figure 6, EH rays must be almost parallel to the layer in 
order to be trapped. 

The index of refraction of the atmosphere does not 
always decrease with altitude. Should n increase with 
altitude, EH rays will be bent upward instead of downward. 

Hore will be said about this latter case, in discussions on 
subrefraction . 

B. REFRACTIVITY 

1 . Refractivity (N) 

The tropospheric index of refraction, n, is a 
characteristic of lower atmospheric properties. Specifically, 
at radio and microwave frequencies, n is a function of 
temperature, humidity, and pressure, and is related to these 
parameters by: 

(n-l)xlO 6 = (77.6 P/T - 5.6 e/T + 3.75xl0 5 e/T 2 ) = N, 

(2.4) 

where P is the barometric pressure in millibars, e is the 
partial pressure of water vapor in millibars, T is the 
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998 = 89 ° 



Figure 6. 



Critical Angle 0 c * 
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absolute temperature in Kelvin, and N is the refractivity , 
a parameter defined for the sake of convenience. 

The second term in Equation 2.4 is often neglected, 
since its contribution to the calculation is very small 
compared with that of the other two terms, and the 
coefficient on the third term is usually modified so that 
the actual equation is generally: 

N = (77.6 P/T + 3 . 73xl0 5 e/T 2 ) (2.4A) 

Beach [Ref. 3: p. 75] states that since the value of 
n at the earth's surface is usually between 1.00025 and 
1.0004, the value of N generally lies between 250 and 400. 

At optical frequencies the calculation of refrac- 
tivity is modified to reflect the fact that water vapor has 
a negligible effect. Thus, for optical frequencies: 

N = 77.6P/T (2.5) 



The significance of refractivity is not in the 

individual value of N for a given altitude, but, rather, in 

the gradient of N, (dN/dZ), where Z is height. This is due 

to the fact that the electromagnetic wave path curvature 

depends, inversely, on the gradient of N. For microwave 

frequencies and a standard atmosphere, N decreases at a 

2 

fairly steady rate of approximately 4x10 N per meter of 
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altitude. It is this characteristic that makes it convenient, 
in computations involving a standard atmosphere, to replace 
the actual earth's radius with a 4/3 earth radius. Using 
this geometry EM waves can be shown to propagate in straight 
lines rather than in curved or refracted lines, as shown in 
Figure 7. The value of using the 4/3 earth radius lies in its 
simplification of distance computations above the earth's 
surface, and in ray tracing. This geometry has been widely 
used in the modeling of radio communication and radar propa- 
gation. However, it must be remembered that the 4/3 earth 
representation is only an approximation of conditions that 
vary with meteorological conditions. The actual value may 
vary from 6/5 to 4/3 earth radius. 

A vertical refractivity gradient causes EM rays to 
be bent with a radius of curvature given by: 

r =-l0 6 / (dN/dZ ) , (2.6) 

where the curvature is in a vertical plane. 

Horizontal N gradients are normally not significant 
in EM propagation and are, therefore, usually assumed to be 
zero. In. some circumstances, however, bodies of water either 
partially or completely surrounded by land may create 
horizontal N gradients which will have some effect on EM 
transmissions through the atmosphere above the 
land-water boundary region. 
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Figure 7a. Bending of antenna beam due to refraction 
by the troposphere. 




Figure 7b. Shape of antenna beam in 4/3 earth 
representation. 
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2 • Anomalous Propagation, Mirage , £ Optical Turbulence 

Refractivity effects are not the same across the EM 
spectrum. As previously shown, the calculation of N for the 
radio and radar portion of the spectrum is different than 
"that for the optical portion (where the effect of humidity 
is ignored). Since the amount of bending is a function of 
the gradient of N, it is to be expected that the amount of 
bending is different for optical and microwave rays. This 
can be seen in Figure 8 . Note the difference in the bending 
of the optical grazing ray and the radio grazing ray 
vis-a-vis each other and the geometric horizon. 

Regions of strong refractivity gradients that are 
much larger in size than EM wavelengths cause anomalous 
propagation. Smaller scale regions with large N gradients 
also exist, in the form of mirages and atmospheric turbu- 
lence. Certain forms of atmospheric refraction, when 
located within the first several meters of the surface of 
the earth, and resulting from rapid changes in the 
temperature, can lead to mirage conditions which have a 
significant effect on optical transmission. Mirage 
formation is a form of anomalous propagation, but, for the 
sake of simplicity, the two will be discussed separately. 
They are best differentiated by relative size. Turbulence 
in the atmosphere creates fluctuations in the air index of 
refraction, resulting from inhomogeneities that vary in size 
from those of almost microscopic dimensions to those 
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Geometric 
(straight) ray 



Grazing ray 
( optical ) 

Grazing ray 
(radio) 



A Geometric horizon 
B Optical horizon 
C Radio horizon 



Figure 8. Comparison of effect of refraction on 
optical ray and radio ray. 
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hundreds of meters in diameter. The wavelengths in the IR 
and visible range are most affected by turbulence, due to 
<-he small size of the refractive inhomogeneities relative to 
the wavelengths of radio and radar signals. Each of these 
three conditions will be individually addressed in the 
following sections, beginning with a discussion on anomalous 
propagation. 

C. ANOMALOUS PROPAGATION 

Anomalous propagation is a misnomer, referring to any 
propagation other than that in a standard atmosphere. 
Anomalous propagation is generally considered to be divided 
into three forms: superrefractive , subrefractive , and 

ducting. Although the name implies rare or infrequent, the 
refractive conditions listed above are far from infrequent. 

C. G. Purves [Ref. 4: pp . 13-14] states that in parts of the 
world, and during certain seasons ducting or trapping occurs 
more than 75% of the time. Although the frequency of 
occurrence of various propagating conditions will be 
addressed in more detail, let it suffice to say that 
anomalous propagation is far from uncommon. 

Another point to address is the fact that superrefrac- 
tion, subrefraction, and ducting occur on a macroscopic 
scale. They are characteristic of relatively large 
seasonal and/or atmospheric conditions, some of which will 
be discussed at the end of this section. This characteristic 
roughly separates anomalous propagation from the conditions 
leading to mirages. 
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1 • Superrefraction, Subrefraction, and Ducting Defined 



Equation 2.4 shows that if water vapor pressure 
decreases rapidly with height or if the temperature increases 
with height, N will decrease with height at a more rapid rate 
than under standard atmospheric conditions. Thus, the value 
of dN/dZ becomes more negative. For small gradients a 
superrefractive condition exists, while when EM rays are 
curved parallel to the earth's surface trapping occurs and 
the phenomenon of ducting exists. Should the index of 
refraction increase, instead of decrease with altitude, a 
subref ractive condition exists. The meteorological conditions 
leading to these forms of anomalous propagation will be 
discussed in detail in Section II. C. 4. of the thesis. 

Superrefraction implies a gradient greater than that 
for a standard atmosphere, but not strong enough to form a 
duct. EM waves propagating in this region will be bent 
downward more than in normal (standard atmosphere) conditions 
but not sufficiently for trapping to occur. Superrefractive 
conditions will result in propagation distances beyond 
those expected in a standard atmosphere. 

Trapping, or ducting, occurs when dn/dZ is less than 
or equal to -157. Signals traveling in a duct propagate in 
a manner similar to propagation in a "leaky” waveguide. 

Ducting is categorized into three subgroups: evaporative 

ducts, surface ducts, and elevated ducts. Evaporative ducts, 
most commonly created by evaporation off of the surface of 
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water, are present only to a maximum height of about thirty 
meters and can be much thinner. These ducts are present over 
the surface of the oceans of the world to some degree most of 
the time. 

J. Beach [Ref. 3: pp. 78-80] states that surface 
ducts, which are usually stronger and generally extend higher 
in altitude than evaporative ducts , vary in thickness but are 
almost always less than one kilometer and usually are 300 
meters or less in thickness. Surface ducts are characterized 
by the lower boundary being the earth’s surface. Surface 
ducts, along with elevated ducts, occur much more frequently 
over land than do evaporation ducts. 

Beach [Ref. 3: p. 80] in his discussion of elevated 
ducts states that elevated ducts begin at altitudes of from 
close to zero meters, absolute altitude, to 6 kilometers but 
generally no higher than 3 kilometers. The thickness of 
these ducts is similar to the surface duct. 

When the index of refraction increases with altitude, 
a subrefr active layer is created. This type of layer is 
normally created by heating of the surface and thus exists 
in only a thin layer close to the earth’s surface as shown 
in Figure 9. EM waves entering such a layer from above will 
be bent upward, as will rays traveling within the 
subrefractive region. Subrefraction is responsible for some 
of the more common forms of optical mirage (see Section D). 
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Surveillance 

aircraft 




Figure 9. Description of Subrefraction. A 

subrefractive layer bends EM energy 
upward, away from the surface of the 
earth. 
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2 . Means of Describing Refraction 

As shown previously, the refractivity , N, is a 
parameter that normally decreases with height and is 
therefore not very useful in depicting refraction 
graphically. Two indices that are more commonly used for 
this purpose are B and M units. 

B is defined by K. Davidson [Ref. 5: pp. 4-2 to 4-4] 
by adding 40 N units per kilometer to all N values 

B = N z + 40Z , (2.7) 

where N z is the value of refractivity at height Z in 
kilometers. Thus, the B gradient is zero for standard 
atmospheric conditions. 

Davidson goes on to describe the "modified 
refractivity" M, as being defined by adding 157 N units 
per kilometer to N, thus 

M = N + 157Z (2.8) 

z 

The M gradient is zero when EM ray curvature equals the 
earth's curvature. Since, in a standard atmosphere, M 
increases with height, layers in which trapping can occur 
will be readily apparent by the negative M gradient. Figure 
10 shows the relationship between N, B, and M units. 
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Figure 11 shows the Navy Integrated Refractive Effects 
Prediction System (IREPS) classification of refractive 
conditions. See also Appendix A for more information on 
IREPS. 

3 . Duct Characteristics 

Before describing the causes of ducting it is 
necessary to explain some of the characteristics of ducts. 

To begin, an important distinction must be made between the 
trapping layer and the duct. As shown in Figure 12, the 
trapping layer is the region where dM/dZ is less than zero. 
Consequently the EM rays are bent downward with a curvature 
greater than that of the surface of the earth, causing them 
to be trapped between the layer and the bottom of the duct. 

A duct is the region below the top of the trapping layer 
where the EM rays are confined or trapped and includes the 
trapping layer itself. The EM rays are trapped because they 
are partially confined between the top of the trapping layer 
(duct) and the bottom of the duct. The trapping layer is 
always at the top of the duct. The duct bottom can be 
located on the earth's surface (surface duct), or can be 
elevated above the ground (elevated duct). 

The top of the duct is where the slope of dM/dZ 
changes from negative to positive. At this height M is no 
longer decreasing with height but begins to increase in 
value, as shown in Figure 13, for surface and elevated ducts. 
The bottom of the duct is located by dropping a line 
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Height In Kilometers 




Figure 10. Profiles for N, B, and M units for a 
standard atmosphere. 




Figure 11. IREPS Classification. 
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Figure 12. Duct Characteristics. 
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Figure 13a. Elevated Duct. 




Figure 13b. Surface Based Duct With' Elevated 
Trapping Layer . 
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vertically from the location of the top of the duct, in 
Figure 13, to the earth* s surface (for a surface duct) or 
to the intersection with the M curve (for an elevated duct). 
The optimum coupling height into the duct is that point where 
dM/dZ goes negative. This is the point where a radar/radio 
transmitter will achieve the greatest extension of range due 
to the trapping effect of the duct. The strength of the 
duct is determined by the most negative value that dM/dZ 
attains and by the thickness of the negative layer. The 
more negative the value the stronger the duct. Note in 
Figure 12 that a surface duct can be caused by a surface 
trapping layer or an elevated trapping layer. The bottom of 
the duct is the location where the downward bent EM rays are 
turned back upward again. It is apparent from this 
discussion that ducts can vary considerably in thickness 
and, depending on the slope of the trapping layer, they can 
also vary considerably in intensity (strength of the duct). 

As will be described later, not all EM frequencies 
are trapped within the duct and only those EM rays entering 
the duct or traveling through the duct at a very shallow 
angle will be trapped. Note that for the EM ray to be most 
effectively trapped the transmitter must be within the duct 
although EM signals can be coupled into the duct and trapped 
within it from an antenna located near but outside the duct. 
In addition, note that the trapped EM ray is either 
refracted back upward after reaching the bottom of an 
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elevated duct or reflected upward off the surface, as in the 
case of a surface or evaporation duct. In both of these 
cases at least part of the energy of the EM ray is lost due 
to surface absorption and diffuse scattering (surface and 
evaporation ducts) or lost through leakage through the 
bottom of the duct (elevated duct). Further, part of the EM 
ray is lost by leakage through the top of all ducts. 

Diffuse scattering and absorption are greater over 
land than over water. M. L. Meeks [Ref. 6: pp. 13-26] states 
that these two surface characteristics result in the surface 
reflection coefficient, the value of which varies as a 
function of frequency, angle of incidence of the reflected 
ray to the surface, and the terrain surface itself. The 
lower the angle of incidence, the smoother the terrain, the 
more moist the surface, and the lower the frequency, the 
greater the amount of reflection of the EM ray. The 
reflection coefficient for water varies from approximately 
.85 to .99 and is largely a function of the sea state. The 
reflection coefficient over land is much more variable, 
ranging from as low as .15 to higher than .9. The 
reflective nroperties of various land surfaces and types of 
vegetation are still not completely understood. 

The amount of EM energy lost as a signal travels 
through a duct is a complicated function of frequency, duct 
characteristics , and in the case of surface and evaporation 
ducts, the terrain’s reflective characteristics. Due to 
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these losses and those resulting from atmospheric scattering 
and absorption it can be seen that there is a limit to just 
how far an EM signal can be propagated in the atmosphere. 

It is apparent from this discussion that ducting of signals 
is less effective in increasing EM signal range over land 
than over water. 

As mentioned earlier, not all frequencies are trapped. 
It is the thickness of the duct itself that is the determining 
characteristic as to which frequencies are trapped. Duct 
thickness is therefore used to calculate the lower limit of 
those frequencies that are strongly trapped, as shown in the 
following formula from Davidson [Ref. 5: pp . 4-15]: 

f = (3.6033 x l0 1:L )d" (3/2) Hz (2.9) 

m 

Frequencies below f are still affected by the duct to a 
- m 

somewhat lesser extent. Table 2 shows values of f for 
various duct thicknesses. 

It is important at this point to reiterate that not 
only does the thickness of the duct limit the frequencies 
trapped within a duct, but so too does the angle of 
incidence of the EM rays entering the trapping layer. As 
mentioned earlier, only EM rays incident at very shallow 
angles relative to the horizon will be refracted enough to 
be trapped. In actuality, angles above 1 degree -from the 
earth’s horizon are not trapped and, according to Skolnik 
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[Ref. 7: p. 450], refraction is troublesome only at angles 
of elevation of 5 degrees or less from the horizon in most 
microwave applications. Figure 5 provides a simplified 
example of the effect of angle of incidence on refraction, 
trapping, and creation of a radar/radio hole. 

4 . Causes of Superrefraction, Ducting, £ Subrefraction 

Superrefractive , subrefractive and ducting conditions 
are created by variations in water vapor and/or temperature 
in the atmosphere, brought on by a multitude of climatic and 
topographic conditions. 

Evaporation ducts are created above the surface of 
large bodies of water by the rapid decrease in humidity above 
the water's surface. The humidity of the air immediately 
above the water is at the saturation point, but decreases 
quickly with height. Within the first thirty meters, an 
ambient level of humidity is reached, dependent on the 
climatic conditions prevailing. This rapid decrease in 
humidity causes M to decrease with altitude to a point where 
M will reach a minimum (the top of the duct) and then 
increase with altitude. Meeks [Ref. 6: p. 9] states that 
over land a similar situation is created after a rainfall, 
when the water evaporates from the ground. Snow, and 
particularly melting snow, can also lead to the creation of 
an evaporation duct over land. Overland evaporation ducts, 
however, are generally of short duration. Evaporation ducts 
will also occur over the continental land masses above 
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larger bodies of water. For the Navy, the occurrence and 
effects of evaporative ducts are very important and largely 
responsible for the development of the IREPS program which 
models the effects of superrefraction and ducting on EM 
communications and on surveillance systems. 

Superrefractive conditions and surface and elevated 
ducts are created by the same climatic conditions. Which 
type of propagation condition occurs depends primarily on 
where the transition between two differing air masses forms 
and on the strength of the resulting refractivity gradient. 
Temperature inversions are the primary cause of superrefraction 
and ducting. As would be expected, superrefractive conditions 
are far more common than ducting. 

A temperature inversion is a region of the atmosphere 
where there is an increase in temperature with height. As 
shown in Equation 2.4, a temperature inversion must be very 
strong to cause superrefraction by itself. If, however, a 
water vapor gradient is combined with the temperature 
inversion, superrefraction and ducting will more easily 
result . Thus , it is no surprise to find that superrefraction 
and ducting are more common over water (particularly oceans) 
than over land. 

Ducting is characterized by the upper air being much 
warmer and drier than the underlying air. There are several 
meteorological conditions that can lead to this situation. 

Over land, the most common ducting is caused by radiation of 
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heat from the earth* s surface on clear calm nights. This 
surface cooling creates a temperature inversion and a sharp 
decrease in moisture with height which will often lead to a 
surface cog in summer. Refractivity gradients are made even 
stronger if the ground is moist, a situation most common in 
summer and in winter over snow. Meeks [Ref. 6: d. 9] 
elaborates on the fact that in desert regions the temperature 
inversion created at night can be significant enough, in and 
of itself , to create superrefractive conditions and even 
ducting. Ducting is therefore most common over land at night, 
and will usually weaken or disappear, along with the weakening 
and disappearance of the temperature inversion, during the day. 

Subsidence, the sinking of air from above, also is a 
common cause of ducting. As the air sinks it is warmed by 
compression capping the cooler air below and a strong 
inversion is formed. If the lower level air has a high 
moisture content, a strong humidity gradient will also be 
formed. Since subsidence is most often associated with high 
pressure areas , this type of ducting is generally a fair 
weather phenomena, characterized by light winds, and clear 
or low stratus skies. The presence of fair weather low 
stratus clouds that show a marked discontinuity at the top 
of the cloud mass will often indicate the presence of an 
inversion and a good possibility that superrefraction or 
ducting is present. Figure 14 shows temperature inversions 
created by radiation of heat and by subsidence. 



Height Height 







Figure 14. Types of Inversion: (A) Subsidence reaching 

the surface; (B) Subsidence above a mixed 
surface layer; (C) Ground Inversion caused 
by cooling of surface air; (D) Inversion 
caused by cooling and mixing of surface 
layer. 
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In temperate regions in the low and middle latitudes, 
ducting occurs most often in summer and does not usually 
occur during bad weather. Storms are characterized by the 
lower atmosphere being mixed by precipitation and winds, 
resulting in a more homogenous atmosphere. In arid regions 
humidity will have a lesser impact on the formation of 
anomalous propagation conditions but in temperate regions 
the air temperature may be high enough, particularly in 
summer, for the water vapor component of the atmosphere to 
be important. Seasonal variations in refractive conditions 
are strongest over large land masses and generally weakest 
over large bodies of water (oceans). 

As described earlier, ducts over land generally are 

weakest or break up entirely during the warmest part of the 

2 

day. This diurnal characteristic of the land does not, 
however, occur significantly over large bodies of water. 

W. Thompson [Ref. 8: p. 3] additionally notes that in 
temperate regions , diurnal variations are more pronounced 
during summer months and less so in the winter. C. Samson 
[Ref. 9: p. 6] states that due to the rapidly changing 
thermal structure of the air near the earth’s surface at 



Z The changes in the thermal structure are due to the 
effects on the atmosphere of the heating and cooling of the 
earth's surface. 
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sunset and sunrise, the most extreme changes in refractivity 
will occur at these times. 

Skolnik [Ref. 7: p. 452] states that the main 
exception to the fair weather rule for ducting is the 
creation of ground ducts by thunderstorm downdrafts. These 
will usually last on the order of 1 hour or less and are 
small in spatial extent. 

Coastal regions develop their own unique ducting 
conditions. As described in Davidson's notes [Ref. 5: 
pp. 7-30 to 7-31] the presence of ducting in coastal regions 
and straits is quite common and is due to the jutaposition of 
marine and continental air masses at the surface. Offshore 
flow of warm dry air would lead to warm, dry, continental 
air overlying the cooler and more moist marine air. Offshore 
flow can be caused.by: 1) local land to sea breezes induced 

by the differential heating of land and water areas and 
sea-breeze recirculation, and by 2) larger scale circulations 
in which an offshore flow occurs with a semi-stationary high 
pressure area located over land. Figure 15 shows both of 
these conditions. 

In the latter situation, extremely strong ducts can 
be formed, some of which will extend well inland from the 
coast. The Santa Ana winds and monsoon circulations are 
characteristic of this type of offshore flow. 

Subrefraction is the rarest of the anomalous 
propagation categories. The only regions where n increases 
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Figure 15a. Illustration of elevated layer occurrence 

due to local differential land-sea heating. 




Figure 15b. Illustration of elevated layer due to 
synoptic scale driven offshore flow. 
Note adiabatic heating with downs lope 
traj ectorv . 
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with height with any regularity are the deserts of the world 
and particularly those that are near large bodies of water. 

The sun heats the desert surface causing it to be very hot 
and dry, the temperature to decrease, and the water vapor 
to increase, with height. This subrefractive condition can 
be made even stronger by an onshore breeze which is 
characteristic of a high pressure area located at sea. The 
onshore breeze carries cool moist air that rises over the 
land . 

Subrefractive conditions tend to follow a diurnal 
cycle and are strongest in the afternoon and usually disappear 
at night. In semi-arid and desert regions it is not uncom- 
mon to see strong subrefractive conditions forming in the 
afternoon, due to solar heating of the surface, and super- 
refractive conditions develop at night, due to radiation 
of heat from the surface. 

5 . Consequences of Anomalous Propagation 

As described earlier, the effects of anomalous 
propagation vary with the wavelength and the gradient of the 
refractive index. Increases in the refractive gradient 
(due to superrefraction or ducting) result in an extension 
of the radar horizon, causing an increase in radar coverage 
beyond that expected for a standard atmosphere, as shown 
in Figure 4. Refraction can also lead to errors in the 
measurement of range and elevation angle by microwave 
systems. An elevation angle error is created by the bending 
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of the EM waves away from a straight path as they propagate 
to the target and back as shown in Figure 16 . Range error is 
created by the increase in the effective path length due to 
refraction of the EM rays from the radar transmitter to the 
target and back to the radar receiver. Figure 17 shows 
range error as a function of the apparent elevation angle. 

An example of the type of range error that can be created by 
atmospheric refraction is that experienced by the FPS-16 air 
search and tracking radar at Mount Lemmon, Arizona. For 
this radar, range errors are generally between 30 and 40 
yards over a 150,000 to 200,000 yard distance. In many 
instances this type of error is of little consequence, but, 
for accurate tracking of a target's location some 
compensation for the range error must be made. 

Communications signals are affected by refraction in 
a manner similar to radar signals , leading to extended ranges 
in "line of sight" communications. However, the terrain 
over which the signals propagate will often have a greater 
effect on communications range than will refraction for 
ground based systems. 

The effect of refraction on Direction Finding (DF) 
systems can, in certain circumstances, be significant. 
Harrington and Callaghan [Ref. 10: p. 3] state that changes 
in the refractive index, particularly if occurring 
horizontally, will lead to the same type of error in 
direction finding as the elevation angle error in a radar, 
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Figure 16. Angular Error Caused by Refraction. 
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Range Error (feet) 




Figure 17. Range error as a function of apparent 

elevation angle, where apparent elevation 
angle is defined as the angle from the 
horizontal that the radar is aimed in order 
to locate the target (as opposed to the 
angle from the horizontal that the target 
is actually located). CRPL reference 
refractivity atmosphere (1958) is 
assumed. Note that the closer the apparent 
elevation angle is to the horizontal, the 
greater the range error. 
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as shown in Figure 18. An example of a situation that might 
result in a horizontal gradient is a signal transmitted 
across a body of water where the air is more humid than over 
land. Similarly, a signal traversing a wide valley, where 
the air is much drier than the surroundings, can result in 
an error in location of a transmitter. 

Multiple modes on an EM signal propagating in a duct 
can result in variations in the strength of the EM signal at 
various points along the path due to interference. This 
spatial fading can be caused by multipath, or, as in the 
above example, by multimode interference, both of which can 
result in a condition where the vector sum of several 
components of a signal are not in phase. The variation in 
phase is a result of the difference in path length induced 
by refraction and/or path route. According to Skolnik 
[Ref. 7: p. 454] fading on the order of 20 dB can occur for 
radar systems experiencing multimode propagation. Fading 
can be greatly enhanced by the presence of ground and/or sky 
waves combining with the refracted signals (multipath) . 

Beach [Ref. 3: p. 75] states that temporal fading will 
frequently occur near the radio horizon due to changes in 
the index of refraction leading to large changes in signal 
strength in the horizontal EM signal path. 

In the presence of ducting, the effects previously 
described become more acute . The angle error in a ground 
based radar will become greater with a larger refractive 
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Denser Air Pocket 




Figure 18 . Signal Refraction because of change of 
humidity over a lake and its effect on 
D.F. accuracy. 
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gradient and can often be more serious than the range errors 
induced. For those VHF, UHF , and microwave signals propa- 
gating nearly horizontal to the earth's surface, trapping 
can also occur, particularly if the transmitting antenna is 
located within the duct. When trapping occurs, the range of 
the EM signal within the duct is greatly extended beyond 
that expeeted in a standard atmosphere as shown in Figure 
19. Thus, for surface based radar signals the coverage can 
be extended well beyond the expected radar horizon and, for 
communications signals, extraordinary ranges have been 
caused by the presence of surface based ducts. Examples of 
VHF signals propagating 2,000 miles without the use of 
satellites are well known. Likewise, electronic support 
measures (ESM) detection of these signals is greatly enhanced 
by the extension of the range of the propagating signal. For 
ground based DF systems, trapped signals can be a mixed 
blessing. Although the range at which a signal can be 
detected is extended, the accuracy of the direction finding 
can be degraded by the very nature of the ducted transmission. 

The existence of radar/radio holes is indicated by 
regions in the atmosphere where EM signal strength is weaker 
in the presence of a duct than it would be in free space. 

The variables that dictate the existence and characteristics 
of a radar or radio hole are: the radar or radio transmitter 

height relative to the duct, the radar target or ESM 
receiver height relative to the duct, the operating 
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frequency of "the transmitter, the duct height, duct strength, 
and the optimum coupling height of the duct. It must be 
realized that these radar/radio holes are not entirely free 
of EH waves . A target can be detected by signal leakage or 
signals propagating in these regions by other means but the 
probability of detection is much less than for free space. 

Surface based radar coverage of the area just above a 
surface duct can be extremely difficult to achieve due to 
the existence of a radar hole. In essence, targets flying 
within the duct can be detected at up to 3 to 4 times the 
normal range whereas a target flying just above the duct may 
not be detected until very close to the radar. 

EW and surveillance aircraft will similarly be 
affected by the presence of ducts, particularly elevated 
ducts. As an example of such effects, John Beach [Ref. 3: 
p. 82] uses studies of airborne early warning aircraft 
flying in and around atmospheric ducts. If the aircraft is 
located within the duct, its detection range within the duct 
is greatly extended as shown in Figure 20. Below the duct, 
its detection is largely unaffected and above the duct there 
will be a hole in the coverage. If the aircraft is just at 
the top of the duct, the coverage within the duct is not 
generally increased but a hole still exists above the duct. 

As the aircraft rises above the duct, the radar hole at the 
top of the duct becomes smaller and smaller and begins 
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Figure 20. Coverage Diagram for a 9000 MHz airborne 
early warning radar transmitting within 
an elevated duct located bewteen 2,590 
feet and 3,220 feet above the surface. 

Note the extended coverage in the duct and 
the partial loss of coverage just above 
the duct. This coverage diagram like those 
in Figure 20 was produced by the Navy IREPS. 
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farther away from the aircraft until a point is reached 
where the hole is no longer apparent. 

Electronic countermeasures (ECM) , ESM, and DF systems 
performance, as well as air to air communications and air to 
ground communications are also affected by elevated ducts, 
with the effect being most likely at higher frequencies. It 
has also been noted that the propagation of EH signals from 
transmitters located well below the elevated duct is greater 
in the duct than expected. One explanation mentioned by 
Skolnik [Ref. 7: p. 453) postulates that the EM signals are 
being scattered into and out of the elevated duct by 
variations in the index of refraction profile due to the 
"wavy" top of the duct, causing it to be similar to a 
nonuniform waveguide. Figure 21 shows examples of the 
effect that the position of the aircraft (whether radar or 
jamming) relative to the elevated duct has on the trapping 
of EM waves and on the effect of the radar/radio hole on EM 
waves . 

For radar systems, the presence of surface or 
evaporation ducts can lead to an increased surface clutter 
return, both in magnitude and in coverage extent, masking 
targets over a much greater range. An example of this 
effect is shown in Figure 22. As described by Beach 
[Ref. 3: p. 85], if the clutter return is stronger than the 
target return at the same range, detection of the target may 
be seriously degraded or impossible. The strength of the 



57 



Radar or Jammer 




Figure 21a. Scenario of Jammer or Radar Source 
3elow Elevated Refractive Layer. 




Figure 21b. Scenario of Jammer or Radar Source 
Above Elevated Refractive Layer. 
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Target 




Figure 22. Air search radar geometry showing possible 
clutter return from rough surface at the 
same range as the air target for a surface 
duct . 
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clutter return is partially dependent on the strength of the 
duct. MTI radars can be seriously affected by the enhanced 
strength of the clutter return which may overcome the radar's 
ability to separate out a moving target from stationary 
clutter . 

The fading phenomena described earlier can also 
become more severe in the presence of ducting and can lead 
to the "skip" effect in the presence of a surface duct. This 
phenomenon is characterized by periodic spatial gaps in radar 
coverage or communications signals. These spatial intervals 
give the appearance of rings on the radar display. Figure 23 
shows the ray paths that characterize skip rings. Note that 
the appearance of skip rings is dependent on the surface 
reflectivity as well as the strength of the surface duct 
and the angle at which the EM ray is launched. Since water 
generally has a higher surface reflectivity than land this 
skip effect is more significant over the oceans of the world. 

Subrefractive conditions can cause the creation of 
blind spots within the subrefractive zone for EM emitting 
aircraft (airborne early-warning and EW aircraft) located 
above the zone, as shown in Figure 9. Aircraft traveling 
close to the earth's surface, below the top of the sub- 
refractive zone, can take advantage of this condition. In 
the case of airborne early-warning systems , penetrating 
aircraft will be able to move in much closer to the radar 
aircraft before being detected. As an example, Davidson 
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Figure 23. Skip Effects. EM rays from the transmitter 
will follow different paths depending on 
the angle they are transmitted at and the 
surface duct characteristics. 
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ERef. 5: pp. 3-10] discusses the effects noted by an E-2C 
squadron. One of the phenomena noted was that for aircraft 
flying at approximately 25,000 feet ray angles of as much as 
5 degrees below the horizontal were required to penetrate 
solar heated subrefractive conditions, seriously degrading 
the early warning ability of these systems. In a similar 
manner, air to ground communications can also be affected by 
subrefractive conditions . 

Optical wavelengths are affected by refraction in a 
manner similar to microwaves, but on a much reduced scale 
due to the relative ineffectiveness of water vapor gradients. 
Figure 8 shows that refraction does extend the optical 
horizon beyond the geometric horizon but not nearly to the 
extent that the radar/radio horizon is extended. Since 
there is usually little measurable difference in the 
refractive effects experienced by IR and visible wavelengths, 
which together form the optical region, they will generally 
be considered together in this thesis. 

Refraction in a standard atmosphere causes the horizon 
seen by the observer to appear higher and more distant than 
the geometric horizon. McCartney [Ref. 11: p. 100] describes 
a formula, given by Bowditch, for calculating the distance 
to the optical horizon in a standard atmosphere as 

Horizon distance = 3.839s km, (2.10) 
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where s is the height of the eye in meters. This relationship 
is shown graphically in Figure 24. Superrefractive and 
ducting conditions will extend the optical horizon even more, 
causing objects beyond and below the horizon to become visible. 
Extremely strong inversions near the surface of the earth will 
form a mirage, known as looming, leading to greatly extended 
optical horizons. This subject will be taken up in more 
detail in the section dealing with mirages. Figure 24 shows 
that optical waves traveling in the atmosphere are also 
subject to downward bending due to refraction, leading to 
range and angle errors in surveillance and tracking systems. 

The degree of bending is dependent not only on the refractive 
conditions prevailing but also on the zenith angle £ . As 
described by McCartney [Ref. 11: p. 112], when the zenith 
angle is 75 degrees or less the effect of refraction in a 
standard atmosphere is largely ignored. For larger values 
of the zenith angle the effects become apparent particularly 
if £ is greater:- than 90 degrees or if superrefraction or 
ducting conditions exist. When the beam path lies at such 
an angle, and is also tangent to the horizon, the refractive 
effect in bending and lengthening the path can be significant. 

The effects of refraction on millimeter wave 
propagation in both a standard and non-standard atmospheres 
is an area not clearly understood at this time. The 
millimeter wave region of the EM spectrum has only recently 
been actively considered for application to communications, 
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Zenith 




True Zenith Angle £ 
Apparent Zenith Angle 



S Height Of Eye 



H - Geometric Horizon 

H’ - Visible Horizon 

H" - Apparent or Optical Horizon 



Figure 24. Refraction Effects on Optical Transmission. 

The optical horizon is the visible horizon 
which appears to be displaced upward. Note 
also the angular and range error induced in 
optically tracking the aircraft by refraction. 
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surveillances and weapons system’s guidance. The majority of 
the information currently available is largely theoretical. 

In a standard atmosphere, there is some reduction of signal 
strength and, like microwave signals, there are errors 
induced in range and angle estimation of radar targets. 

These effects appear to be rather small in standard conditions. 
T. Patton et.al. [Ref. 12: p. 170] states that in non- 
standard atmospheric conditions (anomalous propagation) 
there is some evidence of trapping of millimeter waves by 
surface ducts, but there has been little investigation of the 
effects of elevated ducts. Also, the importance of a water 
vapor gradient in the atmosphere on millimeter waves has not 
been adequately addressed although it does appear to be more 
significant in MMW propagation than in microwave propagation 

and less so in optical propagation. D. Snider, in private 
3 

communications , states that due to the difference in the 
responses to humidity of optical and millimeter wave energy, 
if there exists, simultaneously, a temperature inversion, 
created by heat radiation, and a slight increase in humidity 
with altitude, optical rays would bend downward, due to the 
temperature inversion, while MMW rays would bend upward due 
to the humidity gradient. 



The discussion with D. Snider was in regards to the 
difference in the effects of atmospheric refraction on 
optical and MMW systems. 
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6 . Methods of Predicting Refractive Effects 

Calculating the refractivity of the atmosphere, 
detecting the presence of anomalous propagation and 
determining the effect of refractivity on EM wave propaga- 
tion is very complex and as yet inexact. Characterizing the 
variation of atmospheric refractivity with altitude is 
usually done by making measurements of temperature, pressure, 
and humidity as a function of altitude using a radiosonde 
balloon. Based on these measurements, and utilizing 
Equations 2.4 and 2.5, a determination of the refractivity 
and its gradient at various elevations can be made. This 
information is generally sufficient to determine the 
presence of anomalous propagation phenomena. 

The use of a radiosonde for measurement has a major 
drawback; it is generally a slow response instrument and 
somewhat inaccurate. Consequently the data collected by 
the radiosonde greatly limits the ability to accurately 
calculate refractive index gradients. Ongoing efforts are 
underway to develop and deploy smaller and more accurate 
means of measuring refractivity, directly or indirectly. 
Kaelin [Ref. 13: p. 124] in his thesis states that examples 
of new systems are the mini-refractionsonde and airborne 
dropsonde being developed by the Naval Air Development 
Center, Warmister, Pennsylvania. A direct and accurate, but 
not commonly used method of measuring atmospheric 
refractivity, is the airborne microwave refractometer . 
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A thorough and accurate table of climatological 
records that describe local geographic conditions will also 
provide a valuable means of at least probabilistically cal- 
culating the expected refractive profiles for a region. 
Forecasts of anticipated weather conditions can also be used 
to develop a general idea of expected refractive conditions. 
These two methods, one for prediction and the other for 
forecasting of refractive conditions will both be discussed 
further in Section E of the thesis. 

To determine the effect of the existing refractive 
conditions, the frequency of the EM wave of interest must be 
known. The Navy’s IREPS will take the condition of the 
atmosphere and the wavelength of interest and, using a 
desk-top computer, graphically display the refractive index 
calculated and the effect this condition has on various 
radar and radio systems working in the VHF, UHF, and 
microwave region of the EM spectrum. Likewise, the Navy is 
developing a computer based model to assess atmospheric 
effects (including refractivity ) on electro-optical system’s 
performance. This system, described by F. Snyder [Ref. 14], 
is called PREOS (Prediction of Performance Range for 
Electro-Optical Systems). The Navy eventually intends to 
unite IREPS and PREOS into a single prediction system. 

Figure 25 provides an example of a PREOS output and Figures 
19 and 20 are examples of IREPS outputs. 
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Berore methods can be used to compensate for, correct, 
or take advantage of refractive effects on EM propagation, 
two things are necessary. First, accurate prediction and 
forecasting are needed. Some prediction means have been 
described here and in a later section prediction methods and 
forecasting methods will be described in detail. Secondly, 
an in-depth understanding of the effects of refraction on 
specific wavelengths and even specific systems that utilize 
propagating EM energy is needed. In the second half of the 
thesis several systems will be used to provide examples of 
refractive effects* wavelength dependence. 

D. MIRAGE 

One form of refraction of visible light that occurs in 
the lower atmosphere, with which almost everyone is 
familiar, is the mirage. Before discussing the effects of 
the several types of mirages on EM energy transmitted 
through the lower atmosphere, a general description of the 
conditions leading to mirage formation will be addressed. 

The majority of the information used in this discussion 
comes from an article by A. B. Fraser and W. H. Mach in a 
1980 issue of Scientific American [Ref. 15: pp. 29-37]. 

1 . General Description of the Mirage Phenomena 

As with other forms of refraction, the atmosphere 
creates mirages by the bending of EM energy. As previously 
discussed, the index of refraction in the visible spectrum 
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is largely dependent on temperature, as seen in Equation 2.5. 
The greater the temperature gradient, the stronger the 
gradient of refraction and the greater the amount of bending 
of light. The portion of the atmosphere responsible for 
forming mirages is located in the first several meters above 
the surface. The exact form of the mirage is highly variable 
and is dependent on the solar absorption at the surface, the 
temperature gradient, the position of the object relative to 
the temperature gradient, and the position of the observer. 

Light traveling through the shallow lower layer of 
the atmosphere will, due to the existing temperature 
gradient, travel in a parabola. The curvature of the ray can 
be upward or downward, depending on whether there is a 
decrease or increase of temperature as altitude increases. 

A general way to describe this phenomenon is to draw two 
light rays from an object, one moving directly to the 
observer's eye and the other following a parabolic path, 
as shown in the two-image inferior mirage in Figure 26. 

The first ray presents the light directly from the object. 

The second ray path is bent upward (due to a decrease in 
temperature at higher altitude) and thus appears to be light 
from an object which is located vertically below the actual 
ob j ect . 

The simplest types of mirage to distinguish between 
are superior and inferior. In a superior mirage the 
temperature increases with height, as in superrefraction or 
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Figure 26. 
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General Mirage Effect. Light ray A from 
object P is above the subrefractive condition 
and travels in a straight line to the observer 
C. Light ray B from object P travels through 
the subrefractive condition and thus creates 
a virtual image P* which is seen by the 
observer C. 
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ducting, and a horizontal surface will appear to be concave 
upward, presenting an impression to the observer, of being 
inside a shallow bowl. Figure 27A is a ray diagram of a 
superior mirage, showing that the image of the observed 
object is displaced upward. According to McCartney [Ref. 11 
p. 100], although most mirages only involve relatively short 
distances (1-5 km), a superior mirage created by a very 
strong inversion can make objects located far beyond and 
below the optical horizon visible in a phenomenon known as 
looming. 

An inferior mirage is created by a gradient in which 
temperature decreases with height (subrefraction). As 
described by Fraser and Mach [Ref. 15: p. 35], a horizontal 
surface altered by an inferior mirage will appear to be 
convex upward. The impression is also of a bowl, but this 
time the bowl is inverted and the observer is on top of it. 
There is thus an optical horizon created beyond which the 
surface being observed disappears. This phenomenon is 
referred to as sinking. The decrease of temperature, with 
altitude, which leads to this profile is common over 
enclosed bodies of water in the early morning and in desert 
regions in the late afternoon. Note that, with an inferior 
mirage, light rays from the bottom of the object may not 
reach the observer if this portion of the object is beyond 
the optical horizon, as shown in Figure 27b. 
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Figure 27b. 













Superior Mirage. The virtual image p 1 ' is 
seen by the observer A as being above the 
actual object P. 




This time the virtual image P* is seen by 
the observer A as being below the actual 
object P. 
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If the temperature gradient is constant with height, 
no distortion of the image will occur in a mirage, since 
displacement of all portions of image, upward or downward, 
will be the same. However, such a situation is not the 
general case, since some magnification or reduction of the 
image will usually occur. When an image is reduced in 
height, the condition is referred to as stooping. This 
occurs if the temperature gradient decreases as the tem- 
perature increases with altitude. Such a condition often 
occurs over lakes or inlets on a sunny afternoon, where 
warm air from the land is carried over the cooler water. A 
superior mirage is thus formed with the bottom of the image 
displaced upward more than the top of the image. The result 
to the observer is a squashed image of the object, and this 
image will appear to fluctuate in vertical dimension, as the 
day progresses. 

A magnification of the image, or towering occurs, 
according to Fraser and Mach [Ref. 15: p. 35], in an inferior 
mirage when the image is displaced downward, but the bottom 
of the image is moved downward more than its top is . This 
happens when the temperature gradient increases as the 
temperature decreases with altitude and occurs over enclosed 
bodies of water in the early morning or over sun-heated 
ground later in the day. 
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The familiar mirage of the desert is the two-image 
inferior mirage. The "water" seen by the observer is 
actually an inverted image of the sky seen below the horizon. 
The same atmospheric phenomena that create the inferior 
mirage and towering create the two-image inferior mirage , 
i.e., the temperature and temperature gradient both decrease 
with height. The difference is that the temperature profile 
must be stronger, thus resulting in the ray being bent more 
significantly. Since the object rises so far above the 
surface that it is near the top of the temperature gradient, 
rays will travel both directly from the object and in a 
parabolic manner, similar to the simplified diagram in 
Figure 26. Because the ray of light passing through the 
region of strong temperature gradient is so strongly bent, 
it will not join the eye with the bottom of an object but 
will join the eye with the top of an object to provide a 
second, but inverted, image. As Fraser and Mach [Ref. 15: 
p. 35] state: 

"The image is inverted because as the observer lifts his 
gaze slightly he is looking through a region of the 
atmosphere that has a weaker temperature gradient , so 
that the ray is less strongly curved. It will therefore 
join the eye to a point lower on the object rather than 
higher, as would usually be expected." 

An object moving away from the observer will first appear as 

a single image, then as a double image, and then vanish 

from the bottom up. The presence of a double image will 

make the ground between the observer and the object appear 
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to be a water surface. The distance to the point where the 
"water" begins is the location of the optical horizon. As 
the temperature cools, the optical horizon is extended and 
the " water” recedes and, finally, vanishes. 

According to Fraser and Mach [Ref. 15: p. 36], the 
most common form of the three-image mirage is caused by 
temperature increasing with height, for a distance, and then 
decreasing. This condition will occur over enclosed bodies 
of water on warm afternoons, and results in the surface of 
the water appearing as a "large, flattened letter S" - thus, 
the surface appears to fold over. If the inflection point 
is gentle, the surface will appear to rise up rather than 
fold over. An object at the distance where this occurs will 
be magnified, creating an example of the "Fata Morgana". 

Fraser and Mach [Ref. 15: p. 36] feel that the 
occurrence of gravity waves in the atmosphere, along with the 
conditions leading to a three-image mirage, can result in 
the creation of the classic "Fata Morgana", which will 
sometimes appear as a hugh wall rising up above the surface 
or as a mountain range in the distance. This combination of 
stratified atmosphere and gravity waves causes a periodic 
oscillation of the refractive index. Energy imparted from 
surface winds and a countering force imparted by the earth's 
gravity are responsible for the waves. The oscillating 
bands of varying refractive indices will lead to the 
appearance of bright pilasters and dark windows in the 
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wall created by the "Fata Morgana", and are actually the 
image of the surface below the mirage. 

Although there are other less common forms of mirage , 
those described above show some of the great variety of 
types and the atmospheric conditions leading to their 
creation. Such conditions can potentially disrupt the 
operation of optical systems. 

2 . Effects of Mirages on Optical Devices 

As previously stated, there are a variety of forms 
of mirage, all dependent on the conditions of the atmosphere 
near the surface of the earth. The resulting refractive 
conditions are similar to those of anomalous propagation 
but are on a much smaller scale, though the line separating 
anomalous propagation from mirage is often unclear. The 
basic difference between these two forms of refraction are 
the local nature of the atmospheric conditions leading to 
mirage formation and the larger scale atmospheric conditions 
leading to anomalous propagation. 

The effects of mirages are most apparent on optical 
systems. Generally, there is little difference in the 
effects of mirages on IR and visible wavelengths, even in 
the presence of strong refractive conditions. The effects 
on optical tracking systems can be quite significant, 
resulting in extension (superior mirage) or shortening 
(inferior mirage) of the optical horizon and distortion of 
the image (towering and sinking). Mirages such as the three 
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image mirage and the "Fata Morgana" can appear very solid 
and real, and, therefore, be quite misleading. Just as 
light from the object to the observer is strongly bent, so, 
too, will be any visible energy directed from the observer 
to a target in the mirage. 

The effects of mirages can be significant for optical 
trackers, visible surveillance, and laser systems operating 
at low elevation angles within the mirage. 

In a similar manner, optical systems using the IR 
spectrum will also be affected by refraction in mirage 
conditions. The PREOS model [Ref. 14] incorporates the 
ability to predict the amount of image distortion due to 
mirage effects. 

The effect of mirages on radio wave and microwave 
systems will depend not only on the general type of mirage 
(inferior or superior), but also on the strength of the 
refractivity gradient and the mirage thickness. The effects 
will be the same as those predicted for surface based 
anomalous propagation, whether superrefractive , ducting, or 
subrefractive . The effect on microwave systems, such as 
radar, will only occur at very low elevation angles and at 
frequencies high enough to be above f , (Equation 2.9). 

Prediction of mirage occurrence is difficult, since 
it is dependent on terrain, large scale atmospheric, and 
local atmospheric conditions. Mirages are generally 
associated with partially or fully enclosed large bodies 
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of water, such as lakes, sounds, or bays, and with hot desert 
regions. The probability of occurrence of mirage conditions 
and examples of the effects of mirages on specific systems 
will be addressed later. 

E. TURBULENCE 

Anomalous propagation and mirage conditions are created 
by relatively large scale, steady state variations in the 
refractive index of the atmosphere. Such refractive 
conditions affect EM signals at both the long and short end 
of the spectrum, from radio frequencies to visible light. 

The air is, however, never entirely motionless. Thus, much 
smaller scale variations in the index of refraction also 
exist. These random fluctuations of the index of refraction 
of the air are created by turbulent motion in the atmosphere 
coupled with temperature and humidity gradients. This 
turbulent motion includes inhomogeneities in the atmosphere 
varying in size from hundreds of meters in extent of those 
on the order of 1 mm in length. The effects of turbulence 
are largely confined to the shorter wavelengths (MMW, IR, 
and visible), due to the size of the refractive 
inhomogeneities relative to the size of the wavelengths. 

1 . General Description of Turbulence 

The inhomogeneities in the atmosphere consist of 
varying sized regions (discrete cells) which fluctuate in 
temperature and water vapor content in time and in space. 
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According to Zuev [Ref. 16: p. 348], variation of 1 degree 

centigrade will result in a change in the refractive index 

of approximately 1 x 10 . Fluctuations of the temperature 

of the air at a given location may be as high as tenths of 

a degree and can vary over periods of milliseconds to seconds. 

Zuev further states that the variation of temperature over a 

horizontal path can be several degrees for points on the 
2 3 

order of 10 to 10 m apart. 

These inhomogeneities are themselves constantly 
changing in size and location. They are created and modified 
by wind shear and convection arising from vertical temperature 
gradients and from vertical water vapor gradients over lakes 
and oceans. As stated by F. Lutomirski et . al. [Ref. 17: 

D. 1 ] 

"The most important factors affecting daily variations (in 
the index of refraction) are short wave solar radiation and 
long-wave terrestial and atmospheric radiation, the mean 
wind speed, type of soil and ground cover, and ground 
surface roughness." 

The overall effect of the inhomogeneities present in 
a turbulent atmosphere is to modify the characteristics of 
optical propagation through the atmosphere. The nature of 
the refractive effect is dependent on the operating wave- 
length, pathlength, and nature of the turbulence. Specifi- 
cally, the size of the inhomogeneities is important to the 
overall effect of turbulence on optical and MMW transmissions. 
The size of the inhomogeneities created by wind shear and 
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convection will vary greatly from those hundreds of meters 
in size (outer turbulent scale) L , to those, according 
to Zeuv [Ref. 16: p. 349], on the order of 

"V 3 

1 = / v /e , (2.11) 

o 

where 1 is the inner turbulent scale, v is the kinematic 
o 

viscosity, and e is the rate of dissipation of turbulent 
kinetic energy. The inner scale 1 , is the smallest size 
inhomogeneity that can exist due to viscosity while L q is 
the largest such region that can exist in which temperature 
can be expected to remain nearly constant. According to A. 
Cooper [Ref. 18: p. A-40], between these two scales there 
often exists a region of well established turbulence which 
can be described by the Kolomogorov spectrum as 

8 (K)~ K -11/3 , (2.12) 

where K = 2ir/l is the "turbulence wave number" representing 
the size of the eddy or inhomogeneity and 9 (K) is the 
"intertial subrange" of well established turbulence. 

The structure function of the refractive index between 
two points (r and r^) a me< 3ium wiTh random inhomogeneities 
is 

D (r) = [n(r, + r) - n(r-,)] 3 (2.13) 

n 1 1 
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where r is the distance between the two points. In a 
similar manner Zuev [Ref. 16: p. 3493 defines the temperature 



structure 


function 


D T (r) 


= [T(r 1 + r) - T(r x )] 2 (2.14) 



Atmospheric turbulence is generally described by the 
Kolomogorov model for fluid velocity turbulence. For the 
Kolomogorov turbulence model the two Equations 2.13 and 2.14 
can be written equivalently for small scale turbulence, as 



D T (r) 


= C T 2 r 2/ 3 , (2.15) 


D (r) 
n 


= C 2 r 2/3 , (2.16) 

n 5 


2 

where C 

n 


2 . 

is the refractive index structure value and C T is 



the temperature structure value. These relationships in 
Equations 2.15 and 2.16 only hold for limited separation 
values , when 

1 << r << L (2.17) 

o o 

The two structure constants are related by 



c 2 = 

n 


(79P/T 2 x 10 -6 ) 2 C t 2 , (2.18) 
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and are the usual means of describing the turbulence level. 

According to F. Hall [Ref. 19: p. TuCl-1] C 2 can also be 

c n 

described as the difference in refractive index between two 
points, squared, averaged and divided by the 2/3 power of 
the distance between them, or 

C 2 = (n - n 9 ) 2 /r 2 ^ 3 , (2.19) 

n 1 2 

where r is the distance between the two points. Note that 
2 

has the dimensions of length to the -2/3. Generally 

2 

is the only parameter needed to describe small scale 

turbulence and, for this reason, is the parameter most often 

used in studying optical turbulence effects. Actual 

2 

measurements of can be made by determining the difference 

in temperature between two sensitive, fast response thermo- 
meters located some distance apart or from the temperature 
power spectrum from one thermometer. 

Values of C in the lower atmosphere generally lie 
—16 —12 

between 10“ cm and 10 cm and vary widely with time and 

2 • • 
location. C values, and therefore turbulence conditions, 
n 

are considered to be light around 10~ 15 and severe around 10 1 

2 

Variation in the value of C n occurs not only along a hori- 
zontal path, but also along a vertical path from the earth’s 
surface to the top of an inversion. Turbulence will generally 

be strongest at the surface and at the top of an inversion. 

2 

For such a vertical path normally decreases in value as 
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altitude increases until the inversion region is reached 
2 

where values increase dramatically. Above the inversion 

there can also be significant turbulence depending on the 
wind and large scale weather conditions. 

The usual result of the many factors affecting 
turbulence formation is the turbulence maximum in the early 
or midafternoon when heating of the surface is greatest and 
winds are light, and minima at sunrise and sunset when little 
heating occurs. Turbulence will occur at night, particularly 
if an inversion develops, but only if a light wind is present 
and will generally be weaker than during the day. Thus, 
turbulence shows a characteristic diurnal variation as shown 
in Figure 28, although the diurnal variations will be weaker 
when the sky is overcast. Table 3 provides a rule-of-thumb 
for predicting the occurrence of turbulence given existing 

meteorological and surface conditions. Figure 29 through 31 

2 ... 

are actual plots of values for three different climatic 

conditions. Figure 29 is a desert region in summer. Figure 
30 is a treeless area over snow, and Figure 31 is a wet 
grassland typical of many regions in Central Europe. Note 
the distinct diurnal cycle in each plot and that the 
strongest turbulence occurs over the desert region. 
Significant turbulence is recorded in the other two regions 
as well, particularly over the wet grassland. 
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Figure 28. Diurnal Variation of Turbulence. The graph 
shows the measured values (solid lines) of 
optical turbulence at 1.5 meters above a 
grassy surface in a semi-arid region and 
those from bulk meteorological measurements 
(dashed lines). 
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TABLE 3 



INTERIM RULES-OF-THUMB 

Minimum thermal turbulence in "the boundary layer is 
observed with the following atmospheric and earth-surface 
conditions : 



o Under a heavy overcast, day or night. 

o Over water or wet ground, day or night. 

o With winds exceeding 2 8 m sec - ''" ( s 15 kn), nearly 

regardless of surface heating or cooling, day or night. 

o During clear weather, for periods of 10 to 30 minutes 
shortly after sunrise and before sunset. 

o During generally unstable days, for brief periods in 
the shadows of transient clouds (the same effect maybe 
obtained. during nights that are generally stable). 

o In descending volumes of air, specifically between 

convective plumes in daytime, f ree-convection conditions. 

o On clear nights, with absolutely still air (any air 
motion can drastically change this situation) . 

Maximum thermal turbulence is generally found with the 
conditions : 



o Over relatively dry ground, with strong solar heating 
and light winds . 

o On clear nights with light but nonzero winds (the most 
unpredictable, sporadic turbulence is found under these 
conditions, 

o In the upper atmosphere, in stable strata associated 
with vertical wind shears. 
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Hour 

Figure 29. C vs mean solar time at .5, 2, 8, and 16 meters (top to 
bBttom) for a dry desert in summer. 
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Hour 

Figure 31. C vs mean solar time at .5, 2, 8, and 16 meters (top to 
bottom) for a wet green short-grass grassland. 



2 . Principal Effects of Turbulence on Optical Propagation 

Turbulence in the atmosphere, along the optical path, 

affects the optical properties of focused beams and imaging 

2 

systems. Point measurements of along an optical path 

often show wide variations in values due to temperature and 

water vapor variations. Due to this spatial variation in 
2 

, a path integrated value is usually determined for the 

overall path. Often in calculations of the effects of 

2 

turbulence on optical systems the value of C n used will be 

weighted for the calculation. These weighting functions will 

emphasize the effects of certain parts of the path (the 

center or ends) and may totally ignore other regions. The 

value of weighting will vary with the type of turbulence 

effect being considered, thus, care must be exercised in 

2 

using "weighted" values of in calculations. 

Instantaneous deviations in n due to turbulence can 
lead to a number of effects on optical systems to include: 
scintillation, wave front tilt or fluctuations in angle of 
arrival, beam spread, beam steering or beam wander, and 
degradation of the Modulation Transfer Function (MTF) or 
loss of resolution. Each of these effects will be considered 
separately, although it will be seen that many are closely 
related. E. Crittenden [Ref. 20: p. 2] states that H.T. Yura 
in his analysis of turbulence effects shows that the effects 
on an active focused beam system will be applicable to those 
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imaging systems as well, although the manifestation of 
the effect may be somewhat different, 
a. Scintillation 

Scintillation or beam breakup is due to the 
fluctuation of the intensity of a focused beam or image 
which has propagated through a turbulent medium. This 
characteristic sometimes referred to as "image boil" is 
caused by random bending of optical rays resulting in 
destructive and constructive interference patterns. These 
patterns vary with position and time. Figure 32 shows an 
example of scintillation of a laser beam. The twinkling of 
stars viewed from the earth is another example of 
scintillation . 

The variance of the log amplitude intensity of 

a laser beam, as seen through a pinhole at the receiver, can 

2 

be related to , the propagating wavelength, and the 
pathlength, according to E.C. Crittenden et . al . [Ref. 20: 
p. 16] by 



2 

a LnA 



0.124 



. 7/6 
k 




11/6 



o 



( 2 . 20 ) 



where k = 2ir/X and Z q is the propagation pathlength. As 
expected, the greater the level of turbulence or the longer 
the pathlength, the greater the scintillation observed. It 
should also be noted that the shorter the optical wavelength 
the greater the scintillation effect. Equation 2.20 is only 
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Figure 32. An ’’open camera" photograph of a 
scintillating beam profile. 
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o 

accurate for small amplitude fluctuations, for 2a T A <.8. 

For higher turbulent levels. Equation 2.20 is used only to 
infer the approximate relationships between C n , pathlength, 
wavelength, and log amplitude variance. 

Turbulence levels generally vary along the 
propagation path, rather than being homogeneous as implied 
by Equation 2.20. This fact must be considered in calculating 
the actual log amplitude variance, and Equation 2.20 can be 
rewritten (for spherical waves) as 

°LnA 2 = * 56 k>7/6 f ° C n 2 w dz » (2.21) 

o 

2 

where w is a weighing function that emphasizes the 

values in the middle of the optical path and ignores those 

values at the path ends . 

2 

C values can also be calculated from measure- 
n 

ments of scintillation. According to Crittenden [Ref. 20: 
p. 16], however, in cases of nonuniform turbulence levels 
along the path of propagation such calculations may not be 
of any use in calculating other turbulent effects because 
of the different weighting functions involved, 
b. Wave Front Tilt and Beam Spread 

Wave front tilt or fluctuations in the angle of 
arrival of a beam or image are caused by bending of the 
entire wave front of an image or beam as it passes from one 
inhomogeneity to another. As would be expected, the inner 
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scale turbulent inhomogeneities have little effect on 
creating wave front tilt, but rather, it is the larger, 
outer scale, turbulence regions that lead to this phenomenon 
The result of this effect is to cause "image dance", a move- 
ment of the apparent position of an object's location in an 
imaging system. This is the phenomenon that leads to the 
false apparent location of stars viewed from the earth's 
surface . 

One characteristic of beam spread is the 
increasing of the width of an initially narrow beam as it 
propagates through the atmosphere. This expansion of the 
beam width and resultant decrease in intensity is due to the 
random refraction of the beam as it passes through a turbu- 
lent atmosphere. In essence, relatively small scale wave 
front tilt is occurring at many different locations across 
the beam profile. Long term beam spreading can be separated 
into two components: beam expansion as described above and 

beam wander. Beam wander is the movement of a beam from a 
straight path due to turbulent inhomogeneities that are much 
greater in length than the beam width, thus, refracting the 
entire beam. As expected, beam wander is a function of the 
beam diameter and the turbulence scale. Lutomirski [Ref. 17 
p. 115] states that beam wander is closely related to wave 
front tilt and in cases where the only significant effect on 
a beam propagating through the atmosphere is wave front tilt 
the instantaneous beam wander angle and wave front tilting 
angle will be equal. 
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c. Resolution 



A cumulative effect of the various turbulence 
effects is an overall degradation of the image resolution of 
an object seen through the atmosphere. Characteristics of 
an object are said to be resolved if their images are 
discernible as being separate. This loss of resolution is 
characterized by a degradation of the Modulation Transfer 
Function (MTF) of the atmosphere by turbulence. The MTF 
describes the detail of an image, transmitted from one point 
to another, in terms of the spatial frequency. The MTF for 
a total optical system can be calculated by determining the 
MTF for the components (lens, mirrors, and atmosphere) and 
forming the product of these component values. Knowledge 
of the MTF of the atmosphere is important in predicting the 
effects of turbulence on low, medium and high power laser 
beams leading to beam wander and modification of the beam 
profile. However, in high power lasers the heating of the 
air must also be considered. To be entirely accurate, when 
discussing the MTF as it relates to optical systems, it is 
the Optical Transfer Function (OTF) that is discussed. The 
OTF describes both the phase and amplitude changes in the 
atmosphere. When only the amplitude changes are considered, 
the OTF is called the MTF. 

Knowledge of the MTF can also aid in predicting 
the turbulence contributions to resolving ability, image 
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broadening, and image wander of imaging systems. Using 

weighting functions that emphasize the ends of the path, MTF 

2 

curves can be developed from measured values of C and X. 

n 

For imaging systems , the Mutual Coherence Function 
(MCF) is the equivalent of the OTF and, in fact, the atmos- 
pheric MCF is the same function as the OTF but expressed in 
terms of different variables. The MCF is of considerable 
concern in coherent detection systems. It reflects the loss 
of transverse coherence of an initially coherent wave as it 
passes through the turbulent atmosphere. Knowledge of the 
MCF can also be used to provide information on atmospheric 
effects on laser propagation. 

The OTF for a propagation system, along a path 
where the turbulence varies is described by Crittenden 
[Ref. 10: p. 14] by 

M(f) = exp{ -57.64 x“ 1/3 f 5/3 / ° C r 2 (z/z q ) 5/3 dz} , 

o 

( 2 . 22 ) 

where z q is range, f is the angular spatial frequency, 
defined as f = 1/L, where L is the spatial period of the 
distribution of the modulation of the amplitude (intensity) 
of the optical beam, andxis the wavelength. The same 
function M(f) is used to calculate the MCF of an imaging 
system but using different variables. 
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The most common parameter used to describe FLIR 
propagation through the atmosphere is the Minimum Resolvable 
Temperature Difference (MRTD) . The MRTD implies the use of a 
standard four— bar pattern with an individual bar length to 
bar width ratio of 1:7. The MRTD is the minimum temperature 
difference between the bars and the background for which the 
bar pattern can be resolved at the 50% probability of detec- 
tion level for a given FLIR. Figure 33 shows a four-bar MRTD 
test pattern such as that discussed above. As with other 
means of measuring resolution, the MRTD for a FLIR system 
can be greatly degraded by turbulence. 

Figure 34 provides a comparison of line spread 
functions measuring the atmospheric OTF for two different 
levels of turbulence. Note that several of the effects of 
turbulence are occuring simultaneously. The cumulative 
effect of scintillation, image broadening, and image wander 
(wave front tilt) is to reduce the resolution of an image or 

degrade the shape and move the end point location of a beam. 

2 

As mentioned earlier, C can be calculated from 

’ n 

measurements of scintillation. Likewise, Crittenden et . al. 

[Ref. 21: pp . 13-1 to 13-2] states that measurements of the 

2 

OTF and beam wander can also be used to calculate C and 



with more ease than using scintillation. 



n 
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Temperature 



Figure 33. 




.Temperature 

/ 



T + AT 



Standard Four-bar MRTD test pattern where 
the bars represent a blackbody source of 
temperature T+AT and the spaces represent 
the background temperature T, and AT is the 
difference between the background temperature 
and the targets temperature. 
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I 




Figure 34, 



Line spread spectrum for two levels of 
turbulence. These photos show that associated 
with an increase of C n 2 is a broadening of the 
image related to the OTF, increased detailed 
structure (scintillation) and an increase in 
the wander of the center (image wander). 
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3 . Manifestation of Turbulence in Optical Systems 
In beam forming optical systems, such as laser 
systems, scintillation may result in a few dB to upwards of 
40 dB fluctuation above and below the average irradiance 
value. For a laser communications system this effect can 
be significant in degrading the bit error rate. Scintillation 
is less important at IR and longer wavelengths than at 
visible wavelengths. According to Lutomirski [Ref. 17: p. 230], 
recent studies have shown that at the longer wavelengths (IR) 
a coherent detection laser communications system is preferable, 
while at visible wavelengths, an incoherent detection system 
is best when turbulent effects are concerned. Beam wander can 
cause significant enough movement of a laser beam to degrade 
detection in an optical communications system or to sweep the 
beam off the target in a laser designation or laser weapon 
system. Wave front tilt will also cause a degradation of 
detection in an optical receiver, although this may be the 
simplest turbulence effect to correct. Beam spread can 
result in a reduction of the irradiance of a laser weapon and 
a sensitivity reduction in laser communications systems. The 
MCF can be used to determine the limiting mean irradiance 
from an initially coherent beam and the turbulence limited 
S/N ratio of an optical heterodyne (coherent) receiver. 

Table 4 consists of two tables listing several types of army 
laser systems and the types of turbulence induced degradation 
each system can suffer from. 
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POSSIBLE TURBULENCE-INDUCED DEGRADATIONS FOR VARIOUS SYSTEMS 
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TYPICAL TACTICAL LASER APPLICATIONS 
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Turbulence will create the same effects on imaging 
systems but will often be manifested in a different manner 
than the effects on beam producing optical systems . 
Scintillation will create variations in the intensity of the 
image of the object being viewed and will thereby degrade 
detection and identification. Image wander can cause a 
degradation of detection ability, dependent on the type of 
imaging system used, specifically the nature of the detection 
method used by the imaging system. Image broadening may 
degrade both detection and identification. Wave front tilt 
can present a false apparent location of the object. Tur- 
bulence induced degradation of the resolution of an image 
can be quite severe, affecting detection, discrimination 
and identification. The limiting resolution due to 
turbulence can be determined from the MCF. 

It is apparent from this discussion that turbulence 
effects on optical systems are varied, but generally result 
in a degradation of the system’s operation. Some means have 
been developed to compensate for turbulence effects, a few 
of which will be mentioned. 

Lutomirski [Ref, 17: p. 65] states that perhaps the 
easiest form of turbulence effects to correct is wave front 
tilt to the entire beam or image, if it occurs largely by 
itself. Mechanical or electronic tracking of the image can 
be used in a system to correct the displacement in the image 
plane, producing a nondegraded image. 
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Other forms of turbulence effects are not so easy to 
correct, particularly when several occur at once. Because 
of the small scale of turbulence effects, turbulence 
correction requires a resolution element related to the size 
of the turbulence scale. These elements must be small and 
response times very short, due to the rapidly changing 
nature of turbulence effects. This is the basis of adaptive 
optics, the purpose of which is to minimize turbulent 
effects on optical systems, both imaging and beam forming. 

The adaptive optics system senses phase shifts induced in 
the optical beam by turbulence and introduces a compensating 
phase shift to correct for this. In particular, the Coherent 
Optics Adaptive Technique (COAT) is basically a self-phasing 
phased array system. Errors in the phase of the incoming 
signal are sensed across the wavefront and corrected by the 
COAT system. The COAT output aperature is divided into many 
small segments whose phase can be controlled independently 
using mirrors or splitters, providing the means to correct 
phase errors in the incoming beam. This technique has been 
shown to improve to some degree the degrading effects of 
scintillation, beam spread and beam or image wander but can 
not as yet fully improve the system’s performance to that at 
low turbulence conditions. 
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^ . Turbulence and MMW Propagation in the Atmosphere 

Although turbulence has been largely studied because 
of its performance degradation of optical systems , it has 
been theorized for some time that it may also adversely 
affect millimeter waves. Recent tests have confirmed this 
fact and have also shown agreement with theory in that there 
is a stronger dependence on the atmosphere’s humidity at MMW 
than at optical wavelengths. Thus, in calculations involving 
the effects of turbulence on MMW propagation, both C^, , the 
temperature structure value and Cq, the humidity structure 
value, must be considered. Additionally, the cross corre- 
lation of these two parameters, C T q, must be included in 
determining turbulence effects. Scintillation and 
fluctuations in the angle of arrival of the beam in radar 
applications have been shown to be the primary turbulence 
effects at millimeter wavelengths. McMillan and Bohlander 
et. al. [Ref. 22: pp. 32-39], in their 1982 article on MMW 
propagation, provide a detailed discussion of the formulas 
used in calculating intensity fluctuations (scintillation) 
and angle of arrival fluctuations. 

As described by McMillan and Bohlander [Ref. 22: 
pp . 33-36], results of tests at White Sands Missile Range 
in New Mexico show that scintillation effects, though present, 
are generally not expected to be significant on MMW radar 
systems. However, fluctuations of the angle of arrival of 
MMW beams are on the same order of expected angular accuracy 
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and, thus, can be significant in degrading angle tracking 
ability as shown in Figure 35. The results of the testing 
seem to indicate that turbulence effects on angle tracking 
may be smaller than those due to overall anomalous 
propagation effects. Although more data on this subject is 
necessary before the full nature of turbulence effects on 
MMW is known, the results to date appear to indicate that 
turbulence should be considered when discussing MMW 
applications . 

F. PREDICTION AND PROBABILITY OF OCCURRENCE OF ANOMALOUS 

PROPAGATION, MIRAGE FORMATIONS, AND TURBULENCE 

The effects of atmospheric refraction on EM wave 
propagation, regardless of severity, are of little real 
consequence unless these conditions occur with some 
frequency. Furthermore, to respond to these effects, it 
is necessary to be able to forecast their occurrence and/or 
determine their existence. Some discussion on meteorological 
conditions and characteristics of terrain that lead to the 
three major forms of atmospheric refraction has already been 
dealt with and will not be reiterated here, except, as 
needed, for clarity. Rather, this section will attempt to 
look at the larger scale causes and characteristics of 
atmospheric refraction and focus on the probability of 
occurrence of anomalous propagation, mirage conditions, and 
turbulence in various regions of the world. Finally, two 
specific regions. Central Europe and the Mid-East will be 
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Peak to Peak angle of 




Figure 35. Calculated peak-to-peak angle-of-arrival 
fluctuations for millimeter wave signal. 
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looked at in some detail, not only because of the obvious 
military interest, but also to compare these two distinctly 
different areas. 

Anomalous propagation will be examined first, for two 
reasons: first, because a larger body of information exists 

on this subject; second, because much of what is discussed 
concerning anomalous propagation can be directly related to 
turbulence and mirage forecasting and prediction. 

1 . Anomalous Propagation 

A great deal of study has been conducted towards 
understanding the causes of anomalous propagation in its 
many forms, and in determining analytical methods of 
predicting its existence. Most of these have already been 
discussed. Forecasting anomalous propagation conditions can 
also be performed, given an understanding of the location 
and terrain features of the region of interest, the time of 
year and time of day, and existing large-scale weather 
conditions (air masses, high and low pressure regions, etc.), 
many of which interact closely. 

The general nature of the earth's surface, be it a 
large body of water or a large land mass, contributes to the 
expected refractive characteristics. This is due to the 
physical nature of the surface (mountains, desert, shallow 
tropical sea, etc.) and the effect of the surface on the air 
masses above. The persistence of superrefractive and ducting 
conditions over the oceans of the world is due largely to the 
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constant presence of moisture at low altitudes and the 
relatively slow changes of large air masses over these 
bodies of water. This persistence over water tends to 
decrease with distance (north or south) from the equator. 

Figure 36 shows this fact graphically. 

Over large land masses there generally is no such 
consistency in atmospheric conditions as exists over oceans. 

The humidity content and temperature of the air, as well as 
the large scale weather conditions (fronts and pressure 
regions) change much more rapidly over land than over sea. 

For this reason, persistent ducting is much rarer over land, 
and forecasting of anomalous propagation is generally more 
diff icult . 

The presence of certain terrain features also tends 
to have an effect on the occurrence and characteristics of 
refraction. High mountains will significantly affect the 
movement of air masses across the surface of the earth and 
impact on the creation of trade wind inversions and land-sea 
breeze recirculation. C. Samson [Ref. 9: p. 7] states that 
large river valleys may show patterns of refractivity 
differing markedly from surrounding regions due to 
stratification of the atmosphere by air drainage, above the 
valleys. Heating of large arid surfaces will often result 
in subrefractive conditions. 

Large scale weather conditions are the most significant 
factor in forecasting overland anomalous propagation. Air 
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(a) February 




(b) August 



Figure 36. Regions of persistent oceanic atmospheric 
duct layers. 
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mass characteristics, the presence of high pressure regions, 
and the movement of fronts and cyclones are all important 
in the formation of anomalous propagation conditions . 

Air masses are regions of air in which the properties 
of the atmosphere are relatively uniform. These air masses 
are formed in source regions which include large desert areas, 
tropical seas, and polar regions. The presence of a high 
pressure system over a source region enhances the development 
of the air mass. Source regions used in the classification 
of air masses are polar 'P', tropical ’T', continental 'C', 
and maritime 'M' . Table 5 describes the major air masses 
and their characteristics. Each of these air masses has a 
unique source region and movement pattern that varies 
seasonally. As would be expected, each also has a charac- 
teristic refractivity profile associated with it. Figure 37 
shows the average profiles associated with these air masses. 

The profiles shown in Figure 37 can be modified by 
the presence of subsidence associated with high pressure 
regions. Large, permanent high pressure regions exist over 
the subtropical and mid-latitude oceans of the world, to 
some extent year-round. These high pressure regions shift 
toward the poles during each hemisphere’s summer months and 
tend to drive the movement of frontal regions poleward. 

Frontal regions and cyclones, which are common in the mid- 
latitudes and polar regions, tend to destroy elevated 
layers, and, thus, no significant ducting or superrefraction 



112 



TABLE 5 

AIR MASS CLASSIFICATIONS AND GENERAL DESCRIPTIONS 



Air Mass Source region(s) Properties at Source 



Polar maritime 
(Pm) 


Oceans in latitudes 
greater than 50° 
(approx . ) 


Cool, rather moist; 
unstable 


Polar contin- 
ental (Pc) 


Continents in 
vicinity of Artie 
Circle; Antarctica 


Cold, dry; stable 


Arctic or 
Antarctic (A) 


The Arctic Basin 
and (central) 
Antarctica in 
winter 


Very cold and dry; 
very stable 


Tropical 
maritime (Tm) 


Sub-tropical oceans 


Warm; moist, and 
rather unstable 
near surface , dry 
and stable above 


Tropical 

continental 

(Tc) 


Deserts in low 
latitudes: pri- 
marily the Sahara 
and Australian 
deserts, but also 
south-west U.S.A. 
and Mexico in 
summer 


Hot and dry; 
unstable 
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Figure 37. Typical B-profiles in the source region 
for (A) Summer and (B) Winter. A is the 
ducting gradient. 
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is expected in the vicinity of fronts or of traveling low- 
pressure regions (cyclones). Due to the persistence of 
these traveling weather patterns in extreme northern or 
southern (polar) regions, anomalous propagation is generally 
a rare and intermittent occurrence at and near the poles. 

Coastal regions will reflect the variation in humidity 
and temperature related to the ebb and flow of maritime and 
continental air masses by exhibiting extreme horizontal and 
vertical refractive gradients. 

Over land, and to a lesser extent over water, there 
tends to be a distinct seasonal variation in the occurrence of 
anomalous propagation. This can generally be tied to the 
prevailing large-scale weather conditions and to the 
intensity of solar irradiance. 

Many of the smaller oceanic islands show little 
seasonal change in refractivity patterns, due to the fact 
that their climates are dominated by the maritime air masses. 
Islands large enough to affect these air masses will, on the 
other hand, often develop distinctly seasonal refractive 
gradient patterns. 

Continental regions, particularly those in the 
mid-latitudes, tend to show significant seasonal variation 
in refractive characteristics, with the greatest occurrence 
of superrefractive and ducting conditions occurring in the 
summer months and early fall. This is, in part, due to the 
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shifting of the subtropical high pressure regions, higher 
mean temperatures , and the greater amount of moisture 
present in air masses in summer. 

Seasonal shifts in the prevailing wind patterns, such 
as the monsoon circulations, lead to distinct seasonal 
variations in the refractive patterns of regions thus affected 

Many regions will also show strong diurnal variations 
in refractivity gradients due to surface changes in tempera- 
ture, humidity, and, in some regions, shifts in wind direction 
In general, according to Samson [Ref. 9: p. 9], continental 
and coastal regions tend to show anomalous propagation as 
occurring most often between 2000 and 0700 LST with the 
fewest occurrences between 1200 and 1500 LST. Oceanic 
islands show just the opposite trend in the occurrence of 
anomalous propagation. Subrefraction in desert regions also 
occurs most often in the late afternoon. As was mentioned 
earlier in this thesis, it is not uncommon in semi-arid and 
desert regions to see subrefractive conditions develop in the 
afternoon and superrefractive conditions form at night. 

Figure 38 shows this graphically, along with the resulting 
effects on IR rays. These diurnal characteristics, which 
tend to be strongest in the warmer summer and early fall 
months, are superimposed on the large scale changes in air 
masses . 

Forecasting anomalous propagation conditions involves 
many variables. However, studies of meteorological data 
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Desert Diurnal Refractive Characteristics and their 
effect on Infrared Propagation. 



collected from weather stations located worldwide has led to 
the development of distinct climatic types , which are 
terrain types with similar climates found across the earth’s 
surface. Table 6 is a list of climatic types and their 
associated meteorological and refractive characteristics. 

Meteorological data collected at weather stations 
has been used to prepare climatic charts reflecting the 
percentage of time superrefractive conditions or ducting 
occurs, on the average, across the surface of the earth. 
These charts provide some idea as to the probability of 
occurrence of superrefraction, surface ducting, and elevated 
ducting in any region of interest where sufficient data has 
been collected, and also reflect the seasonal variation in 
anomalous propagation over most of the earth’s surface. 
Figure 39 and 40 are representative of the types of graphic 
output climatic charts can present. These figures show the 
percent of occurrence of surface ducts, elevated ducts, and 
superrefractive layers. Note that Figure 39 is of Northern 
and Central Europe and Figure 40 is of the Mid-East. It is 
readily apparent that the percent of occurrence of surface 
layers and elevated layers varies a great deal between the 
two regions , and that coastal areas tend to show a greater 
occurrence of ducting than do inland regions. Charts like 
these, while helpful, do not show the percentage of time 
that subrefractive conditions occur, nor do they explain 
whether the anomalous propagation that they describe is 
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CHARACTERISTICS OF CLIMATIC TYPES 
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tions greater than 1000 
This group characterize 
dry climate. 



type annual annual characteristics 
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Figure 39A. Percent Occurrence of Surface layers for Northern 
and Central Europe. 
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Figure 39B. Percent occurrence of elevated layers for Northern 
and Central Europe. 



Q. UJ (E U LI Z K 




H 



o 

2 : 



-p 

0) CO 

rC nj 
-P <£ 

P 
0 
4h 

CO 

p 
a; 

>> 



a; 
a 
($ 

4h 
P 
3 
CO 

4h ^ 

o c 
o 

<D *H 

O bO 
G 0 ) 
(D P 
P 

P C 
3 aJ 
O 0) 

u c 

O rtf 
P 

-P P 
C CD 
(D -P 
O -ri 
P '■O 

0) a) 

CU S 



< 

o 

J- 

(D 

P 

3 

*H 

Pm 



123 



to include the Mid- 



CLUJQCULUZH 

I t i : : >15 5 >S5i > > s s > > • : 




Q_ LlJ QC U LaJ Z I— 



X 



124 



Figure 4 OB. Percent occurrence of elevated layers for the 
Mediterranean region, to include the Midr-east. 



intermittent or persistent. Since the measurements are made 
only a few times a day, the validity of the data is also 
somewhat suspect. Thus, it is difficult to judge the 
significance of the data collected on worldwide refractive 
conditions, although, at the very least, it does provide a 
general idea of regional characteristics. 

2 . Mirages and Turbulence 

The occurrence of mirages and atmospheric turbulence 
tends to be of a more localized nature than that of anomalous 
propagation. Collection of information on these two 
refractive conditions has been much less extensive than that 
done for anomalous propagation. 

Mirages are characteristic of large stretches of dry 
barren land or snow and of relatively large bodies of water 
that are located near or are surrounded by land. It should 
come as no surprise that the desert regions of the world are 
most commonly associated with mirages, particularly in view 
of the fact that inferior mirages are associated with 
subrefractive conditions. Prediction of mirage conditions 
requires knowledge of the same information as forecasting 
anomalous propagation but requires focusing down to much 
smaller regions on the earth’s surface. An additional 
problem in attempting to forecast mirage occurrence lies in 
the fact that mirages can be formed in both the absence and 
the presence of anomalous propagation. 
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Forecasting turbulent conditions is, perhaps, more 

difficult than other forms of refraction, since the 

interplay of surface roughness , solar irradiance and wind 

strength is difficult to predict. Likewise, the direct 

2 

measurement of is not always possible, however, Wesely 

and Alcaraz [Ref. 23: p. 9] state that an indirect method 

has been demonstrated. With this method, C 2 is calculated 

n 

using sensible and latent heat flux components of the 
surface energy budgets, which are available in meteorological 
literature. This method is most accurate for heights of 

4 meters or less, where comparison of values of determined 

2 

by the indirect method and by direct measurements of C 

n 

appear to agree closely, and thus provide a reasonably 
accurate means of predicting turbulence over land. Just as 
with mirage conditions and subrefraction, turbulence condi- 
tions occur most often, and with the greatest strength, in 
desert regions in the late afternoon, but, as seen earlier, 
significant turbulence does occur elsewhere and at 
different times of day. 

Although few actual climatic charts of turbulence 
conditions currently exist, over land solar heat flux charts 
and irradiance charts are available, and provide a fairly 
accurate idea of those regions where significant over-land 
turbulence can be expected. Figure 41 is an example of such 
a radiation chart. Note that, in general, the further away 
from the equator, the less solar radiation is received. 
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Figure 41. 



Solar radiation c^art 
measured in k j /cm . 



of eastern hemisphere 
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3 . Comparison of Mid-East and Central Europe 



The Mid-East and Central Europe are two distinctly 
different regions in terms of both terrain and climate. The 
Mid-East generally is typified by a hot, dry climate with 
relatively little seasonal variation in temperature and 
humidity except for its southeastern areas, which are 
affected by monsoons. Much of the interior is rocky desert 
with coastal mountain ranges rising up along the Mediterranean 
Sea. Coastal regions on theMediterranean Sea exhibit mild 
Mediterranean climates, while similar regions on the Arabian 
Sea are significantly affected both by continental air masses 
and by maritime air masses. 

Central Europe has a continental climate, characterized 
largely by cold winters and cool-to-warm summers. There is 
a moderate amount of precipitation year-round and relatively 
high humidity in the summer months. Coastal regions tend to 
experience some modification of the continental climate by 
maritime air masses. The terrain varies from relatively 
flat, open land in the north to wooded, hilly, and 
mountainous terrain to the south. 

It should come as no surprise, judging from the 
differences in climate and terrain, that the refractive 
characteristics of these two regions are distinctly 
different, as shown in Figure 39 and 40. The Mid-Eastern 
desert shows a significantly higher year-round incidence of 
both surface and elevated superrefractive conditions. 
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Figure 42A. Annual percent occurrence of elevated ducts 
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Figure 42B. Percent occurrence of elevated ducts during 
the month of maximum occurrence. This month 
of maximum occurrence varies from region to 
region. 
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Figure 42C. Annual percent occurrence of superrefractive 
layers (SRLRS). 
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Figure 4 2D. Percent occurrence of superrefractive layers 
( SRLRs ) during the month of maximum 
occurrence. This month of maximum occurrence 
varies from region to region. 
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formation is found in Europe, particularly over relatively 
large bodies of water and in the more arid south, but only 
during the warmer summer months. In the Mid-East, these 
conditions occur year-round and are fairly regular phenomena. 
By referring to Figure 41, it can be seen that the Mid-East 
receives a far higher level of solar irradiation annually 
than does Europe, and this provides the primary ingredient 
for the formation of subrefraction, mirages and optical 
turbulence . 

4 . Conclusions 

Although it can be shown, with some accuracy, just 
what percentage of the time superrefraction and ducting can 
be expected to occur worldwide, and, to a lesser extent, 
subrefraction, mirage formation, and optical turbulence, it 
is difficult to determine just where the occurrence of these 
refractive characteristics becomes "significant". Such 
decisions must be based on each specific weapons system (or 
component of a weapons system) affected, and on the nature 
of the effects, as well as on the expected probability of 
occurrence of anomalous propagation. Much research is being 
conducted by the U.S. Army at the Atmospheric Sciences 
Laboratory at White Sands Missile Range and at the U.S. Army 
Electronic Proving Grounds (EPG) at Fort Huachuca, Arizona. 
The research concerns the refractive effects for specific 
regions of the EM spectrum, and on specific weapon systems. 
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What can be determined by the results of this dis- 
cussion is that the various forms of atmospheric refraction 
are not rare, and may even be the norm in certain regions of 
the world at certain times of the year. Due to the wide 
regional variability in refractive forms and in probability 
of occurrence, it becomes apparent that each region of the 
earth poses a unique refractive problem. For this reason, 
when considering military operations in a specific region, 
the potential refractive conditions should not be approached 
in generalities, but, rather, should be addressed 
specifically to the region of interest. 
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III. ATMOSPHERIC REFRACTION AND GROUND WARFARE 



A. INTRODUCTION 

In Section II of this thesis, the three major forms of 
atmospheric refraction were discussed. The characteristics 
of each of these forms and their effects on EM waves 
propagating through the atmosphere were elaborated on at 
length in sub-Sections C through E. Those discussions 
showed the wide diversity of both the forms of refraction 
and their effects across the EM spectrum. Sub-Section F 
described the factors affecting the ability to predict and 
forecast atmospheric refractive conditions, contained a 
general discussion of the probability of occurrence of the 
various forms of atmospheric refraction across the earth's 
surface, and concluded with a description of the conditions 
characteristic of the Mid-East and Central Europe. 

The purpose of this background was to show that 
atmospheric refractive conditions which have a measureable 
effect on EM signals should not simply be considered as so 
rare or unpredictable that they can be "written off" as 
inconsequential. Rather, it has been shown that, at times, 
these conditions are quite common, can have a significant 
impact on EM signal propagation, and are reasonably 
predictable. Of course the probability of occurrence and 
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ability to forecast the three major forms of refraction 
varies from region to region (and often season to season) 
across the earth. 

With this general background as a foundation, the 
remainder of this thesis will now discuss specific examples 
of the effects of atmospheric refraction. This will be done 
in two ways - by use of computer simulations and by 
qualitative study. The examples will cover the full EH 
spectrum and are intended to enhance the reader’s under- 
standing of this subject as important for the employment of 
weapon systems utilizing the EM signals propagating through 
the atmosphere. 

B. EFFECTS OF ATMOSPHERIC REFRACTION ON TANK OPTICAL FIRE 

CONTROL SYSTEMS 

4 

This example is based on private communications with 
Dr. Don Snider of the Atmospheric Sciences Laboratory (ASL) 
in White Sands, New Mexico. This section discusses the 
effects of superref rac tive and subrefractive conditions on 
the accuracy of tank fire. The tank’s optical fire control 
system (using both visible light and IR energy) and current 



The communication with Dr. Snider concerned the effects 
of atmospheric refraction on tank optical fire control 
accuracy . 
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tank gunnery techniques (boresighting, synchronizing, 
zeroing of the main gun, and firing at center-of-mass ) are 
addressed . 

Current tank gunnery doctrine, as described in FM 17-12 
[Ref. 24: pp. 4-9 - 4-10], requires boresighting and 
synchronizing both sights and main gun, followed by zeroing 
the tank weapon systems prior to firing the tank's main gun 
and machine gun(s). Boresighting establishes a convergent 
relationship between the axis of the tube of the main gun 
and its direct fire sight. Synchronization insures that 
main gun and sights remain aligned when they are elevated 
and when they are depressed. The tank crew is expected to 
fire from bores ight using telescope and periscope during the 
day and, if so desired, the IR night vision sight at night. 

At the same time the crew must compensate for drift, cant, 
parallax, tube droop and winds on the battlefield. As would 
be expected, there are a great many potential sources of 
error involved in tank gunnery. The magnitude of the 
sources of error will vary depending on internal and 
external conditions (relative to the tank) and on the weapon 
and ammunition being fired. 

Figure 43 shows a typical tank fire control error budget. 
It is currently thought that the most important cause of 
error in firing the kinetic energy (KE) round is round-to- 
round variation and loss of zero or improper zero. For the 
slower HEAT round, crosswinds and air temperature extremes 
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Figure 43. Typical tank fire error budget. KE is the Kinetic 

Energy round and HEAT is the High Explosive Antitank 



are thought to cause the greatest errors in azimuth and 
elevation. Note that, originally, the figure shows optical 
ray bending causing only minor errors in elevation. 

It was the contention of the ASL report that the full 
effects of atmospheric refraction were not being considered 
in current tank gunnery doctrine, particularly those effects 
created in a hot/arid or semi-arid environment. The ASL 
proposed that refraction was affecting tank gunnery in two 
ways. First, atmospheric refraction could be inducing major 
elevation errors at ranges beyond 1500 meters. Second, 
those elevation errors currently blamed on other sources 
could actually be due to refraction. 

1 . Atmospheric Conditions Creating the Gunnery Errors 

As was previously described in Section II, atmospheric 
refraction patterns tend to follow a distinct diurnal pattern 
over land. During the day, the surface is heated up by the 
sun, and, in hot/arid or semi-arid regions this heating can 
cause an atmospheric temperature structure, such that the 
temperature decreases rapidly with elevation. This condition 
can cause optical turbulence, subrefraction, mirage forma- 
tion, or a combination of all three. At night, the earth’s 
surface re-radiates the heat absorbed during the day into 
the atmosphere, often leading to the creation of a 
temperature inversion. This inversion is the exact opposite 
of the temperature structure that characteristically forms 
in the late afternoon. About one-half hour before sunset, 
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and some time just after dawn, a condition is formed when 
the heat radiating from and into the ground is roughly 
equal . 

Local weather conditions, such as the presence of a 
high pressure region, can have a significant impact on 
atmospheric refraction. Depending on the weather conditions 
present, the diurnal characteristics described above can be 
affected in such a manner as to enhance or degrade either 
the subrefractive condition, the superrefractive condition, 
or both. Thus, atmospheric refractive conditions can vary 
significantly through the day, and from one day to the next. 

The ASL report was primarily concerned with the effect 
of going from a subrefractive condition during the day to a 
superrefractive condition (characterized by the temperature 
inversion) at night, then back again to the subrefractive 
condition. 

2 . Expected Effects on Tank Gunnery 

Since optical signals bend toward the denser medium, 
the direction of bending for the diurnal condition described 
above would be opposite for day (subrefraction) and night 
(superrefraction). This is graphically shown in Figure 38. 

If a target were located at a given range, its image would 
appear to shift upward at night and downward during the day. 
Figures 27A and 27B show such effects, created by a very 
strong temperature inversion (superior mirage) and an 
equally strong subrefractive condition (inferior mirage). 
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It should be noted that a mirage need not be present for the 
upward and downward shifting of an image to occur. 

Bending of optical rays can cause a gunner to aim at 
an image point which is high at night or low during the day, 
as shown in Figure 44. Since most tank crews boresight and 
zero their tanks during the day, the effects of subrefraction 
due to surface heating are largely removed. But this can 
result in a much greater error at night, as shown in Figure 
45. It should be reiterated here that there is little 
difference between the amount of refractive bending of 
visible light and of IR energy. Thus, the same errors occur 
going from an optical fire control system using visible 
light during the day to one that uses IR at night. 

3 . Experimental Results 

An experiment was conducted by the ASL to determine 
whether the theorized refractive effects do, in fact, occur. 
An American tank, using currently-deployed optical fire 
control equipment and current gunnery doctrine, was used in 
the experiment, which took place at White Sands Missile 
Range. An eight foot by sixteen foot target was used. The 
tank was boresighted, synchronized, and zeroed during the 
day, and compensation was made for as many external and 
internal sources of error as possible. Special effort was 
made to maintain the correct boresight while firing. If the 
boresight was off by .1 mils or more the gun was 
re-bore sighted . 
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dusk). The dashed tanks are the images actually 
seen and aimed at. 
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Once the main gun had been satisfactorily zeroed, a 
series of 5-round shot groups was fired over a two-day and 
three-night period. The results of these firings are shown 
in Figure 46. This figure clearly shows that the centers- 
of-impact of the nightfired groups were higher on the target 
than those of the day groups. The group shifts are measured 
relative to the center of impact of the zero confirmation 
group . 

The inversion and temperature lapse strengths can be 
deduced by the relative distance between the day and night 
5-round groups. As can be seen, the greatest distance 
between shot groups occurred between the first night and the 
first day, and the differences between the groups pro- 
gressively became less and less. On the first night the 
air was calm and clear, resulting in the creation of a 
strong superrefractive condition. The second night the sky 
was hazy, weakening the superrefractive condition. On the 
third night a front moved in, bringing winds and clouds 
which all but eliminated the superrefractive condition. The 
correlation between shot placement and weather conditions 
shows that refractive errors and tank fire accuracy are 
greatly affected by weather conditions. 

A statistical analysis was performed by the ASL, 
proving that the day shots and night shots were from two 
statically different populations. The results of this 



145 




S IIW ‘uoxsanoxg; uot^.'BA 0 X 3 



146 



Figure 46. 5-round group centers of impact of tank main 
gun fire. 



analysis showed, at a 95% confidence level, that the distance 
between the two shot centers was at least 46 cm. 

To bring the effect of the difference between shot 
groups into perspective. Figure 47 shows how many of the 
night shots would have missed the target if the gunner had 
actually been shooting at a silouhette of a T-62 tank while 
aiming at the base of the turret. This figure is based on 
the results of a vertical T-test of the 5-round shot groups, 
which determined the one sigma standard deviation in 
elevation and azimuth for each group. 

Although the above statistical analysis does show 
significant differences between day and night shot groups, 
it does not prove that the differences are due to refraction. 
To prove this relationship, temperature measurements at 
various altitudes above the ground were made concurrently 
with the firings. These measurements were made at .5 meter 
invervals from 0-to-4 meters above ground level (AGL) . 
Apparent elevation angle measurements using a theodolite 
were also made just prior to each firing. There measure- 
ments were made so that the line-of-sight of the theodolite 
was roughly parallel to the tank's, and was at the same 
height AGL. A regression analysis was then conducted, which 
showed a definite correlation between the temperature 
gradient, the theodolites observed absolute elevation 
variations, and the shot pattern elevation shifts during 
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Silouhette of T-S2 tank indicating that many night 
shots would have missed the target when using a 
daytime zero. 



firing. These interactions showed a significant relationship 
between refraction and firing errors, as shown in Figures 
48 and 49. 

In order to assess the importance of these refractive 
effects, based on the probability of occurrence of the 
necessary refractive conditions, the ASL investigated the 
frequency of occurrence of those atmospheric conditions 
which led to the most significant firing errors. Figures 50 
and 51 give a general idea of the frequency of occurrence of 
clear skies and light winds in southern Germany and in Iran, 
which would roughly indicate the occurrence of refractive 
conditions. Note that, for both cases, these conditions can 
be quite common at certain times of the year. 

4 . Conclusions 

The results of the ASL experiment showed that a 
definite and systematic day/night shift in elevation angle 
can occur for certain meteorological and surface conditions. 
The values measured agreed closely with what was predicted 
by theory. The magnitude of the measured effect was as 
great as .8 mr . At longer ranges even larger effects would 
be expected. Because the occurrence of temperature gradients 
most conducive to the greatest refractive effects is 
characteristic of hot/arid and semi-arid regions, it is 
expected that gunnery errors of this type will occur most 
frequently in the desert regions of the world. 
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graph compares well with a theoretical slope of .75. 
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knots for Tehran/Mehrabad Iran 
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This type of effect is a characteristic of the 
atmosphere and not of any particular sight design. The ASL 
is developing two refractive prediction models for use with 
tanks. The first model would use a micro-computer, require 
extensive input, and provide an accurate prediction of 
existing refractive conditions. The second model would be 
small enough to be placed within the fire control computer 
of the tank, and would only use input readily available to a 
tank commander. The output of this model would provide the 
tank commander with the guidance needed to correct for 
predicted refractive effects. Another idea would be to 
develop simple rules -of -thumb which could be used by a tank 
commander to predict refractive conditions and to provide 
recommended compensatory actions. 

Theory predicts a similar effect at MMW. If true, 
these effects must also be considered in any research 
currently under way into using this region of the EH 
spectrum. 

C. COMPUTER SIMULATION 1 - THE EFFECT OF REFRACTIVE 
CONDITIONS ON U.S. TACTICAL SYSTEMS, USING THE 
TROPOPLOT COMPUTER PROGRAM 

The purpose of this section is to explore atmospheric 
refractivity effects on actual army systems operating in the 
RF spectrum. A computer program called "Tropoplot" was 
used, which predicts long term median tropospheric 
transmission loss over irregular terrain. This computer 
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program was originally written in FORTRAN by the 
Environmental Science Services Administration (ESSA) research 
laboratories for Control Data Computers [Ref. 25]. The 
program was updated for use on IBM computers by LT James M. 
Callaghan [Ref. 26], A copy of the program is included in 
Appendix B. 

1 . Tropoplot Program Description 

The Tropoplot program uses a single surface refractivity 
value, and, although developed primarily to describe the 
effects of atmospheric refractivity on RF signals, calculates 
other effects, such as knife-edge diffraction. The surface 
refractivity value is related to a refractive gradient through 
a relationship described by Bean and Dutton [Ref. 27: p. 63]. 
However, the relationship does not provide gradient values 
which result in ducting and therefore the program is not used 
to evaluate systems operating in the presence of ducts. 

Other phenomena which affect communications systems, such as 
short and long term fading adjustments to median attenuations, 
are accounted for in the program by using factors taken from 
actual measurements made on various communications systems. 

All variables except the surface refractivity were held 
constant for the data runs in order to investigate only the 
effect of refractivity over the transmission path. 

An interdecile range is used to characterize the 
terrain. The ESSA report and the Callaghan thesis describe 
the interdecile height correspondence of various types of 
terrain. This relationship is shown below. 
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TYPE OF TERRAIN 



Very smooth plains or water 

Smooth plains 

Slightly rolling plains 

Rolling plains 

Hills 

Mountains 

Rugged mountains 

Extremely rugged mountains 



ASYMPTOTIC VALUE 
(meters ) 

OF THE 

INTERCEDILE HEIGHT 
0-5 
5-20 
20 - 40 
40 - 80 
80 - 150 
150 - 300 
300 - 700 
700 + 



The program is usable for frequencies from 20 MHz to 
40 GHz over distances of up to 2,000 kilometers. The surface 
refractivity can be varied from 250 to 400 N units. The 
ESSA research laboratories [Ref. 25: p. l] have made empirical 
measurements for comparison with the program results and have 
found them to be satisfactory. 

2. AN/FPS-16 Radar 

The first system used in the program was the AN/FPS-16 
air tracking radar. This high power system can be found at 
permanently fixed sites. The interdecile range for "smooth 
plains" was used. The transmission and reception line losses 
were not available, but were estimated to be 2 dB . Many of 
the system's parameters are provided in the signal strength 
data input section of Appendix B's data run output. 
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A "generic" radar warning receiver (RWR) with 
characteristics typical of fielded systems was assumed to be 
located on an aircraft at a radar range of approximately 125 
km and an altitude of 610 meters. This example was used due 
to its similiarity to the "standoff" type mission profile 
that might be flown by current army surveillance and signals 
intercept airborne platforms. 

During the data runs, the surface refractivity was 
varied from 250 N units to 400 N units in steps of 30 N units. 
The range at which the receiver could detect a signal varied 
from 117 km (at 250 N units) to 138 km (at 400 N units). 

This represented a fairly significant increase in the 
distance at which an RWR will detect the radar signal. 

The next step involved use of the IREPS program to 
compare the same radar in the presence of a strong duct and 
in standard atmosphere. Although IREPS is not intended for 
use over land, the comparison is helpful in exploring the 
similarity of conclusions for Tropoplot and IREPS. The 
outputs to the Tropoplot program is provided in Appendix B. 

The IREPS results, shown in Figures 52 and 53, cannot 
be meaningfully compared with the Tropoplot results for 
several reasons, including the fact that reflective 
coefficients are different for land and sea water. However, 
the loss diagram for IREPS in a standard atmosphere versus 
that for a surface based duct is interesting. It shows a 
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*»** LOSS DISPLAY **** 
FPS-16 TRKING RDR ( DET > 



location: NOT SPECIFIED 
DATE/TINE: STANDARD 




FREE SPACE DETECTION RANGE EASED ON 50>: POD OF A 1 SCI M AIR TARGET. 



DASHED LINE INDICATES DETECTION, COMMUNICATION, OR INTERCEPT THRESHOLD 



FREE SPACE RANGE: 163.3 NAUTICAL MILES 

FREQUENCY: 5325 MHZ 

TRANSMITTER /RADAR HEIGHT: 65.0 FEET 
RECEIVER/TARGET HEIGHT: 2000.0 FEET 



Figure 52. 



IREPS loss display for the FPS-16 tracking 
radar in a standard atmosphere. 
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**** LOSS DISPLAY **** 



FPS- 1 6 TRKING RDR < DET > 



LOCATION: 38 N 008E 
DATE-'TIME: SB DUCT 540 FT 




FREE SPACE DETECTION RANGE BASED ON 50M POD OF A 1 SO N AIR TARGET. 
DASHED LINE INDICATES DETECTION, COMMUNICATION, OR INTERCEPT THRESHOLD 



FREE SPACE RANGE: 163.8 NAUTICAL MILES 

FREQUENCY: 5825 MHZ 

TRANSMITTER /"RADAR HEIGHT: 65.0 FEET 
RECEIVER/TARGET HEIGHT: 2000.0 FEET 



Figure 53. IREPS loss display for the FPS-16 tracking 
radar operating in a surface based duct. 



loss which causes the radar’s minimum detectable signal 
threshold to be reached at about 74 nm in a standard 
atmosphere. This is actually the horizon limited value for 
the given heights of the transmitter (65 ft.) and target 
(2000 ft.). In a strong surface based duct, the path loss 
is much less and the threshold is reached at about 350 nm. 
The loss diagrams are included in Appendix B. 

3. AN/VRC-1 2 Radio 

The second system examined was a tactical radio with 
parameters similar to the AN/VRC-12 series used extensively 
by the U.S. Army for tactical VHF communications. The 
parameters were obtained from the U.S. Army Intelligence 
School, at Ft. Devens, Massachusetts, and from FM 24-24 
[Ref. 28: pp . 1-12] and FM 24-25 [Ref. 29: pp. 2-16], 
Parameters that were not available were estimated, and are 
considered to be close to those which actually exist. 

For this scenario, two radios are assumed to be 
operating over a fairly smooth surface with their antennas 
placed on terrain which allows for an antenna height of ten 
meters above the immediate land in order to enhance 
communications. The refractivity was again varied from 250 
to 400 N units, although the increments were increased to 
50 N units. 

In this example, the results were less impressive 
than in the previous case. The maximum range for reception 
only varied from 24 km to 26.1 km. Although not shown, a 
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simulation of a ground surveillance radar was run, using 
actual parameters, which also resulted in little or no 
change in the maximum detection range of 3000 meters, 
regardless of surface refractivity . 

The results of these latter computer runs suggest that 
refractivity may not be critically important for shorter 
range, lower power systems, a fact that is not totally 
un expected. 

4 . Conclusions 

Although it does appear that refractivity plays a 
critical role in the quantitative effects on RF signal 
propagation, further study, using a program for irregular 
terrain with strong vertical refractivity gradients , as 
opposed to a surface refractivity value, is necessary to 
define more precisely the critical ranges for refractivity. 

At present, there do not appear to be any programs which 
have this capability. 

D. COMPUTER SIMULATION 2 - THE EFFECT OF HORIZONTAL 
REFRACTIVE GRADIENTS ON AIRBORNE DIRECTION FINDING 
ACCURACY 

Airborne Direction Finding (DF) systems provide 
critically important support to U.S. Army tactical commanders. 
Systems which will provide enhanced accuracy are currently 
being planned or fielded. In general, accuracy is dependent 
on the system’s ability to measure individual lines of 
bearing (LOB) to the target emitters and on the precision 
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with which the navigation system can determine the airborne 
platform's location. Five major causes of error, which are 
independent of a DF system but which affect its ability to 
measure LOB’s are reflections, re-radiations, diffraction, 
refraction, and fading signals [Ref. 10: p. 2], The purpose 
of this section is to investigate the impact of refraction 
on DF accuracy through the use of computer simulation. 

It is important to recognize typical situations where 
refraction could cause problems. These include geographic 
areas where rapid changes in temperature and humidity (and 
sometimes pressure) occur in the horizontal plane. Such 
areas include broad river valleys, or large bodies of water 
surrounded, at least partially, by land, and coastal 
regions. In such areas the geometry involved in taking LOBs 
may vary greatly, depending on where the refractive gradient 
occurs in relationship to the emitters and the DF platform. 
This simulation used a fixed situation from which it is 
possible to extrapolate to other circumstances. This 
procedure is explained further in sub-Section 4. 

This section is organized into several sub-sections. 
Section 2 provides a brief general description of DF and of 
the tactical scenario used in the simulation. Section 3 is 
a basic description of the program, including the 
assumptions and algorithm. A copy of the program is 
included in Appendix C. The final section describes the 
simulation runs, the results of the simulations, and the 
conclusions . 
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It was necessary to reduce the amount of output from the 
simulation since a large number of LOBs , fixes, and other 
data was generated on each "flight" across the tactical 
area. The thesis contains examples of the pertinent data 
and graphic representation of the output in this section and 
in Appendix C. In addition, it would be possible, but not 
practical, to do a number of sensitivity analyses by changing 
the variables in different combinations. For example, the 
depth of the refractive area, the horizontal gradient, or 
the depth of the emitters could be changed, either singly or 
in combinations. However, the purpose here was to make 
representative changes in the variables from which it would 
be possible to draw definitive conclusions. 

1 . General Description of Direction Finding and the 

Tactical Scenario - 

Airborne DF systems operate, basically, like ground 
based systems, except that they must also determine the 
aircraft's location at all times. This is normally done 
using an inertial navigation system (INS) with some type of 
external update. The aircraft flies a preplanned course, 
often a "racetrack" pattern, while taking LOBs on target 
emitters. These are entered into emitter LOB files and are 
eventually correlated into fixes, a term used to describe 
the target's calculated location. The emitter files may be 
sorted by frequency or by other parameters. Most systems 
use further correlation or refinement to "throw out" some of 
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the fixes which do not fit the main distribution of fixes. 
Reports on emitters usually include an elliptical error 
probable (EEP) to allow for the expression of a confidence 
factor in the fix accuracy. The simulation computes an 
"average" fix based on all of the correlated fixes. 

A simple type of DF system uses a "two-bearing" cross 
(TBC) to produce fixes. In this system the point where two 
LOBs cross is a fix. The accuracy of the fix improves as 
the angle formed by the two LOB’s approaches ninety degrees. 
In this simulation, a TBC routine is used to produce fixes. 

The simulation used a corps level DF mission along a 
corps front of 200 kilometers. This is reasonable for a 
three division front. A single aircraft was used, although 
some systems use two or more in order to provide geometry 
that will allow almost instantaneous fix determination 
without having to fly long distances. The aircraft flew at 
a cruise speed of 200 knots and required a little over 30 
minutes to fly the front. To simplify the simulation, the 
aircraft flew the same baseline on each sweep of the corps 
area. 

Four emitters were used in the scenario and their 
distance from the corps FLOT could be varied (see Figure 54). 
The deeper locations are particularly applicable to the 
tactical concerns of a corps. Since emitters located close 
to the corps boundary would be difficult to accurately 
fix, due to poor geometry, the emitters in the simulation 
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Figure 54. Diagram of the Tactical Area 



were placed well within the corps lateral boundaries. Of 
course, in an actual situation, some compensation could be 
made for this geometry problem, such as an interface with DF 
assets of any adjacent corps. Three of the emitters (A, B, 
and C) are set up to evaluate errors from refraction and 
Emitter D is used to compare other types of error, such as 
DF system error or platform location reporting error, with 
refractive error. Emitter D, therefore, is not affected by 
refractivity . All of the emitters can move about within the 
200 kilometer area which corresponds to the aircraft 
baseline. The aircraft was assumed to be flying parallel to 
an area with a refractive gradient similar to that found 
along a coastline, taking fixes on land based emitters. The 
gradient area has a uniform horizontal refractive gradient 
(see Figure 55). The depth of this area can be varied in 
order to evaluate the effects on lines of bearing which are 
bent over long or short distances. Any refractive gradient 
value can be used, except zero. 

With the exception of Emitter D, an LOB departs the 
emitter in a straight line, passes through an area where it 
is refracted, and departs that area in a straight line 
which, like the previous straight line, is tangent to the 
arc. The LOB arrives at the baseline, where the angle of 
arrival and the position are recorded (see Figure 55). 
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Straight LOB 
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Figure 55. Geometry of lines of bearing from an emitter. 



The program does not examine the effects of vertical 
changes in refractivity since the ray tracing problem 
becomes extremely involved when working with both horizontal 
and vertical ray tracing. 

2 . Algorithm Description 

Throughout the program, a concerted effort was made 
to eliminate as many assumptions as possible, or, at the 
very least, to make assumptions which closely resembled real 
world circumstances. In addition, a graphic output was 
chosen over other output forms, such as statistical tables, 
as graphs are most effective in quickly evaluating where the 
fixes were computed. 

The following assumptions are pertinent to the 
understanding of the simulation: 

1. Flat earth geometry. 

2. LOBs are taken every 15 seconds and emitters are 
selected using a uniform distribution. 

3. LOB angles are taken with respect to the aircraft’s 
longitudinal axis. 

4. All emitters can be received at all times during the 
flight. 

5. The aircraft track is along the Y = 0 axis. The 
emitter depths assume a standoff range in order to provide 
aircraft survivability in a high-threat air defense 
environment . 

6. Emitter positions are fixed for each mission. 
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In reference to Emitter D, a set of preset flags 
controls the "on-off" condition of the variables. These are 
described in the first part of the program (see Appendix C). 
Basically, these flags control introduction of LOB error, 
navigation error (drift or random) , and navigation system 
update. During the majority of the simulation runs, LOB 
system error was used only with Emitter D. This allowed 
comparisons with the other Emitters (A, B, and C) which had 
refracted LOBs, but which were not subject to system error. 
Navigation error was not used as a variable in any of the 
runs, since earlier simulations showed it to be less of a 
factor than LOB system error. One simulation run was made 
in which Emitters A, B, and C were affected by both 
refraction and system error, as noted in the data (see 
sub-Section 3). 

During the simulation, the flags were turned on only 
when a comparison of system error versus refractivity error 
was desired. 

The algorithm which follows can be described quite 
briefly, even though the program itself is fairly large. An 
aircraft is started on a track with constant updating of its 
actual and reported positions. It selects an emitter, 
calculates the LOB to the emitter, and stores both the LOB 
and the position where the LOB was received. In the case of 
Emitter D, the actual and reported aircraft positions are 
recorded. After each LOB is taken, position on the track is 



169 



checked. If the aircraft is at a boundary of the corps 
area, a subroutine is called in to correlate the LOBs into 
fixes, calculate an "average" fix, and plot the output or 
summarize the output statistically. 

In order to validate the results, simulation runs 
were made using identical programs , except that variable 
values were reported after each step. A check was then made 
of selected data, both qualitative and quantitative. For 
example, the LOBs were checked qualitatively to ascertain 
whether their values approached that of a normal drawn from 
the baseline to the emitter, as the aircraft approached the 
normal. They were checked quantitatively by using a hand- 
held calculator to verify the angles. Many other checks 
were made, including a run with the dN/dX (refractive 
gradient) value close to zero in order to determine if any 
fix error was present. As a final check, large variations in 
gradient area and refractivity were introduced to insure they 
had the expected effects on the fix accuracy. 

Figure 55 illustrates the basic geometry used in the 
simulation. The LOB on the left would be placed into one 
file and that on the right into a second file. Later, arter 
the aircraft had reached the end of the baseline, the two 
files would be checked and the two LOBs would be correlated 
into a fix. The correlation routine checks every possible 
combination of LOBs which produce fixes. 
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The program can calculate a dN/dX value from inputs 
of temperature, pressure, and humidity for two different 
areas separated by a "gradient" area, as used in the 
simulation. In addition, a dN/dX value can be inserted into 
the program to facilitate sensitivity analysis. Different 
types of data can be produced by use of "write" statements 
located throughout the program. For example, the program 
will print out all of the LOB angles along with the locations 
where they were taken. 

An outline of the basic algorithm follows : 

Initialize variables. 

Calculate a dN/dX value. 

Use a random number routine to vary aircraft start on 
track . 

Establish a loop to control round trips of the corps front. 
Update actual aircraft position. 

Check to see if aircraft has completed left to right 
sweep; otherwise continue. 

Option for random nav. error, update reported position. 
Option for drift nav. error, update reported position. 
Option for LOB error. 

Select emitter from uniform distribution. 

Check aircraft position, if past point where a normal 
drawn from the emitter to the baseline would hit the 
baseline, set storage flag (Emitter D only). 

For Emitters A, B, and C use a slightly different routine 
to calculate the angles - one which depends on the aircraft 
location on the baseline. 

Calculate angle to selected emitter, in radians. 
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Convert angle to degrees. 

Calculate location where LOB was taken (Emitters A, 

B , and C ) . 

Check to see if the aircraft position is past the normal 
points (A, B, and C) . 

Call a subroutine to store the LOB, actual and reported 
aircraft position (Emitter D) . 

Call a subroutine to store the LOB and the position 
where the LOB was received (A, B, and C). 

Check to see if a right to left sweep is complete. 

If complete, call a correlation subroutine. 

Otherwise, continue flying and taking LOBs. 

Update aircraft position. 

Go to emitter selection. 

End of loop. 

Subroutine to store LOBs (4). 

Select one of two LOB files based on flag. 

Store LOB, actual position, and reported position in 
proper file (D) . 

Store LOB and position received in proper file (A, B, or C). 

Return to main program. 

Subroutine to correlate LOBs. 

Establish 2 loops to check every possible combination on 
LOBs in the storage files for each emitter. 

Establish a counter. 

Add 2 LOBs together and make a check for acceptable 
geometry. This is optional: for further information, 

see text. 

Determine vertex angle and convert to radians . 
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Calculate distance between both LOBs at baseline. 

Convert LOB angles to radians. 

Calculate length of one side of two-bearing cross from 
vertex to baseline. 

Calculate the Y coordinate. 

Add in Y nav. error (D). 

Calculate the X coordinate. 

Add the X coordinate to reported aircraft X position. 

Keep a running total of the X and Y coordinates for all 
fixes . 

End of loops. 

Call plot routine for fixes. 

Calculate average X, Y, and average fix error. 

Call plot routine for X, Y, and average error. 

Return to Main Program. 

3 . Results and Conclusions 

The graphs and tables in this section and in 
Appendix C, represent the data resulting from the multiple 
runs of this program. Since one graph for each emitter 
could be produced on each sweep of the corps area, it was 
necessary to limit the amount of output considerably. This 
was done through condition statements which allowed for 
sampling of data via collection from selected sweeps. Due 
to the uniform distributions used and the geometry, there 
were differences from sweep to sweep, as would be expected 
in the real world. Therefore, data used for comparison 
purposes were taken from the same sweep. 
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One of the first problems encountered concerned "angle 



gating" . This problem is related to the geometry of the 
two-bearing cross used to determine fixes. Ideally, an 
angle of ninety degrees formed by two LOBs (known as the 
vertex angle) would give the best fix. However, because of 
various factors, including width of baseline, emitter depth, 
emitter proximity to corps boundaries, and the necessity to 
have several fixes for better averaging, it was necessary to 
specify an angle gate which would provide acceptable 
geometry and yet allow for several fixes on each sweep for 
each emitter. This was done by using a hand-held calculator 
to determine the maximum and minimum possible vertex angles 
for each emitter at all ranges. Throughout the data runs, 
the accuracy was highly sensitive to the angle gates. This 
results from the fact that narrowing the gates eliminates 
some of the outliers in the data. 

The first simulation run checked the program by using 
a dN/dX value close to zero to indicate negligible fix error. 

The results showed that double precision was required in the 
program (due to division errors by numbers close to zero). 

The results are shown in data Table C-l in Appendix C. 

The initial simulation runs compared the effects of 
different amounts of refractivity on the fixes. The data 
is presented in Graphs 1, 2, and 3, on the following pages, 
and summarized in Tables C-2, C-3, C-4, and C-5, in Appendix C. 
The size of the gradient area was held at a constant 15 kilometers 
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during these runs. The output shows that the refractive 
index did have an effect on fix accuracy, but there were no 
cases where the error was as high as fixes produced with an 
LOB system error of 2.0 degrees (Emitter D). Even with 
dN/dX = -150.0, the errors for Emitters A, B, and C did not 
exceed 500 meters. Part of this is due to a "cancelling" 
phenomenon, which will be explained further. 

A gradient of -150 N units could not occur over a 15 
kilometer area since the range of N values would greatly 
exceed the norm. However, there is a case where the geometry 
could cause individual LOBs to travel for a fairly long 
distance through a strong gradient. This situation is 
illustrated in Figure 56. The individual LOB departs the 
emitter and travels in a direction which roughly parallels 
the refractive gradient. This results in the ray being bent 
over a longer distance, thus introducing greater error. 
Although extensive changes to the program would have been 
required to reproduce this particular geometry exactly, 
simulation runs were made with large gradients over areas of 
15 to 30 kilometers in order to investigate the effect of 
strong gradients and longer distances on LOB accuracy. The 
results are shown in Graphs 4 , 5 , and 6 , on the following 
pages, and are summarized in Tables C-6 and C-7 (Appendix C). 
As shown, this situation produced the greatest error from 
refractivity , particularly for the deepest Emitter (C), 
where the error exceeded 700 meters. The data appeared to 
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Figure 56 . 



Diagram of a Retracted Line of Bearing 
Using Alternate Geometry. 
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reflect some "cancelling” effect. This resulted when LOBs 
refracted in one direction were correlated with LOBs refracted 
rn the opposite direction, causing some lateral error to be 
cancelled out. The error, therefore, occurred mainly in 
range (Y-coordinate) and is shown in the fix data by the 
reported fix being deeper in range than the actual location. 

In order to investigate this further, a run was made where 
refracted LOBs from Emitter A were correlated with straight 
LOBs in a situation like that depicted in Figure 54 (i.e., 
refraction occurs over a large area) . An LOB system error of 
1.5 degrees was also introduced. The results are shown in 
Graphs 7 and 8 on the following pages, and Table C-8 of 
Appendix C. Note that Emitter A’s fix error greatly 
increased from previous runs, exceeding Emitter D’s error in 
fifty percent of the cases. 

The primary conclusion to be drawn from the simulation 
is that horizontal refractive gradients can be a problem under 
certain conditions, but that, generally, system errors will 
induce greater LOB inaccuracy than will refractive effects. 

The worst case is the problem shown in Figure 56, where the 
LOB is refracted over a long distance and the refracted LOBs 
are correlated with straight LOBs. However, if the 
possibility of these conditions exists, measures can be 
taken to minimize their impact. For example, the correlation 
software routine could ignore fixes produced by LOBs which 
came through an area of suspected strong refractivity , or a 
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correction factor could be introduced for refracted LOBs. 

Since good refractivity data would be very difficult to 
obtain in the emitter area, the software routine could employ 
a predictor system for LOBs to check for any that consistently 
arrive at the baseline at a different angle than would be 
expected. This would require close operator interface and an 
assumption that the baseline of the aircraft is long enough 
to allow sufficient LOBs to be taken which do not pass through 
an area of strong refraction. Although correction tables are 
sometimes used to correct DF errors arising from known 
refractive conditions , they are not a practical solution for 
mobile DF systems [Ref. 10: pp. 4-7]. 

This simulation pointed out several areas which could 
be studied further. These include a computer simulation which 
exactly reproduces the geometry of Figure 56 or one which 
looks at a time-of -arrival system working in an area of strong 
refractive gradients. The latter could involve some extensive 
work, since vertical gradients could not be considered 
independently. 3oth vertical and horizontal bending would 
affect the time of arrival. 

The magnitude of the effect of horizontal refractive 
gradients is dependent on many factors, including the 
emitter depth, the strength of the gradient, and the distance 
the LOB travels through the gradient. In many cases, the 
effects of inherent system errors will be greater than any 
refraction induced errors. Since refractivity errors are 
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often random and difficult to predict, attempts to correct 
or compensate for them may be hard to implement. Although 
improved software routines and other techniques may be 
effective in improving accuracy, the best procedure may 
be to recognize possible areas of strong refractivity and 
to use the fix results accordingly. 

E. SCENARIO 

The purpose of this Section is to show, by means of a 
general ground combat scenario, the large number and variety 
of U.S. Army systems that could be affected by atmospheric 
refractivity and the wide range of the resulting effects. 

1 . Scenario Background 

To provide an extreme but not unusual atmospheric 
refractive situation, a Mid-Eastern coastal scenario will be 
used. A mechanized division, with the necessary corps 
supporting assets, has deployed to this region during a 
period of increased tension, but before any actual combat 
has occurred. 

The mechanized division has three brigades on line, 
a tank battalion held in reserve, and a strong covering 
force forward of the front line of troops (FLOT). Since 
this description is not meant to be a lesson in ground 
combat operations, little further effort will be spent on 
describing the disposition of or organization for combat of 
the elements of the division. 
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The division has one flank anchored in a coastal 
mountain range, adjacent to another U.S. division, and the 
other flank is anchored in the sea. The division is currently 
located in a blocking position across a major high speed 
avenue of approach, as shown in Figure 57. 

The division is deployed in an arid, desert region, 
characterized by large areas of sand and rock of generally 
good traf f icability , interspersed with impassable semi- 
forested mountains, sand dunes, and ravines. Very little 
vegetation exists in the region except on the upper slopes 
of the coastal mountain range. The temperatures at this 
time of year range from the 70's at night to 110 degrees 
during the day. Days and nights are often clear, but, 
vision can quickly be obscured by violent sandstorms. 

In this scenario, two days and one night will be 
considered. On the first day the extreme daytime tempera- 
ture leads to the creation of a strong subrefractive 
condition over the coastal strip which is strengthened by an 
inferior mirage from the coastal range to the sea. The high 
temperatures and onshore breeze also lead to high 
turbulence, throughout the day, along the coastal strip. 

That night, a high pressure area begins moving 
slowly into the region, enhancing the superrefractive 
condition created by the radiation of heat from the desert 
surface. Some turbulence exists but not nearly so severe as 
that which existed during the day. 
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Figure 57. Example Scenario. 
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During the second day, the high pressure area leads 
to the creation of an offshore flow of dry, hot air and 
results in the creation of a strong elevated layer and 
surface duct that extends inland almost to the coastal range 
2 . Atmospheric Refractive Effects 
a. Optical Effects 

The presence of the subrefractive condition, 
inferior mirage, and optical turbulence on the first day 
will significantly affect many of the division’s optical 
surveillance, fire control, and tracking systems. 

The distance at which the mirage begins will 
limit visual observation, by ground observers, of areas to 
the front of the division’s covering force. Mirages can 
significantly affect optical systems as close as 1000 meters 
In any combat, this near boundary of the mirage can reduce 
the detection and engagement ranges and accuracy of 
direct-fire weapons, particularly if optical turbulence is 
also severe. In a like manner, the presence of the mirage 
condition and optical turbulence could also reduce the 
effectiveness and range of laser rangefinders and laser 
designators . 

Some relief from these adverse refractive 
effects can be gained by placing observers as high above the 
surface of the coastal strip as possible. Those observers 
located above the top of the mirage, either on the 
surrounding high ground or in aircraft, will be unaffected 
by the mirage. 
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Another adverse refractive effect would result 
if the division 1 s tank main guns are boresighted and zeroed 
during the first day and fired the first night or second 
day. As explained in Section III.B., the result of such a 
situation could be significant reduction in tank main gun 
accuracy during night firing and during the second day. 

Some bending of IR paths will occur at night due 
to the superrefractive condition, but the effect of this on 
FLIR systems would probably be insignificant. 

The ducting on the second day would result in 
a weakening or total absence of the inferior mirage seen on 
the first day. Thus, visual observation and laser designation 
and rangefinding systems would be more effective on the 
second day, although optical turbulence might still be quite 
strong . 

b. Effects on Communications and Electronic Warfare 

Systems 

Desert regions, in general, have an adverse 
effect on RF communications. Some of the effects include 
difficulty in adequately grounding communications systems, 
existence of a weak HF ground wave, and a generally low 
surface reflectivity, all of which lead to a significant 
shortening of HF, VHF, and UHF communication ranges. 

The presence of the subrefractive condition on 
the first day, although causing some bending of ground-to- 
ground communications signals and ECM jamming signals, will 
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not result in any additional significant loss of range. This 
is not necessarily the case for air-to-ground communications 
and jamming. The existence of the subrefractive layer over 
the surface of the coastal strip will set a limit on the 
minimum angle, from the aircraft to the ground station, at 
which communications will be dependable between the two 
stations . This latter effect will also place some limit on 
friendly airborne communications intercept and direction 
finding CDF) capabilities for platforms located behind the 
PLOT and attempting to acquire and DF hostile communications 
emitters. The height of the aircraft above ground, the 
distance to the ground communications system, the power 
output of the airborne system, and the strength, height, and 
extent of the subrefractive condition are all critical in 
determining the overall effect. 

The accuracy of airborne and ground based DF 
systems can also suffer from a horizontal atmospheric 
refractive gradient created by the juxtaposition of land and 
sea air masses along the coastal zone. These errors can be 
quite significant, as shown in Section III.D. 

The night time superrefractive condition and the 
ducting occurring during the second day, both of which are 
surface based, will tend to extend communication ranges of 
those signals operating within these conditions and above or 
near the minimum trapping frequency, f . This will serve to 
increase the distance between ground based units , for which 
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communications will be considered reliable, and also the 
distance from j ammer-to-target at which jamming will be 
effective. Communication intercepts of distant hostile 
emitters will similiarly be enhanced. In a like manner, 
hostile Electronic Warfare (EW) units will be able to jam, 
intercept , and DF friendly communications from greater 
distances. Terrain features may, however, override the 
effect of the superrefractive condition, and, in any event 
the actual enhancement may be only a few kilometers for 
ground based systems, as shown in Section III.C. 

c. Effects on Noncommunications Emitters 

The division's noncommunications emitters, 
primarily radar, are spread throughout the division and 
serve a variety of purposes. 

The accuracy of the division's counterbattery/ 
countermortar radars may be degraded by the subrefractive , 
superrefractive, and ducting conditions. The extent and 
nature of such atmospheric refraction produced errors in 
the Firefinder radar is currently a subject of interest. 

Ground surveillance radar systems will suffer 
very little from atmospheric refractive effects, since these 
systems are operated over relatively short ranges and are 
directed along the surface of the earth. 

Ground based air defense radars, however, can be 
adversely affected by all forms of atmospheric refraction. 
The first day's subrefractive condition will result in 
errors in both angle and range, but will be unlikely to 
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affect early warning or tracking radars significantly. The 
first night’s superrefractive condition will cause radar 
signals to bend in the opposite direction from those 
propagating during the day, and will thus cause angle and 
range errors in air defense radars; but, again, because of 
the nature of the radar’s purpose, these errors will 
probably be insignificant. The strong ducting occurring on 
the second day will result not only in larger angle and 
range errors than those occurring at night, but will also, 
in all probability, result in the development of radar 
holes. These radar holes could cause reduced detection 
ranges for early warning radars. 

Airborne radars, such as the side-looking 
airborne radar (SLAR), will suffer a reduction in range as a 
result of the subrefractive conditions occurring on the 
first day, if the aircraft is flying above the condition. 

The superrefractive and ducting conditions of the first 
night and second day may result in an extended range for the 
SLAR system, but this would depend on the height at which 
the aircraft is flying and the height of these refractive 
conditions . 

3 . Cone lusions 

It can be seen that a multitude of combat systems as 
well as surveillance and intelligence gathering systems are 
affected by atmospheric refractive conditions. An ability 
to predict the existence of such conditions can go a long 
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way towards explaining irregular or "unusual" effects on 
army ground systems. The ability to forecast atmospheric 
refractive conditions will also provide an opportunity to 
take advantage of those effects that can enhance certain 
capabilities and to compensate for those that will degrade 
other capabilities. 
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IV, CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

It has been the purpose of this thesis to highlight the 
e^.ects of atmospheric refraction on ground force* s systems 
that utilize EM waves propagating through the atmosphere. 

To this end, the various forms of atmospheric refractivity 
in the troposphere have been described, and their 
characteristic effects on EM waves have been examined. 
Additionally, specific examples and simulations were used to 
expand and to reinforce some of the points brought out 
earlier in the thesis. 

From the discussions, examples, and simulations presented 
herein, a number of points should be apparent. First, 
atmospheric refraction in several related but different 
forms does affect the propagation of EM waves traveling 
through the atmosphere. Second, the specific effects of 
atmospheric refraction on the different regions of the EM 
spectrum will vary, not only with wavelength, but with the 
form of the refractive condition. Third, to some degree the 
forms of atmospheric refraction can be forecast, given 
certain meteorological parameters. In a like manner, 
knowledge of the earth region of interest can be used to 
predict what form or forms of atmospheric refraction to 
expect. Fourth, specific military systems using the EM 
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spectrum are affected by refractive conditions in the 
atmosphere, and these effects result in the associated 
weapons system’s performance also being affected. Fifth, 
the significance of atmospheric refraction on any given 
system is dependent on the form of the refraction, the 
strength of the effect, the form and purpose of the affected 
system, and the tactical scenario. There is, therefore, no 
easy answer to the question of whether atmospheric refraction 
is significant to a system’s performance or whether it occurs 
often enough to be considered a significant effect. This 
latter consideration is dependent on the system affected and 
on the region of interest. Sixth, tactics can be developed 
to compensate for or to utilize certain atmospheric 
refractive effects. Finally, it should be apparent that 
atmospheric refractive effects on ground warfare should be 
given consideration by the tactical commander. 

The ability to forecast, predict, and understand the 
effects of atmospheric refraction allows the commander to 
avoid being surprised by many of the "unexpected" and 
"unexplained" problems that occur in systems using EM 
radiation. This ability also provides the commander with 
opportunities to use certain refractive effects to his own 
advantage . 

The tank gunnery results and the simulations in this 
thesis indicate that the U.S. Army’s primary problems with 
atmospheric refraction will occur in the electro-otpics 
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spectrum. There are some cases in the RF spectrum where 
refractive effects are fairly significant, as when lines of 
bearing travel ror long distances through a refractive 
gradient. However, the effect on electro-optics (and 
possibly millimeter waves) over relatively short distances 
is more significant. Due to the shorter operational ranges 
of most EO systems, as compared with RF systems, the 
possibility of detecting and correcting refractive problems 
in this region may also be greater. 

In short, it is hoped this thesis has shown atmospheric 
refraction to be worth further consideration and study. In 
the next sub-section a list of recommendatiors for the 
direction of further studies, is provided. 

B. RECOMMENDATIONS 

Following are recommendations for further study of the 
effects of atmospheric refraction on ground warfare, and 
additional recommendations as to when and where this 
subject should be considered. 

1. Further study of the wavelength dependence of 
refractive effects is needed. 

2. Atmospheric refraction should be a consideration in 
system design, particularly in system testing. 

3 . Study is needed concerning atmospheric refractive 
effects on specific systems (e.g. AH64 FLIR, Firefinder 
Radar, JSTARS, etc.). 
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4. As a result of studying the effects of atmospheric 
refraction on specific systems, methods should be developed 
to compensate for and/or take advantage of these effects. 

5. An effort should be made to relate the effects of 
dispersion, extinction, reflection (multipath), and 
diffraction into the overall picture of atmospheric 
propagation . 

6. Further study is needed on determining the causes of 
atmospheric refraction in its many forms , and into those 
meteorological parameters needed to forecast atmospheric 
refractive conditions accurately. 

7 . A greater understanding of the reflection 
coefficients for various vegetation and terrain types and 
their impacts on EM propagation is needed before a 
prediction model can be developed. 

8. Development of a systems prediction model along the 
lines of the Navy’s IREPS model is a must. 

9. Once a prediction model is under development, a 
study will be needed to determine who, within the Army's 
organization, will measure the meteorological parameters 
needed to predict atmospheric refractive effects. 
Determinations must also be made as to whom the responsibility 
for operating the prediction model will be given, to whom 

the data will be disseminated, and to whom the data will be 
disseminated . 
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APPENDIX A 



IREPS 

The U.S. Navy has been aware of the effects of the 
atmosphere, particularly of atmospheric refractivity , on EM 
waves since at least World War II. However, until recently 
the exact nature of these atmospheric refractive effects , 
and therefore an ability to accurately predict them, has 
been non-existent. 

With the introduction of an Integrated Refractive Effects 
Prediction System (IREPS) by the Naval Ocean Systems Center 
(NOSC) of San Diego, the Navy has now developed the 
capability to predict atmospheric refractive conditions at 
sea. IREPS uses a desk-top programmable computer (HP 9845) 
with a graphic display to provide evaluations of the effects 
of refraction on surveillance, communication, electronic 
warfare, and weapon systems. This system is currently 
being installed on all Navy aircraft carriers, to be used as 
an aid in developing tactics to assess, compensate for, and 
take advantage of refractive effects at sea. Eventually, 
other atmospheric effects will also be integrated into the 
IREPS program, along with the PROCAL-PREOS airborne FLIR 
performance model, in order to develop a more comprehensive 
prediction system for future use. 
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It is important to note that the effects of concern in 



IREPS are only significant at EH frequencies above 100 MHz, 
are considered only for the lower atmosphere ( troposnhere) , 
and only over water. 

1 . Products of IREPS Model 

There are eight basic products that can be generated 
by IREPS, according to the IREPS Revision 2.0 User's Manual 
[Ref. 30: p. 30], Each product is dependent on a variety of 
data that must be entered into the program. For the first 
product, only the latitude and longitude are needed, but for 
the remaining seven products, a variety of on site environ- 
mental data are required. Each product is displayed on the 
minicomputer screen, and can be printed on an 8 1/ 2-by-ll-inch 
page. The eight products are: 

a. The Historical Propagation Conditions Summary 
Figure 58 is an example of the Historical 

Propagation Conditions Summary. This product provides a 
description of the refractive average conditions for the 
geographic location specified by an input longitude and 
latitude. A refractive profile can also be generated for 
this region by season and by conditions of darkness. 

b. The Environmental Data List 

Figure 59 is an example of the Environmental 
Data List, used primarily for checking values of radiosonde 
data entered into the program, such as pressure, temperature, 
and relative humidity. Values that are produced and printed 
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Figure 58. Historical Propagation Conditions 
Summary Product. 
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IREPS REV 2.3 

**** ENV I RONMENT AL DATA LIST 

LOCATION: 31 56N 1 1 S 36W 

DATE/TIME: 17 JUN 0045Z 

WIND SPEED 12.0 KNOTS EVAPORATION DUCT PARAMETERS: 

SEA TEMPERATURE 18.2 DEGREES C 
AIR TEMPERATURE 15.1 DEGREES C 
RELATIVE HUMIDITY 89 PERCENT 
EVAPORATION DUCT HEIGHT 28.0 FEET 

SURFACE PRESSURE » 1008. O 
RADIOSONDE LAUNCH HEIGHT » 60.0 FEET 
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SURFACE REFRACTIVITY: 341 --SET SPS-48 TO 344 



Figure 59* Environmental Data List Product* 
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out in the Environmental Data List are the dewpoint 
depression, altitude, N unit, N unit gradient, M unit and 
a description of the refractive condition. 

c. The Propagation Conditions Summary 

Figure 60 is an example of the Propagation 

Conditions Summary, showing the existence of refractive 
conditions for a specified region at a given date and time. 
The summary provides a plot versus altitude diagram of N 
units and M units. The extent of any duct present is shown 
by the use of shaded areas in the block on the right hand 
side of the summary. A brief statement on the performance 
of several types of EM systems operating in the specified 
region is provided at the bottom of this summary. 

d. Surface-Search Radar Range Tables 

Figure 61 is the Surface-Search Radar Table 
format. This product provides detection range predictions 
for naval surface-search radars. An actual example of the 
output is not provided since the predicted range values are 
classified CONFIDENTIAL. 

e. The Coverage Display 

Figures 19 and 20 are examples of the Coverage 
Display product. This product shows the area of coverage of 
a specific type of radar or radio emitter on a curved earth 
range-versus-height plot. The shaded area is the area of 
detection or communication for the system represented. The 
impact of atmospheric refraction on a system can be seen 
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F> OFAGAT I 0 r 4 C ON D I T I OH *. SUMMARY 



LOCATION: 31 56N 113 36U 

DATE TIME: 1? JUN 0045Z 




WIND SPEED= 12.0 KNOTS 
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: UP FAC E REFRACT I VI TV: 341 --SET SPS-43 TO 344 



Figure 60. Propagation Conditions Summary Product. 
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LOCATION: 

TINE: 



***** SURFACE SEARCH RADAR RANGE TABLE ***** 



SURFACE SEARCH RADAR: 

RADAR ANTENNA HEIGHT: FEET 



U S SHIP TYPE'CLRSS 


DETECT] 

MIN 


rUTT RANGE 
AVG 


nH'HH 

MAX 


CVF'CVN 

CG-'CGN 

DD-'DDG 








FP-'FTTT 

LCC 

LHR 








LPH 

LKA 

LPD 








LSD 
LST 
AE- AF 

AO - AO E-' AO R 









SOVIET SHIP TYPE/CLASS 


DET ECT ] 

m 1 n 


OHRHNbE 

AVG 


■ IN NH 
MAX 


KIEV CLASS 
MOSKVA CLASS 
CLG 








CG -cC - CA 

DD-DDG 

FRIGATE 








CORVETTE 

OSA-STENKA CLASS 
PRIMORYE CLASS AG I 








LENTRA CLASS HOI 
OKEAN CLASS AG I 









EVAPORATION DUCT HE I GHT = METRES. 

= FEET 



Figure 61. Surface Search Radar Range Table Format. 
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graphically by comparing printouts of the radar coverage 
with and without refractive effects present. Note also that 
the presence of a radar hole can be seen in this product by 
the shape of the coverage pattern above the duct. 

f. The Loss Display 

Figure 62 is an example of a Loss Display 
showing one-way path loss in dB versus range. The dashed 
line on the printout is the threshold for detection, 
communications, or intercept. The solid line is the 
system's signal. Thus, a determination of where the signal 
would be strong enough to be detected can be estimated from 
this product. 

g. The Radiosonde Observation Analysis 

This product is provided largely to assist 
meteorological personnel, and consists of a radiosonde 
listing, a skew T, a log P thermodynamic diagram, and an 
encoded upper air observation message. 

h. The ESM Intercept Range Table 

Figure 63 is the format for the ESM Intercept 
Range Table. This table provides the maximum intercept 
range for several different Naval ESM systems versus a 
variety of hostile emitters. An actual example is not 
provided since much of the information on the hostile 
emitters and on the predicted values is classified. 
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***» LOSS DISPLAY *»■** 
SURFACE UHF COM 



LOCATION: 31 56N US 36W 
DATE/TIME: 17 JUN 0045Z 




BASED ON A 200 NM FREE SPACE COMMUNICATIONS CAPABILITY 
TO A SURFACE SHIP PLATFORM 

DASHED LINE INDICATES DETECTION, COMMUNICATION, OR INTERCEPT THRESHOLD 



FREE SPACE RANGE: 200.0 NAUTICAL MILES 
FREQUENCY: 300 MHZ 

TRANSMITTER /RADAR HEIGHT: 13O.0 FEET 

RECEIVER/TARGET HEIGHT: 100.0 FEET 



Figure 62. Loss Display Product. 
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LOCATION 

TIME: 



***** ESM INTERCEPT RANGE TABLE ***** 



ESM RECEIVER: 



EMITTER CLASS: 



EMITTER 


— Frrcr ' 

CM Hz> 


TTff?n WTEBCEPT 
RANGE (nmi > 


KWIFE REST 
KNIFE REST B 
CROSS BIRD 






squhrE hEAD 
HIGH POLE 
FAN SONG E MG 






Top trough 
BIG NET 
TOP SAIL 






With LUNE 
SCOOP PAIR 
HEAD NET 






slimHet 
LOW SIEVE 
BALL END 






"HICH'SIEVE 
FRONT DOOR 
TRAP DOOR 






STRUT PAIR 
STRUT CURVE 
FAN SONG E MT 







WLR-1 CV 



SOVIET 



EMITTER 


FREO 
< MHz ) 


MAX iNtErcETPT 
RANGE (nmi> 


HEADLIGHT 
MUFF COB 
POP GROUP 






DA5S TlLt 
DRUM TILT 
OWL SCREECH 






SQUARE TIE 
SNOOP TRAY 
PEEL GROUP 






HHWK SCREECH 
TOP BOW 
SNOOP PLATE 






DONETS 

DONETS-2 
POT HEAD 






LOU THOUGH 
SUN VISOR 
NEPTUNE 






normt? 

DON/DON-2 







EMITTER CLASS: 



EMITTER 


— FPECT 

< MHz ) 


"WBK" INTERCEPT 
RANGE <n mi ) 


SPS-43A 

SPS-29 

SPS-37 






SPS-37H 

SPS-32 

SPS-40 






“SP5-49 

IFF INT 
TACAN 






SPS-39 

SPS-42 

SPS-48 






SPS-32 

MK-26 

SPS-39A 






MK-33/M0D 0 

SPS-33 

SPS-30 






SPW-43 

SPN-S 

SPS-10 






SPG-49 RCQ 

SPG-51 

MK*-37 






tfFS-5,9, 11-15 

MK-13 

MK-34 







U.S. NAVY 



EMITTER 


FREQ ’ 

CMHz> 


'MAX' "INTERCEPT " 
RANGE (nmi ) 


SPC-33A 

MK-68 

SPG-34 






5PG-T0 
SPQ-9A 
MK-23/M0D 3 






HH-35^H0D 2 
MK-36 

MK-23/M0D 2 






MK-&7 

SPN-35 

SPS-46 






5P“S^3 
CPR 1300 
CPR 2900 






LH 5b 

RAYTHEON 2302 
RAYTHEON 2840 






RAYTHEON 1 900 
DECCA 202 
DECCA 914 






HEL-h "18/9 

SPS-33 

SPN-12 






— MK-1 13 

SPG-33B 

SPN-41 







EVAPORATION DUCT HEIGHT- METRES 

- FEET 



Figure 63. ESM Intercept Range Table Format. 
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2 . Limitations of the IREPS Model 



There are a number of limitations in the IREPS 
model, some of which are being corrected. These limitations 
include : 

a. Limited frequency range. 

b. Neglect of the effect of land and sea clutter. 

c. Assumption of horizontal homogeneity in the 
refractive structure. 

d. Limits on the maximum and minimum ground 
based antenna height. 

e. Exclusion of the effects of sea-reflection 
interference on airborne systems. 

f. Assumption that all EM systems are horizontally 

polarized . 

g. Absorption is not accounted for. 

h. Limits exists on the consideration of low 
level elevated ducts. 

i. The coverage diagram is not applicable to all 
radar applications. 

3 . Applicability of IREPS to Over-land Refractive 
Prediction ~ 

While IREPS is proving to be of considerable value 
to the Navy even with its wide range of limitation, its 
applicability over water does not imply that it will work 
equally well over land. IREPS in its present form, because 
it was designed for over-water use, does not account for 
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terrain features, vegetative effects, soil type, or 
variation in the surface reflective coefficient, all of which 
have great effect over land. In addition, the nature of 
atmospheric refractivity over land is much more transient 
and localized than over water, and this must also be taken 
into account in developing an over-land model. The basic 
design of IREPS, however, may prove to be viable in the 
development of an over-land model and, at the very least, is 
worthy of further study in this regard. 
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APPENDIX B 



COMPUTER SIMULATION 1 LISTING AND OUTPUT 

This Appendix contains the program and output for 
Simulation 1. Each computer run is summarized in four 
pages: Summary of the Parameters, a Transmission Loss Table, 

A Loss Diagram, and a Signal Strength Table. The latter 
contains the refractivity value. The type of system (either 
radar or communications) can be identified by the frequency 
at the top of the first output page for each run. 
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WRITE "DISTANCE, TRANS-LOSS" HEADING AND DATA 
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WRITE "DISTANCE AND SIGNAL STRENGTH" HEADING AND DATA 



cOX 

too 



ZUJ 






•k 




UJX 






<M 




>M 






• 




MCO 






r* 




O 1 






o 




O 






UL •* 




X** 






••CM 




a co 






v • 




X» 






r^ 




•» 






U.O 




IUX 






OH 




0<M 






* 




z - 




in 






<- 




* 


ujm 




MlU 




m 


>o0 


Nk 


COO 




it 


LU< 


X. 


mz 




z 


OLU 


N* 


Q< 




* 




* 


H- 




— 


H- 


« 


cOcO 




-* 


3h 


1 


33 




M 

1 


X< 

Z 


1 


X- 




z 


MlU 





tu •* 
>x 
m 



* 

— o 



K I 



CM 



I" '"2 



>*— V.UJV. 
•*H» *< * 
</)• 



CO 

x 

LUZ 



00 

X 

LU 



-J< 



— Q LU 
UJ *-*COO 
I X H- cOh-c/O X< 

00 | LU LU-»^w 2TO01UQ. 

| — -I X-»XOO^ <<H -00 
S<* 1 t/)S 22 ?tOMZ LU 
* • •ULO— 'C5— OLUOZXUJ 

— • oOwuj x I -j<lu i ujcn 



X^ 


X 


CM— 


>< 






M 


Oh- 


-j^j-j x>x>zoxv-xxa 


i-m 


M- 




mO 






UJ 


xz 


a<Q-J M K z z X 


o •* 


* 


X - 


LU 






X 


— o 


xo **<X M »-• < < 


zx 


m 




000 






< 




• • 12 • 


ium 




cOX 


LU< 






X 




hQO i *>Z "Z •‘Q •» 


x - 


r 


Q 43 


XZ 






< 


ax 


Zo£ *m m m a lu «j 




c 0 


•k «k 


z 






X 


hO 


a LULU 00 O O — • UJ LU 


cO>v 


a 


—04 


<LU 








H-X 


MX (\ 43 in<C\JH-in 304 CQ 04 


Vk 


'W 


— • 


M 






X 


< 


1 — I * 41 * «Q %c£ •<X. • • 


JS 




MOO 


XZ 






< 


NJ 1 


< 1 h-ON 2 MLN>MUN 2 N 


<s 


• 


1 LL 


o< 






X 


1 


1^4 I 00 lLGU.lua.liJlLCOU.Oli. 


zs 


X 


z - 


LL 




O' 


CJ> 


X 1 


M | M>» »X "J •* MH •* 


os 


ao 




z 


00 


43 


OH 


< 1 


X I QU* • • UJ* UJ* M* 


HS 


•k 


* 1 


X*LU 


X 




X 4 D 


-J 1 


<1 *Z i uj i uj | lOKl 


00 ^ 


• 


O 1 


<*UJ 


UJ 


a 


X 


a i 


-J 1 UJUJ lOIOULlZlGl 






ni | 


X# 35 


l~ 


h- 


o 


X 1 


o i H-a i<i<iai<izi 


1 •• 


X 


+ 1 


Q#H 


UJ 




1 h- 


- 1 


xi -aiu.iu .1 ihiuji 


1 • 


X 


X- 


** LU 


X 


O 


1 


- 1 


- 1 C 0 UJ 1 X 1 X 1 X 1 00 1 h- 1 


1 1 




-k-k * 


o- to 


< 


O 


i a 


X 1 


- | ZX IGIGIGlHihl 


- 1 


* 


Q <f 


O - 


X 




• o 


m i 


XI -LL| 00 | 00 | 00 lOt<l 


- 1 


* 


00 • 


XXUJ 


< 


— 


•k 


M | 


(A| 00 - | - | . | . 1 - | - | 


X 


X 


W GO 


xoo 


X 


in 


X— 


-1 


M I * 1 1 1 1 1 1 


moo 


— mo a. 


-z 




• 


040 


V* | 


HU. 1 | | 1 | | 


cox 




-*N*<— 


h- 


CJ 


in • 


>> 1 


^1 1 f 1 1 1 1 


> HU 


> 4 - m 


•tnin 


f^\h» 


o 


UJ 


— -o- 


- 1 


^S 1 — •> 1 1 1 1 1 1 


>• M 


O MO •* 


ZOS^ 


X 


•>$•• #a 


> 1 


CM 1 ^X 1 1 | 1 1 1 



ro 

*muj mx it mx 

43- X •‘CVJIUX •x.a# 

— *< 43 -> 43 -' 3 X 43 ~ * 

hd^ | ^ 

iu«— uj<ouj<m ii t-u<x 
Hia-HzmHXHOi-Zcnr 
*-«X X*-»X mQJZLDmQCO^ 

ocaTHCtooaiooxaio 
3 TU.UJ •*iu.QZU,Uoi 3 U.*» 
** ** 



3 

O 

LU 



x 

X 



c^r^i-M—vO 
Q£ *• O •* 



I 



I 



I I 



I 



I I 
I I 



I 



I 



I 



so i 43 in i 

cn * j *m i 

LU>0^ •43—’ | 434)—’ I 'O— i 

I I I 

*^LU<OLU< I OLU< I LU< I 

'-HXQihS I hhJ I hX I 

hoc-ho: i mo; i moc i 

iLXOLLXn l oxo i xo t in i in t in I in t in 

mXlLmZU- I OSU. I 3 LL MMMI HlH 

# * ****#**###* 



• I 
X I 



I I 



• I 
X I 



I 

I 

I 

• I 
X I 



L)UUU 



't 


M 


insO 


r- 


00 


** 


o 


MCM 


co>t 


o 




oo 


o 


o 


p^- 


43 


4343 


43 43 


CO 


o 


COCO 


CO 


cO 












CO 






CJCJCJOO 











218 



1 1 


1 


| 


UJ 


(M 






LL 


LL 


CO 


CO 


»-4 


LU 


| 


t 


1 


1 


1 

1 


t t 


1 


t 


< 


X 


Q 


UJ 


O 


a 


G 


X 


H- 


> 


1 


i 


1 


1 


1 t 


1 






w 


X 


o 






CO 


a 


< 




1 


i 


1 


1 


1 


1 t 


1 


— 






G 


< 


UJ 


UJ 


X 




G 


G 


— 1 


i 


1 


y 


1 


t t 


1 


LU 


LU 


Q 


a 


a. 


> 


> 


UJ 


i— 


X 


< 


'N. 1 


1 


1 


i 


1 


t t 




w 


G 


X 


X 


CO 


al 


X 


> 


< 


iu 




s t 


l 


1 


i 


1 


1 1 


— 




<* 


G 


G 




G 


G 




G 


h— 


X 


^ 1 


l 


1 


i 


1 


1 


G 


h— 


a. 


G 




UJ 


O 


O 


CO 


X 




a 


N* | 


t 


1 


i 


1 


1 — 


u 


X 


CO 


X 


UJ 


UJ 








LU 


< 


— * 


^ 1 


l 


1 


i 


1 


1 G 


< 


<f 




G 


> 


(X 


X 


X 




►- 




h— 


S 1 


i 


1 


i 


1 


-J 


*** 


h- 


UJ 




a 


LL 


h* 


h- 


X 




a: 


< 


— 1 


l 


1 


i 


1 


— Q 




cO 


UJ 


LU 


CO 




o 


O 


a 


< 


LU 


G 




i 


1 


i 


1 


X 


CO 


X 


X 


> 


< 


2 


a 


O 


»-4 




H— 


X 


1 I 


i 


1 


i 


1 


G 


u 


G 


LL 


G 




a 


X 


X 


h- 


a£ 


H- 


UJ 


1 1 


t 


1 


i 


1 


w 00 


Q 


O 




CO 


h- 


u 


CO 


00 


<t 


UJ 


< 


h- 


t 1 


l 


1 


i 


1 


LU 






3 


< 


X 


UJ 






G 


H 


O 


h» 


t 


l 


1 


i 


1 


H* O 


UJ 


G 


G 




G 


TO 


< 


c 


X 


L 


G 


<f 


h* 1 




1 


t 


1 


X -X 


•* > 


*-4 


-u 


•4— 




•* 


— 


— 


UJ 




— 


— ^ 


-o 1 


G - 


► 1 — 


— i 


•» | 


C3\<\2. 


— X>*cU>*X'N*LU^*X^LLV.lL 


1— — 








a. t 


GX 


1 '“v 


^ i 1 


^ 1 — 






















i^NZsOs 


X 1 


-ON* 


1 N* 




N, 1 X 


UJ -G 


•4— 


- o 


— 


— ■ 4 


— 


— *— « 


— 


- 


s< - 




•^ct 


—♦-4 — 


M | 


N*U - 


1 — 


-G 


— | — 


X- -•« 


» G‘ 


■ LU 


S' 


■ LU. 


• <« 


- h— • UJ— UJ ■ 


o< 


■ — H • 


! 


N* - 


| — 


• G* 


■ | — 


Q 


—4 


-j 


X 


X 


< 


CL 


a 


*a: 


X 


X 


o 


< 1 


-LU 


1 


o 


1 


X 


Q 


LU 


M 




X 


G 


a 


a« 


■ UJ 


< 


a 


< 


^ 1 


- X 


1 


u 


1 


■-» Zt/) 


H- 1 


h- 


< 


UJ 


X 


-j 


u 




cO— i 


oc 


< 1 


•—4 


1 




1 


< QQSH- 


Q 


< 


x< 


}—» 


UJ 


CO 


CO 


1— 


XQC>- 


LL 


GO 1 


u 


1 




1 


agcoi'vtuK 




G 


LOXGH- 






< 


OUJO 


LL 


< 1 






1 


a£a£*-»>- 


LU 


h* 


LUQC<QCH- 


ooXlo 


OX 


-•»—< 


h-»X 


XX 1 


X 


X 


—4 


1 


ujuioauz 


g> 


2 


h-UJ 1 


LU^ 


UJ^LU 


COX 


H— LUOC 


ax 


aCh- l 


o 


—4 


-1 


1 


HhQZO 


Xh 


UJ 


ZKoCH 


2SZS S 


OXLL 


N> 


-G 


G— » 


< 


X 




LUX 1 ►-•CO -h- 


h- CO <t LU UJ LU X CO • CO • X LL CO 


< | LL COLL CO 


_J2t-t-oOcOOcO 


cazcazco 


LLX 


H-Q : 


ShSQQu.Qu.saa 


O^^hQQQ 


JLU3J-COQ 


Q>QOQhQ 


a x 


< 


XU 


< 


X 


h- 




UJ 


LUCO 


LL 


Q 




xoca<i— • 


<t h* — 4 


< 


H— 


G 


-LU 




UJ 


— • 




Q 


QQLU 


LL 




LU 


CtfH- 


2 X 


X — • G 


G 


UJ- X- X' 


• *HX 


x< 


• > 


- Z 


. o- 


■ 


>- 




» LU— >• 


-GaC 


G- 


X- >- G< 


■ — 


G -< 


-UJ 


-X 


LU 


-♦-4 


-CO 


-< 


-H- 


— H— 


QC - 


a 


-0C 


-QC - 


G UJ 


x -LU — » — *— » 


-< - 


2 UJ 


h-» 


*X 


h- 


UJ 


X 


a: 


< 


<t- G 




LU 


G 


GO 2 


- < 


h— h- X 


X 


< 1 


h- 


oo 


h* 


o 


< 


LL 


X 


X 


-O 


UJ 


X 


O 


x<to 


-oC 


X — < G 


X 


xc\iX(\l<c\JLU c\j<<\jLU(\ixc\JLLc\l^inh- 


in 


acc\j2nj in 


XXa. 


h-C\J<C\JG(\jX<\JUJC\J 


•H- 


• 


•<> 


•O 


•X 


•h- 


• M 


• 


• 


LU • 


UJ 


• 


• LU • 


-G 


• 


• X •< 


• h- • 


ujr-ah-xh- r-Gr~ i 


r*» i 


^NhNhvQXN 


Xr-.2fh»ir» 


G oc (\i qc r- oi r- uj r*. a; h- X f- 


JlLOlLOltOMU. LLULL JlL LLXU.X 


• 1— LL 


2U.OLLHU. 


<GUJ 


• UJ LL UJ LL GLL h- LL <t LL 



*-* i —*-4 — Q **<t •»< — O HU — UJ O' — — *-4 — — X hOh •*<-*> — — — — 

U* (/)• h» • X — 4 — LLJ «■ CX • Q£« LU» —4 • —hLL LL — LU «■ LL *» h* | h"HH • h* • CtZm CX • Q£ — 

UJ | I O I -H|h|3nih|(J|0*Ol UKIOI #.|MU-M|M|UJ|UJ|UJ| 
OtllKI Of- I < | h- t t- I < I — i | -H- I 02 IOI I OU IX»XtXl>j>t>| 
QC | G t cx I CNJH-I ISiOIOlSIUhllUJl X< | X I Of | 'OJ- oo. oo | 00 | »-i j H | m | 

lu i t u. I 12 |M| 0 | 0 |M|q,|u.io.i h- i uj a. <xx2 I z i Z i uj i uj i uj i 
uir iu,i -«f ihirfiQ^ihiuiilujioi i h-*a-»Ohvt<i i<ji<j luiuiui 
ZI3|N| GuJ I^Ihlhl^lOiai Jl 0—» I hXJOh -sf aC t OX t OX t UJ I LU I LU I 
M|(OIO| UJQ. | LU | </) | (/) | (il I O I O I (/) I <Q <*OGX a*h|hlhlQ^IXlQCI 



1 - 1 


- i <- i - 


1 - 


1 • 1 - 1 - 


« - i - 


»• 


<t- 


'-aa.s- l • 


1 • 1 • 1 - 


i - I 


t t 


i •» i 


| 


t t 1 


t t 


— X X 




GXS 1 


t t 1 


I t 


I l 




1 


1 ( | 


i t 


GO O 




UJ h-N, 1 


I ! 1 


i | 


i t 


i ! 


1 


i i i 


1 i 


X — 


u 


U— >s | 


i 1 i 


i ! 


l ) 


i -* *i 


1 


i i i 


i i 


■— 1 — — 


< 


X*H— X. | 


t I i 


i t 


t ! 


l inx 1 


1 






sOXU 


G 


< »C0\ | 


i I I 


i i 


t t 


I sOin 1 


1 


i t t 


— — 


uj'0m<t 


Or 


h-UJsOX> | 


i i i 


t t 


i i 


1 --4 1 


1 




-H C\J 


O »hG 


UJ 


GX | 


l t i 


1 l UJ 


i i 


1 'O— 1 


1 




X X 


Xvo — a 




—1 •sO'— | 


i t i 


1 1 


t i 


1 — H- | 


1 


i 1 1 


— » 1 


<wh»UJ 


UJ 


OCQ— H* 1 


t t i 


1 1 X 



I I I LU < | It—— h-LU<t Ot QLU<f I I I I I I 

• I - I - I -I -O -Q ••OX HX -Gh-XUJ -< •* I * I -I -I - I -I 

X I X I X I *+CC I X x XHXUJ XOXOXnmoCqCX XOwhqCX I X I X | X I X I X I 

in i in I in I xo i in— in— in:nin<f in< in<inoxo<rin:^tnGiL xrom i in t in i m i in i in i 

H I H M I 3 LL I r-H «H 2U. HQi-IQ^M2U,H |H|HMM|H| 

****** *************** ***** ************ 

m 'O r*» co 

sO 'O 'O'O 



OG 



219 



69 CONTINUE 
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PROGRAM PARAMETERS 
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TRANSMITTER IRANSMISSICN LINE LOSS 2*00 00 

TRANSMITTER ANTENNA GAIN 45.40 00 

RECEIVER SENSITIVITY - 40.00 00 M 

RECEIVER TRANSMISSION LINE LOSS 2.00 00 

RECEIVER ANTENNA GAIN 10.00 00 



FREQUENCY OF SIGNAL -F- INEG-HERT l I 5825.00000 
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TRANS-LOSS (OBI 01STANCE (K-HETERSI 



SIGNAL STRENGTH VERSUS DISTANCE FOR GIVEN SYSTEM PARAMETERS — 
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♦*♦*♦ FOR A RECEIVER INPUT LEVEL OF -40.00 O0M, DISTANCE 6ETREEN ANTENNAS CAN 0E AT LEAST 117.00 AM ♦♦♦♦♦ 



PRUGRAM PARAMETERS 
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TRANSMITTER POWER OUT 1000000*00 WATTS 

TRANSMITTER TRANSMISSION LINE LOSS 2.00 Ob 

TRANSMITTER ANTENNA GAIN *5.60 06 

RECEIVER SENSITIVITY 06 M 



FREQUENCY OF SIGNAL -F- 4HEG-H6RTZI 5825.00000 
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RECEIVER TRANSMISSION LINE LOSS 2.00 

RECEIVER ANTENNA GAIN 1 0-00 
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PLOT OP TRANS. -LOSS VERSUS DISTANCE 
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TRANS-LOSS 4DBI DISTANCE (K-NETERSI 



SIGNAL STRENGT H VERSUS DISTANCE EUR GIVEN SYSTEM PARAMETERS 
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•♦**• FOR A RECEIVER INPUT LEVEL OF -AO. 33 08N, OISTANCE BETWEEN ANTENNAS CAN 8E AT LEAST 138.00 KN 
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POLARIZATION IPULl. VERTICAL 

FREQUENCY OF SIGNAL IF I 40.00 HE G -HERTZ 

SURFACE CONDUCTIVITY I SI 0.00500 HTQ/METER 

SURFACE REFRACT! VITV I NS I 300.00 N-UN1TS 

SUM OF ELEVATION ANGLES I T E I - .00288 RADIANS 
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NOTE 10 FIND SIGNAL STRENGTH FROM GRAPH, SUBTRACT TRANS. -LOSS FROM 61.44 



© 

o 

3 






O 

3 



6 

o 



m 

A 

CD 



<6 

O 



o 

3 



6 



i 



< 



o« 
o ♦ 


• 












t3 


• • 


• 












♦ • 


o* 


• 












♦ O 




• 












♦ rn 


• 


• 












♦ 


♦ 


• 












* 


» 


• 












♦ 


• 


• 












♦ 


♦ 


• 












« 


♦ 


• 












♦ 


♦ 


• 












* 


♦ 


• 












♦ 


• 


• 












« 


« 


• 












« 


♦ 


• 












♦ 


f 


. 












• 


# 




• 










* 


• 




• 










♦ 


♦ 




« 










* 


o« 






• 








♦ o 


m ♦ 






• 








* in 








• 








• • 


<M # 






• 


















. 






• <NI 


♦ 








• 






# 


• 








• 






# 


• 








• 


• 




* 


• 










• 




♦ 


♦ 










• 




1 


« 










• 




♦ 


♦ 










• 




♦ 


♦ 










• 




♦ 


4 










• 




* 


♦ 










. 




• 


* 










• 




• 


♦ 












• 


♦ 


• 












• 


• 


♦ 












• 


« 


♦ 












• • 


* 


O ♦ 












• 


♦ o 


o ♦ 












• 


★ o 


• # 












• 


♦ * 














• 




-*♦ 












• 




• 












• 


• 


♦ 












• 


♦ 


♦ 












• • 


♦ 


♦ 












• 


* 


♦ 












• 


» 




o m 

O m 



f O 

'O © 



o 

m 



o 

3 



m 

m 

« 

CD 



>o 

<o 



O 

o 

o 



261 
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TRANSilfT ER ANTENNA GAIN 10.00 06 

RECEIVER SENSITIVITY 80.00 06 

RECEIVER TRANSMISSION LINE L3SS 2.00 06 

RECEIVER ANTENNA GAIN 10.00 08 



FREUUE NCV OF SIGNAL -F- JMEG-HERTZI 40.00000 
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PLOT UF TRANS. -LOSS VERSUS DISTANCE 
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TRANS-LOSS (OBI DISTANCE (K-HETERSI 



SIGNAL STRENGTH VERSUS DISTANCE FOR GIVEN SYSTEM PARAMETERS — ■ 
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PROGRAM PARAMETERS 
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TRANSMITTER ANTENNA GAIN 10.00 00 

RECEIVER SENSITIVITY ** 00.00 00 H 

RECEIVER TRANSMISSION LINE LOSS 2.00 00 



FREQUENCY OF SIGNAL -F- (NEG-HERTZI 40.00000 
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TRANS-LOSS 4 06 I 01STANCE IK-NETERSI 



SIGNAL STRENGTH VERSUS DISTANCE FOR GIVEN SYSTEM PARAMETERS — 
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APPENDIX C 



COMPUTER SIMULATION 2 LISTING AND OUTPUT 

This Appendix contains the program and output for 
Simulation 2. The output consists of Tables which are 
referenced by number in the main text of the thesis. At 
the top of each data Table, is an explanation of the 
parameters used in the run. Each Table represents one 
run of the program. 
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CALL MESSAG<»DNDX=-40.0 GRADIEIMT=15 KM$ • , 100, 0 . 5» 3 . 51 



CALL CURVEl XRA Y A, YRAYA .COUNTA » -1 ) 
CALL ENDPL(-l) 
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CALL CURVE ( XRA Y B» Y RAY B t COUNTB 
CALL ENCPL ( -2 ) 



FOLLOWING SUBROUTINE IS USED TO PLOT THE FIXES FOR EMITTER 
SUBROUTINE PICHC ( XRAYC, YRAYC , COUN TC.MAXERR ) 
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CALL CURVE! >RAYC»YRAYC »COUNTC • -1) 
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SUBROUTINE GRAFIK PQSRAY, AVGRA Y) 
REAL* 8 F0SRAY(4I»AVGRAYI4I 
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CALL CURVE! FOSKAY, AVURAY,4,1) 
CALL ENCPL ( -5 I 



TABLE C-l 

Simula t ion 2, Run 1. Simulation to check program. 
Gradient area equal to 15 Km. All values are in 
kilometers. Emitter D has system error of 2.0 
degrees . 



DNDX= - 1.300 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100 • 00101 2 AVERAGE X= 80.000079 
AVERAGE MISSED DIST ANCE= 0.001015 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.000441 AVERAGE X=119. 999938 
AVERAGE MISSED DIST ANCE= 0.000445 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 002003 AVERAGE X=149. 999896 
AVERAGE MISSED DIST ANCE= 0.002006 



EMITTER D AVERAGE ERROR 

AVERAGE Y= 99.999820 AVERAGE X= 79.999996 
AVERAGE MISSED DIST ANC E= 0. 000180 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100. 000853 AVERAGE X= 80.000006 
AVERAGE MISSED DIST ANCE= 0.000853 



EMITTER 
AVERAGE 
A V ERAGE 



B AVERAGE ERROR 

Y= 70.000413 AVERAGE X=119. 999964 
MISSED DIST ANC E= 0. 000415 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 002263 AVERAG E X = 149 .999909 
AVERAGE MISSED DIST ANC E= 0. C02264 



EM ITTER 
AVERAGE 
AVERAGE 



D AVERAGE ERROR 

Y= 99.999889 AVERAGE X= 79.999996 
MISSED DI ST ANC E= 0. 000111 
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TABLE C-l 



Simulation 2, Run 1 (.continued). 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100. 000803 AVERAGE X = 80.000032 
AVERAGE MISSED DIST ANCE= 0. C00804 



EMITTER 8 AVERAGE ERROR 

AVERAGE Y= 70.000464 AVERAGE X=119. 999922 
AVERAGE MISSED DISTANCE* 0.000471 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 002500 AVERAGE X=149. 999883 
AVERAGE MISSED DISTANCE* 0.002503 



EMITTER D AVERAGE ERROR 

AVERAGE Y= 99.999870 AVERAGE X= 79.999997 
AVERAGE MISSED DISTANCE* 0.000130 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100. 000808 AVERAGE X= 80.000025 
AVERAGE MISSED DISTANCE* 0.000809 



EMITTER 8 AVERAGE ERROR 

AVERAGE Y= 70.000431 AVERAGE X=119. 999934 
AVERAGE MISSED DISTANCE* 0.000436 



EMITTER 
AV ERaGE 
AVERAGE 



C AVERAGE ERROR 

Y=250. 002532 AVERAGE X=149. 999904 
MISSED DISTANCE* 0.002533 



EMITTER D AVERAGE ERROR 

AVERAGE Y= 99.999909 AVERAGE X= 79.999997 
AVERAGE MISSED DISTANCE* 0.000091 
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TABLE C-2 



Simulation 2, Run 2. Gradient area equal to 15 Kn. 
All values are in kilometers. Emitter D has system 
error of 2.0 degrees. 



DN DX= - 


40.000 


EMITTER 
AV ERAGE 
AV ERAGE 


A AVERAGE ERROR 

Y=100. 043380 AVERAGE X= 80.002240 
MISSEO OIST ANC E = 0.043438 


EMITTER 
AVERAGE 
AV ERAGE 


B AVERAGE ERROR 

Y= 70.019390 AVERAGE X=119. 998662 
MISSED DIST ANC E= 0.C19436 


EMITTER 

AVERAGE 

AVERAGE 


C AVERAGE ERROR 

Y=250. 127474 AVERAGE X=149. 999134 
MISSED DIST ANC E= 0. 127477 


EMITTER 

AVERAGE 

AVERAGE 


D AVERAGE ERROR 

Y=102. 759927 AVERAG E X = 81 . 703879 
MISSED OIST ANC E= 3.243517 


EMITTER 

AVERAGE 

AVERAGE 


A AVERAGE ERROR 

Y=100. 040231 AVERAGE X= 80.001522 
MISSEO DIST ANC E= 0.040260 


EM ITTER 
AVERAGE 
AVERAGE 


B AVERAGE ERROR 

Y - 70.020793 AVERAGE X = 119 .997969 
MISSED DIST ANC E= 0.020892 


EM ITTER 
AVERAGE 
AVERAGE 


C AVERAGE ERROR 

Y=250. 126468 AVERAGE X= 149. 999207 
MISSED DISTANCE= 0.126471 


EMITTER 
AV ERAGE 
AVERAGE 


D AVERAGE ERROR 

Y= 98.649566 AVERAGE X= 79.368710 
MISSED QI ST ANC E = 1.490704 
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TABLE C-2 

Simulation 2, Run 2 (continued) 



EM ITTER 
AV ER AGE 
AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 
AV ERAGE 
AVERAGE 



EMITTER 
AV ERAGE 
AV ERAGE 



EM ITTER 
AV ERAGE 
AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 
AVERAGE 
AV ERAGE 



EMITTER 
AVERAGE 
AV ERAGE 



A AVERAGE ERROR 

Y=100 *042220 AVERAGE X* 80.002135 
MISSED 01 ST ANC E = 0.042274 



8 AVERAGE ERROR 

Y= 70.020391 AVERAGE X=119. 997753 
MISSED DISTANCE* 0.020514 



C AVERAGE ERROR 

Y=250. 130829 AVERAGE X=149. 998735 
MISSED DISTANCE* 0.130835 



D AVERAGE ERROR 

Y=102. 725045 AVERAGE X* 80.692302 
MISSED DISTANCE* 2.811610 



A AVERAGE ERROR 

Y=100. 042224 AVERAGE X* 80.002115 
MISSED DISTANCE* 0.042277 



B AVERAGE ERROR 

Y= 70.019693 AVERAGE X=119. 998251 
MISSED DISTANCE* 0.019770 



C AVERAGE ERROR 

Y=250. 127309 AVERAGE X=149. 999286 
MISSED DISTANCE* 0.127311 



D AVERAGE ERROR 

Y=102. 410316 AVERAGE X* 79.863946 
MISSED DISTANCE* 2.414153 
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I 









TABLE C-3 

Simulation 2, Run 3. Gradient area equal to 15 Km. 
All values are in kilometers. Emitter D has system 
error of 2.0 degrees. 



DNDX = - 80*000 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100. 086917 AVERAGE X* 80.004486 
AVERAGE MISSE0 DISTANCE* 0.C87032 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.038854 AVERAGE X=119. 997338 
AVERAGE MISSE0 DISTANCE* 0.038945 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 256023 AVERAGE X=149. 998342 
AVERAGE MISSED DISTANCE* 0.256029 



EMITTER D AVERAGE ERROR 

AVERAGE Y=102 • 7 5992 7 AVERAGE X* 81.703879 
AVERAGE MISSED DISTANCE* 3.243517 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100. 080602 AVERAGE X* 80.003047 
AVERAGE MISSED DISTANCE* 0.080659 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.041659 AVERAGE X=119. 995953 
AVERAGE MISSED DISTANCE* 0.041855 



EMITTER 

AVERAGE 

AVERAGE 



C AVERAGE ERROR 

Y=250. 253971 AVERAGE X=149. 998479 
MISSED DISTANCE* 0.253975 



EMITTER 

AVERAGE 

AVERAGE 



D AVERAGE ERROR 

Y= 98.649566 AVERAGE X* 79.368710 
MISSED DISTANCE* 1.490704 
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Simulation 2, 


TABLE C-3 

Run 3 (continued) . 



EMITTER 
AVERAGE 
AV ERAGE 


A AVERAGE ERROR 

Y=I00. 084615 AVERAGE X= 80.004267 
MISSED DIST ANCE= 0.C84722 


EMITTER 


B AVERAGE ERROR 



AVERAGE Y= 70.040865 AVERAGE X=119. 995525 
AVERAGE MISSED DIST ANC E= 0.041X09 

EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 262486 AVERAGE X=149 .997556 
AVERAGE MISSED DIST ANCE= 0.262498 



EMITTER 

AVERAGE 

AVERAGE 


D AVERAGE ERROR 

Y=102. 725045 AVERAGE X= 80.692302 
MISSED DI ST ANC E= 2.811610 


EMITTER 

AVERAGE 

AVERAGE 


A AVERAGE ERROR 

Y=100. 084590 AVERAGE X= 80.004232 
MISSED DIST ANC E= 0. 084696 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.039468 AVERAGE X=119. 996522 
MISSED DI ST ANC E= 0. C39621 


EMITTER 
AV ERAGE 
AVERAGE 


C AVERAGE ERROR 

Y=250. 255841 AV ERAG E X = 149 .99 8652 
MISSED DI ST ANCE = 0. 255845 


EM ITTER 
AVERAGE 
AVERAGE 


D AVERAGE ERROR 

Y=102. 410316 AVERAGE X= 79.863946 
MISSED DI ST ANC E= 2.414153 
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TABLE C-4 



Simulation 2, Run 4. Gradient area equal to 15 Km. 
All values are in kilometers. Emitter D has system 
error of 2.0 degrees. 

r 

ON D X= -120.000 



EMITTER 

AVERAGE 

AVERAGE 


A AVERAGE ERROR 

Y=100. 130485 AVERAGE X* 80.006726 
MISSED 01 ST ANC E = 0* 130658 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.058320 AVERAGE X=119. 996020 
MISSED DI ST ANC E= 0.058456 


EMITTER 
AV ERAGE 
AVERAGE 


C AVERAGE ERROR 

Y=250.38471 9 AVERAGE X* 149. 99 7550 
MISSED DISTANCE= 0.384727 


EMITTER 


D AVERAGE ERROR 



AVERAGE Y=102. 759927 AVERAGE X* 81.703879 



AVERAGE 


MISSED DISTANCE* 3.243517 


EM I TTER 
AVERAGE 
AV ERAGE 


A AVERAGE ERROR 

Y=100. 121006 AVERAGE X* 80.004569 
MISSED DISTANCE* 0.121092 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.062523 AVERAGE X=119. 993945 
MISSED DISTANCE* 0.062815 


EMITTER 
AVERAGE 
AV ERAGE 


C AVERAGE ERROR 

Y=250. 381620 AVERAGE X=149. 997752 
MISSED DISTANCE* 0.281627 


EMITTER 
AVERAGE 
AV ERAGE 


0 AVERAGE ERROR 

Y= 98.649566 AVERAGE X* 79.368710 
MISSED DISTANCE* 1.490704 
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TABLE C-4 

Simulation 2, Run 4 (continued). 



EM ITTER 
AVERAGE 
AV ERAGE 


A AVERAGE ERROR 

Y=100. 127041 AVERAGE X= 80.006394 
MISSED DISTANCE* 0.127202 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.061337 AVERAGE X = 119 . 993306 
MISSEO DISTANCE* 0.061702 


EMITTER 

AVERAGE 

AVERAGE 


C AVERAGE ERROR 

Y=250. 394293 AVERAGE X=149. 996377 
MISSED DISTANCE* 0.394310 


EMITTER 
AV ERAGE 
AVERAGE 


D AVERAGE ERROR 

Y=1D2. 725045 AVERAGE X* 80.692302 
MISSED DISTANCE* 2.811610 


EMITTER 

AVERAGE 

AVERAGE 


A AVERAGE ERROR 

Y=100. 126989 AVERAGE X= 80.006345 
MISSED DISTANCE* 0.127147 


EM I TTER 
AVERAGE 
AVERAGE 


B AVERAGE ERROR 

Y= 70.059245 AVERAGE X=119. 994801 
MISSED DISTANCE* 0. 059473 


EMI TTER 
AVERAGE 
AVERAGE 


C AVERAGE ERROR 

Y=250. 384520 AVERAGE X* 149 . 998019 
MISSED DISTANCE* 0.384525 


EMITTER 

AVERAGE 

AVERAGE 


D AVERAGE ERROR 

Y=102. 410316 AVERAGE X* 79.863946 
MISSED DISTANCE* 2.414153 
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TABLE C-5 

Simulation 2, Bun 5. Gradient area equal to 15 Km. 
All values are in kilometers . Emitter D has system 
error of 2.0 degrees. 

ON DX= -150.000 



EMITTER 

AVERAGE 

AVERAGE 


A AVERAGE ERROR 

Y=100. 163182 AVERAGE X* 80.008403 
MISSED DISTANCE* 0.163398 


EMITTER 


8 AVERAGE ERROR 



AVERAGE Y= 70.072923 AVERAGE X=119. 995035 



AV ERAGE 


MISSED DISTANCE* 0.073092 


EMITTER 
AVERAGE 
AV ERAGE 


C AVERAGE ERROR 

Y=250. 481339 AVERAGE X=149. 996957 
MISSED DISTANCE* 0.481348 


EM ITTER 
AVERAGE 
AVERAGE 


D AVERAGE ERROR 

Y=102. 759927 AVERAGE X* 81.703879 
MISSED DISTANCE* 3.243517 


EMITTER 

AVERAGE 


A AVERAGE ERROR 

Y*100. 151332 AVERAGE X* 80.005709 



AVERAGE MISSED DISTANCE* 0.151440 



EM ITTER 
AVERAGE 
AVERAGE 


8 AVERAGE ERROR 

Y= 70.078171 AVERAGE X = 119 .992447 
MISSED DISTANCE* 0.078535 


EMITTER 
AVERAGE 
AV ERAGE 


C AVERAGE ERROR 

Y=250. 477453 AV ERAG E X* 149 . 99 7206 
MISSED DISTANCE* 0.47746 2 


EMITTER 

AVERAGE 

AVERAGE 


D AVERAGE ERROR 

Y= 98.649566 AVERAGE X* 79.368710 
MISSED DISTANCE* 1.490704 
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Simulation 2, 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 
A V ERAGE 
AV ERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 
AVERAGE 
AV ERAGE 



EMITTER 
AVERAGE 
AV ERAGE 



EMITTER 
AVERAGE 
AV ERAGE 



TABLE C-5 
Run 5 (continued). 



A AVERAGE ERROR 

Y=130. 158883 AVERAGE X= 83.007987 
MISSED DISTANCE 3 0.159084 



0 AVERAGE ERROR 

Y= 70.076691 AVERAGE X=119. 991647 
MISSED DISTANCE 3 0.077144 



C AVERAGE ERROR 

Y=250. 493245 AVERAGE X=149. 995494 
MISSED DISTANCE 3 0.493266 



D AVERAGE ERROR 

Y=102. 725045 AVERAGE X= 80.692302 
MISSED DISTANCE 3 2.811610 



A AVERAGE ERROR 

Y=100. 158809 AVERAGE X= 80.007927 
MISSED DISTANCE 3 0.159007 



B AVERAGE ERROR 

Y= 70.074079 AVERAGE X=119. 993514 
MISSED DISTANCE 3 0.074363 



C AVERAGE ERROR 

Y=250. 481126 AVERAGE X=149. 997545 
MISSED DISTANCE 3 0.481133 



D AVERAGE ERROR 

Y=102. 410316 AVERAGE X= 79.863946 
MISSED DISTANCE 3 2. 414153 
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TABLE C-6 

Simulation 2, Run 6. Gradient area equal to 30 Km. 
All values are in kilometers. Emitter D has system 
error of 2.0 degrees. 



ON D X* - 


80.000 


EM ITTER 
AVERAGE 
AVERAGE 


A AVERAGE ERROR 

Y=100. 146546 AVERAG E X= 80 .007539 
MISSED OI ST ANCE* 0.146740 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.044350 AVERAGE X = 119 .996988 
MISSED DISTANCE= 0.044452 


EMITTER 

AVERAGE 

AVERAGE 


C AVERAGE ERROR 

Y=250. 494021 AVERAGE X=149. 996880 
MISSED DISTANCE* 0.494030 


EMITTER 

AVERAGE 

AVERAGE 


D AVERAGE ERROR 

Y=102. 759927 AVERAGE X= 81.703879 
MISSED DISTANCE* 3.243517 


EMITTER 
AV ERAGE 
AVERAGE 


A AVERAGE ERROR 

Y=100. 135913 AVERAGE X= 80.005123 
MISSED DISTANCE* 0.136009 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.047533 AVERAGE X=119. 995420 
MISSED DISTANCE* 0.047753 


EM ITTER 
AV ERAGE 
AV ERAGE 


C AVERAGE ERROR 

Y=250. 490033 AVERAGE X=149. 997136 
MISSED DISTANCE* 0.490042 


EMITTER 
AVERAGE 
AV ERAGE 


D AVERAGE ERROR 

Y= 98.649566 AVERAGE X* 79.368710 
MISSED DISTANCE* 1.490704 
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Simulation 2, 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EM ITTER 
AVERAGE 
AVERAGE 



EM ITTER 
AVERAGE 
AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 

AVERAGE 

AVERAGE 



TABLE C-6 

Run 6 (continued) . 



A AVERAGE ERROR 

Y=130. 142684 AVERAGE X = 80,007169 
MISSEO DIST ANCE = 0. 142864 



B AVERAGE ERROR 

Y= 70.046624 AVERAGE X=119. 994932 
MISSED DISTANCE= 0.046898 



C AVERAGE ERROR 

Y=250. 506230 AVERAGE X=149. 995379 
MISSED DISTANCES 0.506251 



D AVERAGE ERROR 

Y=102. 725045 AVERAGE X= 80.692302 
MISSED DISTANCES 2.811610 



A AVERAGE ERROR 

Y=100. 142623 AVERAGE X= 80.007113 
MISSED DISTANCES 0.142801 



B AVERAGE ERROR 

Y= 70.045045 AVERAGE X=119. 996062 
MISSED DISTANCES 0.045217 



C AVERAGE ERROR 

Y=250. 493808 AVERAGE X= 149 .997484 
MISSED DISTANCES 0.493814 



D AVERAGE ERROR 

Y=102. 410316 AVERAGE X= 79.863946 
MISSED DISTANCE* 2.414153 
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TABLE C-7 



Simulation 2, Run 7. Gradient area equal to 30 Km. 
All values are in kilometers. Emitter D has system 
error of 2.0 degrees. 



DNDX= -120.000 



EMITTER 
AV ERAGE 
AV ERAGE 



A AVERAGE ERROR 

Y=100. 220026 AVERAGE X= 80.011285 
MISSE0 0IST ANCE* 0.220315 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.066523 AVERAGE X=119. 995517 
AVERAGE MISSED DIST ANCc= 0.066679 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 742346 AVERAGE X=149. 995360 
AVERAGE MISSED DISTANCE* 0.742361 



EMITTER D AVERAGE ERROR 

AVERAGE Y=102. 759927 AVERAGE X* 81.703879 
AVERAGE MISSED DISTANCE* 3.243517 



EMITTER A AVERAGE ERROR 

AVERAGE Y=100. 204084 AVERAGE X* 80.007670 
AVERAGE MISSED 0ISTANCE= 0. 204228 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.071281 AVERAGE X = 119 .993181 
AVERAGE MISSED DISTANCE* 0. €71606 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 736341 AVERAGE X = 149. 995739 
AVERAGE MISSED DISTANCE* 0.336353 



EMITTER 

AVERAGE 

AVERAGE 



D AVERAGE ERROR 

Y= 98.649566 AVERAGE X* 79.368710 
MISSED DISTANCE* 1.490704 
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TABLE C-7 

Simulation 2, Run 7 (continued) 



EMITTER 

AVERAGE 

AVERAGE 



EMITTER 
AV ERAGE 
AVERAGE 



EM ITTER 
AV ERAGE 
AVERAGE 



EMITTER 
AV ERAGE 
AV ERAGE 



EMITTER 
AV ERAGE 
AVERAGE 



EMITTER 
AV ERAGE 
AV ERAGE 



EMITTER 
AVERAGE 
AV ERAGE 



EMITTER 

AVERAGE 

AVERAGE 



A AVERAGE ERROR 

Y=100. 214249 AVERAGE X= 80.010728 
MISSED DI ST ANC E= 0.214518 



B AVERAGE ERROR 

Y= 70.069917 AVERAGE X=119. 992449 
MISSED DI ST ANC E= 0.070324 



C AVERAGE ERROR 

Y=250. 760544 AVERAGE X=149. 993117 
MISSED DI ST ANC E= 0.760575 



D AVERAGE ERROR 

Y=102. 725045 AVERAGE X= 80.692302 
MISSED DIST ANCE= 2.811610 



A AVERAGE ERROR 

Y=100. 214143 AVERAGE X= 80.010646 
MISSED DI ST ANCE= 0.214408 



B AVERAGE ERROR 

Y= 70.067568 AVERAGE X=119. 994138 
MISSED OIST ANCE= 0.067822 



C AVERAGE ERROR 

Y=250. 742099 AVERAGE X=149. 996269 
MISSED DISTANCE^ 0.742109 



D AVERAGE ERROR 

Y=102. 410316 AVERAGE X= 79.863946 
MISSED DIST ANCE= 2.414153 
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TABLE C-8 

Simulation 2, Run 8. Gradient area of 30 Km. System 
induced error of 1.5 degrees for all LOB's. Emitters 
A, B and C have refracted LOB's correlated with 
straight LOB's. All values are in kilometers. 



DN DX= -100.000 



EMITTER 

AVERAGE 

AVERAGE 



A AVERAGE ERROR 

Y 3 96.136822 AVERAGE X= 80.311562 
MISSED DISTANCE 3 3.875721 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.055440 AVERAGE X=119. 996250 
AVERAGE MISSED DISTANCE 3 0.055567 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 618112 AVERAGE X=149. 996120 
AVERAGE MISSED DISTANCE 3 0.618124 



EMITTER D AVERAGE ERROR 

AVERAGE Y=101.991913 AVERAGE X 3 81.272969 
AVERAGE MISSED DISTANCE 3 2.363930 



EMITTER A AVERAGE ERROR 

AVERAGE Y=1D1. 483913 AVERAGE X= 80.127291 
AVERAGE MISSED DISTANCE 3 1.489362 



EMITTER B AVERAGE ERROR 

AVERAGE Y= 70.059409 AVERAGE X=119. 994296 
AVERAGE MISSED DISTANCE 3 0.C59683 



EMITTER C AVERAGE ERROR 

AVERAGE Y=250. 613116 AVERAGE X=149. 996437 
AVERAGE MISSED DISTANCE 3 0.613126 



EMITTER D AVERAGE ERROR 

AVERAGE Y=100. 731179 AVERAGE X 3 79.514574 
AVERAGE MISSED DISTANCE 3 0.877645 
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TABLE C-8 

Simulation 2, Run 8 (continued). 



EMITTER 
AVERAGE 
AV ERAGE 


A AVERAGE ERROR 

Y= 98,926232 AVERAGE X* 80,760453 
MISSEO DISTANCE* 1.315779 


EMITTER 

AVERAGE 

AVERAGE 


B AVERAGE ERROR 

Y= 70.058276 A VERAG E X* 119 .993687 
MISSED DISTANCE* 0.058617 


EMITTER 

AVERAGE 

AVERAGE 


C AVERAGE ERROR 

Y=250. 633317 AVERAGE X=149. 994248 
MISSED DISTANCE* 0.633343 


EMITTER 
AV ERAGE 
AVERAGE 


D AVERAGE ERROR 

Y=101. 946699 AVERAGE X* 80.515853 
MISSED DISTANCE* 2.013887 


EMITTER 

AVERAGE 

AVERAGE 


A AVERAGE ERROR 

Y=101. 203482 AVERAGE X* 79.673178 
MISSED DISTANCE* 1. 247070 


EMITTER 
AV ERAGE 
AVERAGE 


8 AVERAGE ERROR 

Y= 70.056309 AVERAGE X=119 .995098 
MISSED DISTANCE* 0.056522 


EMITTER 

AVERAGE 

AVERAGE 


C AVERAGE ERROR 

Y=250. 617883 AVERAGE X=149. 996876 
MISSED DISTANCE* 0.617891 


EMITTER 

AVERAGE 

AVERAGE 


D AVERAGE ERROR 

Y*10 1.64986 8 AVERAGE X* 79.896641 
MISSED DISTANCE* 1.653102 



310 



I 









; 






LIST OF REFERENCES 



1 . National Oceanic and Atmospheric Administration Document 
NOAA-S/T 76-1562, U.S. Standard Atmosphere, 1976, t. 20, 
October 1976. 

2. Lake, J. S., "Atmospheric Refraction: A Neglected 

Fundamental Comes of Age", Defense Electronics, v. 14 
no. 10, pp. 22-28, October 1982 . 

3. Beach, J. B., "Part II, Atmospheric Effects on Radio 
Wave Propagation", Defense Electronics , v. 12 no. 10, 
pp. 75-83, October 1980. 

4. Naval Research Laboratory Report 7725, Geophysical 
Aspects of Atmospheric Refraction , by C. G. Purves, 
p. 14, 7 June 1974. 

5. Davidson, K. L., Atmospheric Factors in Electromagnetics/ 
Electro-Optics , material for classroom presentation. 

Naval Postgraduate School, Monterey, California, 

pp. 3-1 - 7-31, 1982. 

6. Meeks, M. L. , Radar Propagation at Low Altitudes , 
pp. 9-26, Artech House, Inc., 1982. 

7. Skolnik, M. I., Introduction to Radar Systems , pp. 448- 
454, McGraw-Hill, Inc., 1980. 

8. Transportation Systems Center - National Aeronautics 
and Space Administration Report 71-6, Atmospheric 
Transmission Handbook , by W. I. Thompson III, p. 3, 
February 1971. 

9. Office of Telecommunications Report 76-105, Re fr activity 
and Rainfall Data for Radio Systems Engineering , by 

C. A. Samson, pp. 6-9, September 1976. 

10. Harrington, J. B. and Callaghan, T. G. , "UHF/VHF 
Direction Finding", Watkins - Johnson Company Technotes , 
v. 8 no. 1, pp. 1-7, February 1981. 

11. McCartney, E. J., Optics of the Atmosphere , pp. 100-112, 
Wiley S Sons, Inc."^ 197 6 . 

12. Patton, T. N., Petrovic , J. J., and Teti, J., "Propagation 
Effects of Millimeter Wave Fire Control Systems", S. J. 
Johnston, editor. Millimeter Wave Radar, p. 170, 1980. 



311 



13. 



Kaelin, C. M. , EA-6B Mi ssion Planning: Tactical 

Considerations for ESM/ECM Employment Under Non-standard 
Refractive Conditions , Masters Thesis. Naval 
Postgraduate School, Monterey, California, p. 124, 1980. 

14. Naval Ocean Systems Center Technical Note 810, PROCAL- 

PREOS: A Fast-Running Airborne FLIR System Performance 

Model, by F. P. Snyder, p. 1, February 1980. 

15. Fraser, A. B. and Mach, W. H., "Mirage”, Scientific 
American, pp. 32-36, 1980. 

16. Zuev, V. E., Propagation of Visible and Infrared 
Radiation in the Atmosphere , op. 348-349, Wilev 5 Sons, 
Inc., 1974. 

17. Defense Advanced Research Projects Agency Report 1171, 
Degradation of Laser Systems by Atmospheric Turbulence , 
by R. F. Lutomirski, R. E. Huschke, W. C. Meecham, 

and H. T. Yura, pp. 1-230, June 1973. 

18. Cooper, A. W. , Electro-Optic Principles and Devices , 
material for classroom presentation. Naval Postgraduate 
School, Monterey, California, p. A-40, 1982. 

19. Hall, F. F., "Index of Refraction Structure Parameter 

in the Real Atmosphere - an Overview", Technical Digest, 
Optical Propagation through Turbulence, Rain, and Fog , 
p. TuCl-1, August 1977. 

20. Naval Postgraduate School Report 61-78-003, Optical 
Resolution in the Turbulent Atmosphere of the Marine 
Boundary Layer , by E. C. Crittenden, A. W. Cooper, E. A. 
Milne, 6. W. Rodeback, R. L. Armstead, S. H. Kalmback, 

D. Land, and B. Katz, pp. 2-16, February 1978. 

21. Crittenden, E. C., Cooper, A. W. , Milne, E. A., 

Rodeback, G. W. , and Kalmback, S. H., "Multiwavelength 
Extinction and Index Fluctuation Measurements" , Advisory 
Group for Aerospace Research and Development Conference 
Proceedings - 300, Special Topics in Optical Propagation, 
pp. 13-1 - 13-2, July 1981. 

22. McMillan, R. W. and Bohlander, R. A., "Millimeter Wave 

Atmospheric Turbulence Measurements: Preliminary 

Results and Instrumentation for Future Measurements" , 
Optical Engineering , v. 22 no. 1, pp. 33-36, January/ 
February 1983. 



312 



IV> 



23 . 



USA Ballistic Research Laboratories Memorandum Report 
No. 2243, Diurnal Cycles of Refractive Index Structure 
Function Coefficient , by M. L. Wesely and E.C. Alcaraz, 
p. 9, November 1972. 

4. U.S. Army Field Manual FM 17-12, Tank Gunnery, pp. 4-9 - 
4-10, 21 March 1977. 

25. Environmental Science Services Administration (ESSA) 
Research Laboratories Technical Report 79-ITS 67, 
Prediction of Tropospheric Radio Transmission Loss over 
Irregular Terrain - A Computer Method-1968 , by A. G. 
Longley, and P. D. Rice, p. 1, July 1968. 

26. Callaghan, J. M. , Tropoplot: An Improved Fortran 

Computer Program for Prediction of Long-Term Median 
Tropospheric Radio Transmission Loss Over Inrregular 
Terrain, Masters Thesis, Naval Postgraduate School, 
Monterey, California, pp. 10-42, 1973. 

27. U.S. Department of Commerce - National Bureau of 
Standards Monograph 92, Radio Meteorology , by 
B.R. Bean and E.J. Dutton, p. 63, March 1966. 

28. U.S. Army Field Manual FM 24-24, Radio S Radar 
Reference Data , pp. 1-12, May 1977. 

29. U.S. Army Field Manual FM 24-25, Wire and Multichannel 
Re ference Data , pp. 2-16, June 1977 . 

30. Naval Ocean Systems Center Technical Document 481, 

IREPS Revision 2.0 User's Manual , by H. V. Hitney , 

A"! Paulus , C~. P^ Hattan , W. D. Anderson, and G. E. 
Lindem, p. 30, September 1981. 



313 



BIBLIOGRAPHY 



GTE Sylvania Contract: MDA-904-78-C-0511, Radiosonde Data 

Analysis III : World Contour Maps , G. K. Miller, W. B. 

Moreland, and L. N. Ortenburger, December 1979. 

GTE Sylvania Contract: MDA-904-87-C-0511 , Radiosonde 

Data Analysis III: Summary Maps of Observed Data , L. N. 

Ortenburger , S"I B~^ Lawson , and G^ IT! Miller , December 1978 . 

National Bureau of Standards Monograph 22, Climatic 
Charts and’ Data of the Radio Refractive Index for the United 
States and the World , by B. R. Bean, J. D. Horn, and A. M. 
Ozanich Jr. , November 1960. 

Van Brunt, L. B., Applied ECM , EW Engineering, Inc., 1978. 

U.S. Army Atmospheric Sciences Laboratory, Effects of 
Atmospheric Turbulence on FLIR Performance , by H. A. Quevado 
and D. E. Snider, January 1983. 

B. Khular, K. Thyagarag jan , and A. K. Ghatak, "A note on 
Mirage Formation”, American Journal of Physics , v. 45 no. 1, 
January 1977. 

U.S. Army Field Manual FM 100-5, Operations , 1 July 1976. 

U.S. Army Field Manual FM 34-10, Military Intelligence 
Battalion (Combat Electronic Warfare Intelligence) (Division) , 
3 July _ 19 81 . 

Environmental Science Services Administration (ESSA) Research 
Laboratory Technical Report 65-ITS 58-2, Tabulations of 
Propagation Data over Irregular Terrain in the 230-to-9200 
MHz Frequency Range - Part II Fritz Peak Receiver Site , by 
P. L. McQuate, J. M. Harman, M. E. Johnson, and A. P. Batsis, 
December 1968. 

Environmental Science Services Administration (ESSA) Research 
Laboratory Technical Report 65-ITS 58-3, Tabulations of 
Propagation Dave over Irregular Terrain in the 230-to-9200 
MHz Frequency Range - Part III North Table Mountain - Golden , 
by P. L. McQuate, J. M. Harman, M. E. McClanahan, A. P. 
Batsis, July 1970. 



314 



Richter, J. H., Hughes, H. G. , and Rose, R. B., "Electro- 
magnetic Propagation Assessment", Electromagnetic 
Propagation Division, Naval Ocean Systems Center Technical 
Document 260, Proceedings of Conference on Atmospheric 
Refractive Effects Assessment, June 1979. 



315 



INITIAL DISTRIBUTION LIST 



No. Copies 
2 



1. Defense Technical Information Center 
Cameron Station 

Alexandria, Virginia 22314 

2 . Deputy Under Secretary of Army of Operations 
Research 

Room 2E261, Pentagon 
Washington, D.C. 20310 

3. Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93943 

4. Chairman, Electronic Warfare Academic 

Group, Code 73 
Naval Postgraduate School 
Monterey, California 93943 

5. Professor Gordon E. Schacher, Code 61Sq 
Department of Physics and Chemistry 
Naval Postgraduate School 

Monterey, California 93943 

6. Professor Kenneth L. Davidson, Code 63Ds 
Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93943 

7 . Commander/Director 

U.S. Army Atmospheric Sciences Laboratory 
Attn: DEDAS-AE-E (Dr. Don Snider) 

White Sands Missile Range, New Mexico 88002 

8 . Commander 

U.S. Army Electronic Proving Grounds 

Attn: STEEP-MT-T 

Ft. Huachuca, Arizona 85613 

9. Lieutenant Colonel Thomas Houser, USA 
Dep. C-of-S for Ops £ Plan(DCSOPS-FDI ) 
Washington, D.C. 20310 

10. Captain Theodore Mouras , USA 
Post Office Box 320 
Ft. Huachuca, Arizona 85613 



1 

2 

1 

1 

1 

3 

2 

1 

1 



316 



Thesis 

M8635 

c.l 



onoi: 70 

202572 



Mouras 

Effects of atmo- 
spheric refraction on 
U. S. ground warfare. 



5 dec ee 0 0 0 3 1 



Thesis 

M8635 

c.l 



202572 

Mouras 

Effects of atmo- 
spheric refraction on 
U. S. ground warfare. 



